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some Comments on Current Aviation Topics 


The Twentieth Wright Brothers Lecture 


SIR ARNOLD HALL* 
Hawker Siddeley Group Limited 


SUMMARY 


The lecture falls into two parts. The first contains comments 
on administrative matters affecting the research and develop 
ment process; it discusses particularly the time taken in the 
process, the place of scientific criticism, the weapon system con 
cept, technical education, and the relations between industry 
and government. 

The second part of the lecture contains comments on some 
current technical issues particularly affecting civil aviation 
The possibilities of long-range supersonic flight, the problems of 
structural fatigue, and the suppression of noise are among the 


topics discussed 


INTRODUCTION 


M* I THANK You for inviting me to deliver this lec 
ture on the anniversary of the first airplane flight. 
Paying one’s tribute to the memory of the Wright 
brothers, and to your Institute, is not an easy task in 
times like these when considerations of security must 


limit one’s discussion in a public lecture of advanced 
technical matters. That I, like others, have had to 
work within these limitations you will readily accept, 
and I make no further apology for suitably restricting 
ny comments. 

What is contained in this lecture will not be new to 
many of you, but I hope that it will not always be the 
same part of the audience that feels itself on too familiar 
ground. Since aviation isa matter of great public inter- 
est and since it disburses considerable public funds, it is 
not surprising that it should be the subject of inquiry 
from time to time, and I have on many occasions been 
called to answer questions in such circumstances. I 
have found that processes which are well known to engi- 
neers are often by no means fully understood by the 
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Formerly, 


people who become involved in these matters, and 
Part I of this lecture embodies statements and opinions 
on points I have found to be obscure to intelligent but 
nonprofessional inquirers. Part II is an attempt to set 
down some of the technical factors which will control 
aviation development in the coming years; it is by no 
means a complete statement for the reasons I have 
mentioned above. 

The opinions expressed are my own and do not 
necessarily reflect the views of the official establish 
ment with which I was, until recently, associated or the 
group of companies with which I am now associated. 


I) RESEARCH AND DEVELOPMENT 


THE RESEARCH AND DEVELOPMENT PROCESS 


Technological progress comes about in both an un 
organized and an organized manner. The organized 
approach is largely a twentieth century invention, and 
aviation has had a lot to do with introducing it. The 
steam engine was not invented because a government 
department set out to have it invented—nor was the 
first airplane—but most of the rest were the products of 
an organized approach to an extent which has increased 
as time has gone on. In organized research and devel 
opment, progress is sought as a matter of policy either 
by the state or by large private undertakings, and they 
bring together the people and the equipment thought 
necessary to achieve it to an extent far exceeding the 
financial capacity of any individual worker. 

This has brought about the large-scale employment 
in many countries of scientific and technical people, 
loosely described as scientists. These large numbers of 
scientists, except on rare occasions, do not contribute 
to the advancement of science in the classical meaning 
of the term; they add little to the broad sweep of our 
understanding of the universe in which we live. They 
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are essentially concerned with detail well within the 
outer boundaries of scientific knowledge, but detailed 
mastery of the fine structure is vital to efficient engi- 
neering exploitation. The broad science on which aero- 
dynamics is based was to a large extent complete be 
such men as 
But 


the application of the established principles of under 


fore the Wrights ‘..v the first airplane 


Newton, Reynolds, and Mach had seen to that. 


standing presents problems of the greatest detailed 
difficulty, and it is with these problems that the mass of 
modern scientific workers are concerned. 

Advance in technology comes in several ways, but 
some specific trends are recognizable. There are radi- 
cal departures from the situation that exists; typical 
of these are the change from wood to metal in construc 
tion, the move from the 350-m.p.h. regime of aviation 
during World War II to the 600-m.p.h. regime of the 
present high-subsonie military and civil aircraft, and 
the jump now taking place from the subsonic to the 
supersonic regime. These jumps produce a new 
position which is then usually the subject of develop 
ment. This produces very effective, if not radical, ad- 
vances involving the exploitation by detailed change of 
a device already achieved. This process is often taken 
further by continuity of development, which carries on 
the process of simple development by the creation of a 
new device essentially similar to that from which it was 
derived, and based on the same knowledge, but freed 
from the limitations —perhaps of size or strength —which 
brought the exploitation of the original to an end. 
Development by a series of short steps on a single basis 
of knowledge will rarely produce a radical change, but 
it can pave the way, and if previous development has 
not been taken far enough, the achievement of a radical 
change is made more difficult since the gap to be bridged 
is so much the greater. 

There are four phases in the attainment of a radical 
departure, which then forms the basis for further de- 
velopment. The first is that of undirected research in 
which detailed knowledge over a wide front is accumu- 
lated. 
that a basis of knowledge sufficient for building a useful 


The next phase starts when it becomes evident 


engineering device is in sight, and at this stage the re 
search becomes more pointed in the direction of the 
project visualized. The third stage starts when a de- 
cision is taken to attempt the project; then the normal 
mechanism of engineering design comes into play, re 
search becomes more pointed still in those directions 
where the knowledge upon which the pro‘ect must be 
founded appears to be weak, and detailed technical 
work specific to the particular geometry of the project 
There is then the fourth stage of making the 
Usually there are a 


starts. 


project work in every detail. 
variety of snags to overcome, some of which may well 
involve further detailed scientific work. Meanwhile, 
arrangements for production have been in hand; this 
process also involves very considerable research since 
the manufacture of new and advanced devices always 
poses new problems in workshops. So what might be 


described as the ‘Mark I” of a series is produced. 


FICAL 
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In this process the choice of time for the initiation , 
the engineering phase is crucial. If the project is high] 
commercial so that the main consideration is that it wil 
giving ever 


sell successfully and perform in a way 


satisfaction to the customer, the engineering work mug § 


not be started until the basis of detailed knowledg 


upon which its success depends, is very complete—s 


complete, in fact, that there are no doubts of an evey | 


mildly serious nature. Not many aeronautical proj 
ects have been of this type, though there have bee 
some. If, however, there are military matters involved 
then the main issue is that there shall be a counter t 
any threat the potential enemy may mount, at or be 
fore the date on which he mounts it. If the political 


situation is such that this matters a great deal, the ele 


ment of risk accepted in the way of incompleteness of | 


detailed knowledge upon which to base the project 
must be greater; technical advantage will go to thos 
who balance risk and time most successfully. If com 
plete success in the sense that it meets everything every- 
body had hoped for is demanded from the Mark I of 
new military project, it will almost certainly be too lat 
to counter the threat it was intended to meet because 
to ensure success, it will be started at a later stage in the 
build-up of background knowledge. In highly compet 
itive aviation development of the kind we have see 
in recent years, technical risks must be run; the alter 
native is to accept the fact that the potential enemy 
taking a greater technical risk, may get through and 
secure a dominant lead. 

What is learned during the project stage is not s 
much new basic knowledge but how to exploit the orig 
inal basis of knowledge in terms of what is often called 
“hardware,” and much of this learning is specific to the 
particular geometry involved. During the time taken 
in bringing the project to fruition, the basis of knowl 
edge upon which it is founded rarely changes much 
Occasionally there is a radical intrusion—such as when 
the concept of sweepback as an antidote to subsonic 
drag rise made its appearance —and if this happens, the 
project should be reassessed and perhaps dropped. But 
otherwise it is best to beware of ‘‘taking hands off the 
wheel.”’ 

Generally, by no means all that there is to learn can 
be applied in the Mark I stage, and Mark I projects 
usually do rather less well than their designers had 
hoped. They form, however, the essential foundation 
for development; it is usually the Mark IT design, an 
amended version of the first, which brings the exploita 


tion of the scientific foundation within full view. Ther 


are at least two reasons for this: (1) No designer 1s s0 


gifted as to be able to see on the drawing board ever) 


possibility of detailed improvement in his provect 
things are very different when a Mark I has been mad 
and experience gained in its use. (2) What is thought 
to be most important often changes with time, and a 
different and more advantageous balance of the give 


and take between the various performance characteris 


tics becomes acceptable when it is seen against the 
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SOME COMMENTS ON Crt 


background of the achieved performance of the Mark 
| machine 

The long course of development and the gains it 
brings about, though well known to engineers, are less 
appreciated by those without direct experience of it. 
Fig. 1 illustrates a typical example of the development 
of an airplane, showing what has evolved from the orig- 
inal Avro Manchester, which was designed in outline 
in 1936. The plan forms of the various aircraft in the 
series are shown, and their common derivative is evi- 
dent. The curve shows how the all-up weight in- 
creased; it more than doubled while the basic structure 
Each of the aircraft that evolved 
The 


Manchester itself did not go into service in numbers for 


remained the same. 
has given great service, and many are still in use. 


a variety of reasons, but it was the essential basis for the 
Lancaster, which was one of the most successful bomb- 
ers of all time. It is unlikely that those who placed 
the order for the original Manchester foresaw what 
would result from it, but they saw far enough to start 
the process off. Had they not done so, it is unlikely 
that the rest could have followed because it is much 
more difficult to reach a particular point in the evolu- 
tion directly. When great issues hang on technical 
superiority, the maintenance of a reasonably steady 
continuity of development is essential to success. 

The great increase in weight lifted by the same basic 
structure, extended and altered in detail, demon- 
strates the point that what is found to be possible with 
detailed development is usually far more than could 
reasonably be anticipated at the beginning. This is be- 
cause, in any engineering device, the ultimate limita- 
tions are not usually reached by all parts at once, and 
relatively detailed changes can greatly increase the 
capabilities of the device. Furthermore, it is almost 
always the case that some parts are found in practice 
to be capable of far more than it would have been 
prudent to assume at the beginning. Knowledge that 
this is so presents a problem of great subtlety in design. 
To what extent in the initial stages should provision be 
made—and how —for long development? Though this 
isa problem in all engineering, it is a particularly diffi- 
cult one in aviation since the margins between success 
and failure—in terms, for example, of weight or drag 
are often quite narrow. It is on this point that many 
scientific assessments of projects fall down, for it is 
not necessarily the design that is thought to be best on 
immediate criteria that will prove to be so as it proceeds 
on its evolution. Design is, and is likely to remain, an 
art greatly aided by science but by no means wholly 
contained in it, and this is one of several reasons why 
scientific assessment is not the whole basis for decision. 

Fig. 2 shows how the thrust obtainable from a typical 
jet engine has increased with time and again shows what 
can be achieved once a device is available for develop- 
ment. The long-term influence of continuity of de- 
velopment is illustrated in Fig. 3, which shows how the 
Weight of engines per unit power developed has fallen 
with the passage of time. The much more rapid ad- 
vance of the jet engine compared with the other types 
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reflects the fact that its development has been contin 
uous and that the resources put behind the task have 
been great. 

It is interesting to consider whether the scope for de 
velopment in future aircraft and engines will be as great 
asin the past. At first sight it appears that this might 
not be so, but I think that this is always the case when 
a radical advance is started. What is true is that, as 
design methods become more exact, the full potential is 
more readily achieved at the first attempt. But it is 
generally the case that, when a radical advance is 
started, new materials must be used and new construc 
tional methods adopted so that the accuracy attained in 
established design arrangements is less transferable 
than might at first be supposed; as more is learned 
about the new methods, scope for development again 
appears. 

THE TIME TAKEN IN THE RESEARCH 
AND DEVELOPMENT PROCESS 


The time taken to bring a radically new device into 
being is the sum of these elements: the time for re 
search to reveal the possibility, the time taken to make 
the decision to attempt to realize the possibility, the 
time in learning how to apply what is known in the 
broad to the detail of the particular project—to find 
out how to make it a practical article and to embody 
this knowledge in drawings from which it can be made 
the time it takes to build the prototype and make it 
work, and the time to produce in numbers. By com 
mon consent, the total time is getting longer and longer 
in aviation and, under highly competitive and unsettled 
international conditions, is becoming too long. I have 
no easy answer to this but would like to draw attention 
to one or two points which arise from considering it. 


Applicable Research 


The broad search for detailed knowledge, from which 
the possibility of a radical engineering change emerges, 
is the essential element from which all else springs. 
There can be little doubt that this applicable research 
is as good an investment as modern society can make. 
It is usually the case that it proves applicable over a 
wider field than is ever foreseen; aviation research has 
been responsible for advances in the properties of ma- 
terials to an extent far greater than any other indus 
try, and these improvements have benefited engineer 
ing at large. Whereas the application of knowledge, 
when acquired, is a matter of choice and the exercise of 
that choice depends on all sorts of factors such as, in 
aviation, the military threat, the state of the exchequer, 
and the policy of the political party in power, it is al 
most certain that the chance to make the choice will 
not be available unless an adequate carpet of knowledge 
is continuously created. But this does not mean, as 
some seem to argue, that undirected research should 
be free of all restraint. Much useless so-called scientific 
work is done under cover of the idea that anything called 


scientific research isa goodidea. There is no golden rule 
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which enables the right choice of research to be made. 
It is an essential qualification for those engaged in the 
direction and administration of research that they 
should have sufficient awareness of the kind of engineer- 
ing project that it might eventually become important 
to attain and sufficient ability to foresee applicability of 
the right kind to be able to make the right choice of 
work. It is a long and difficult process to create really 
good research laboratories, but they can deteriorate 
rapidly. Furthermore, they can run down a long way 
without it being apparent to those who are not skilled 
in research and its administration; the reports still 
come out, the paint on the door is just as fresh, but the 
knowledge being acquired is not of the right quality or 
applicability. It is essential to relieve these places of 
too much attention from people who know little of 
science and research and who therefore confuse the sub- 


stance with the shadow, thereby diverting the atten 


AERONAUTIC 


AL SCIENCES MARCH, 1957 


CIVIL TRANSPORTS [ [ [> 
(STILL IN SERVICE) 





> 
tuoor | 
MK.2 ¥ 
f \ } 
fT ~ 
i 
> f cy 
L , a a |__1—t 
~ — “ T + ——— 
> | p 
p> > 
| 
ef voor} 
MK ‘ MK4 
T-\ 
} 
Pp &> 
ADAPTATION iW | > 
OF BOMBER __ <5 + 
FOR PASSENGERS l | ad 
‘a 
h P 
ANCASTRIAN | 
| 
t » 
ji _ 
17 <=> 
+ > 
[> 
f 
INCOLN || 
1956) 
TFIMIAIM] J] JTATS]O[N]O [J] F [MJA]M] J] J] ATS]O[N]O | J|F [M[A[M] J] JTA]S[OIN]O] J]F]MTA]M] J] J] A[S]O[N]O) 
1943 1944 1945 1946 


tion of those who can detect the difference from doing 
so. 

The fact that aviation has such wide national im 
plications has led inevitably to the creation of large re 
search laboratories operated under the auspices ol 
governments. There have been many experiments 
made to try to find the best way of operating such 
places; my impression is that systems under which the 
establishments are shielded from the direct influence 0! 
a government department, through the medium of 4 
professionally knowledgeable governing committee re 
porting directly to a high political level, have much t 
be said for them. Though the test of the system 1s 
whether it attracts the right men, it is also true that the 
But 


the applicability of the work done is so crucial an issu 


right men will make any tolerable system work. 


in assessing its value that it is essential to create 4 
position of trust and close contact with the industn 
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Fic. 1 (Continued). 


that makes use of the knowledge and the user who oper- 
ates the engineering product. There is only one basis 
on which this position can be created: It is the recogni- 
tion that the official establishment should not seek to 
discourage research enterprise in industry in order to 
keep the work for itself and that it should not attempt 
to do project work for which the industry’s facilities are 
much more suited. There is a continuous and perhaps 
natural tendency on the part of the staff of some of these 
establishments to do both these things. The latter is 
particularly damaging; not only does it create a sus- 
picion that government sources are attempting to take 
work from industry but it aids a tendency in research 
establishments to evade the issue, which is to carry out 
good applicable research. This is not an easy thing to 
do if project work of a serious kind is allowed to intrude. 
Such project work creates a flurry of day-to-day execu 
tive work which clogs the minds of those who are not 


particularly suited by training to handle it, and it pro- 
vides an excuse for not facing the really hard task of 
breaking new ground. Those who work in a govern- 
ment-operated research unit should seek new applica 
ble knowledge, create the position of trust with indus- 
try and users which will greatly aid their assessment of 
applicability, and recognize that, if this is to be 
achieved, they must not try to do the other fellow’s job. 
They should also recognize that, if they do try to do 
industry’s job, it is likely that they will not do it well, 
but even if they do, it is certain to reduce their ability 
to create new knowledge of the right kind; it is by the 
applicable knowledge they produce that they will in the 
end be judged, and rightly so. A few odd per cent of 
gain should not be the target when means of securing 
more radical changes is the real aim; it is another ex 
cuse for avoiding the issue, but the few per cent—and 


more—will come in any case from normal development 
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Fic. 2. Increase of thrust from typical jet engine due to develop 
ment 


If the reason people turn to project work is because they 


without authority, the system will not work well 
decisions are delayed or avoided, and morale, enthy 
siasm, and determination, which are the elements essen 
tial to the speedy achievement of technically difficult 
tasks, are depressed. If I had to find a single statement 
that would encompass such difficulties as we have had ir 
the relationship between government departments, on 
the one hand, and the aviation industry and research 
establishments on the other, it would be “‘the interrela 
tion of authority and responsibility.’’ In a situation 
where each of the parties is essential to the other, ther. 


is room for further mutual study on how best to achieve 


cannot really do research well or because they prefer the proper balance of authority and responsibility on 


engineering application, then it is better that they both. 





should leave the research laboratory and 


where. 


Strengthtof Resources 


Few will dispute that the time taken in research or in 
project design relates directly to the strength of man- 
But this is not 


power and facilities that is deployed. 
to say that there is not ample room for individual bril 
liance and managerial skill. 


qualities are at work is, of course, that more seems to 


come from the resources used than might at first be ex- 


pected; one has only to examine the facts to be aware 
that there is a factor of two or three to one in the effec- 


tive technical productivity of research and design units, 
and the methods which produce this “gain from no- 
where” are of far greater value and importance than 
most other things in this process. 


resources” argument so overburdens the minds of some 


who have to decide, that they are apt—and it has hap- 


pened in England—to ‘“‘write off’ in their minds, and 
actively to discourage, those places that produce re- 
sults out of proportion to their apparent strength on 
the grounds that they are too small to be able to handle 
There is a lot to be said for ‘“‘backing the 
stable.’ Much harm is done by too blunt a use of the 
“number of heads’’ and the “technical planning”’ argu- 
ments. The fact is that individual flair and method can 
make a very big difference indeed, and it is of great im- 


the job. 


portance to recognize and encourage these qualities. 
But with a given level of skill, time is governed by 
strength of resources. 

A substantial part of aviation research and develop- 
ment is paid for by the state in support of military pro- 
grams. In the West the economic system is based on 
private enterprise, and aircraft are supplied in the main 
by privately owned companies, with a government de- 
partment having a good deal to say on what they do. 
This is an excellent arrangement—the best—but, un- 
less both parties are wise, it is possible to ‘‘fall between 


two stools.’ It seems to me that the point that is 


go else- 


The evidence that these 


The “strength of 


I would like to say a word about the British position 
on resources. The cost of aeronautical equipment, ex- 
cluding effects of monetary inflation, rose about ten 
times between the W orld 
War II and the high-subsonic and transonic eras after 
it. The British nation was not in a position at the 


low-subsonic era_ before 


end of the war to finance immediately all the new 
resources needed for modern aviation; constructional 
effort had to go into the building of houses and indus- 
trial plants to replace those destroyed by bombing 
Whether, in the circumstances, more might have been 
done is a political question which is not for me to dis 
cuss here; the view was that a large aeronautical re 
equipment program could not be started at the time 
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SOME COMMENTS ON C1 


the conditions left by the war. The present 


hecause O1 
of British aircraft have, therefore, had to be 


generation 
| without as much benefit as we would have 


levelope: 
liked from 
.s being such an immense help in achieving performance 
The critics 


the technical resources which are recognized 


nd in reducing the time of development. 
{ the British industry, both at home and abroad- 
ind there have been the United States 

have perhaps forgotten that, notwithstanding these 


some in 


iificulties, the present generation of British aircraft 
nd engines are as good as, or better than, their compet 
itors to the extent that 
buving them at the rate of 300 million dollars worth 


the world outside Britain is 
. year, and during the period in question the gas-tur- 
bine-propelled civil air liner, in both the jet and pro 
peller forms, has been added to the already long list of 
engineering firsts which came from Britain. That there 
is room for criticism of the position I would not deny; | 
lo not know how people can attempt and achieve diffi- 
cult tasks in difficult circumstances without being open 
to comment. The more curious situation is the one 
in which everything appears to be so good that nobody 
wishes to suggest that it might be improved. But the 
task of the genuine critic in aviation is not easy be 
cause the subject is complex. I cite an example 
the breaking of the world’s speed record calls for 
sunilarly, military and civil 


later 
much more than speed; 
aircraft are not adequately assessed on a single per- 
formance parameter. Aeronautical resources in Britian 
have now improved considerably; in the last 5 or 
6 years it has been possible to make a large invest- 
ment in equipment from both national and private 
funds, and this should have a significant effect on Brit- 
ish development in the next decade. 

One of the big factors in resources is technical man 
power. This has been talked about a good deal in re- 
cent times, and I think that the heart of what has been 
said is contained in Fig. 4. Figures such as these can 
be misleading because the definition of a graduate 
varies from country to country. The figures for Great 
Britain are not comparable with the others for this 
reason: In Great Britain substantially more people 
come into technology through nonuniversity channels 
than through the universities, and many of these people 
would count as graduates in the United States or the 
U.S.S.R. Nevertheless, too few able young men and 
women come into science and technology in Britain and 
probably in the United States. Undoubtedly one of 
the factors at work, anyway in Britain, is the feeling 
that science and technology are slightly odd. My own 
view is that there is no particular virtue in being a 
scientist, but nor is there a virtue in not being one; it 
ison the latter point that there is controversy, for some 
people still seem to think that there is. Such views 
should now disappear, and science and technology 
should be recognized as neither odd nor particularly 
Virtuous but as essential elements of modern society, 
capable of great influence and service. Centers of 
learning would act wisely by embracing these fields 
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fully rather than, as has regrettably been the case, mak 
ing them a matter for argument instead of action 

It is often interesting, and sometimes moderating, to 
look back on what has happened before. On doing this 
in the matter of scientific education and training I find 
that, over 100 years ago in England, there was a great 
outcry that there would be a disastrous shortage of 
scientists, and the terms used were very similar to those 
employed in the current argument on the matter; the 
present position is not new! But it was this outery of 
more than a century ago that made England a consider 
able pioneer in scientific and technological education, 
and this effort served her well. Now many steps are 
again being taken in Britain, and one of the most inter- 
esting is the creation of a fund to which all British in- 
dustry is invited to contribute (and has already done so 
to the tune of many millions of dollars) for the purpose 
of creating, in the best schools in England, scientific 
laboratories and libraries worthy of the place that sei 
ence should have in school education. 

There is, I think, one matter which comes out of an 
examination of the Russian position in aviation that 
bears on the subject of strength of resources and on the 
subject of the next section of this paper. Fig. 5 con 
tains outline drawings of a variety of Russian aircraft 
Until the point marked “‘first break,’ Russia has sup 
plied all her military and civil needs on one basic aero 
nautical formula—sweepback and boundary-layer fences 
and in addition there has been a considerable con 
sistency in tail-plane arrangements. (There are several 
other Russian types, in addition to those shown, based 
on the same formula.) This means that the basic data 
on the configuration and experience from any aircraft in 
flight was readily transferable from one to the other 
The economy that this represents in research effort and 


in detailed technical work is considerable and goes some 





16S TOURNAL OF THE ABRRONAUVUTICAL SCIENCES MARCH, 14 








—= — rh 
q& — ~ —_— 
; 2 — } cat 
S\ —O— YN * yy 9 
— , g ‘6 = 
Sorina <p } SO _—-_> —— > > F 
L = bi / : = ‘ Pony 
+ fs lad a SLa 
hong Z  : a fay ees 
> i Z — LZ = —_—== = 
« | SE ~ 
— hee ee col 
/ be 
Scit 
MIG.17. (BAY FIGHTER) SUPER FARMER (bdAY FIGHTER) SUPER PLASHLIGHT. (ALL WEATHER FIGHTER AND atracy m0 
pri 
Be 
ing 
) bee 
_ X \ | ire 
\ \ 
aa — — \ \, ™ oo Ly ——--9 SF ¢ = ane cer 
\ \ =O a ee oe 
‘ — \ wh 
id a > 
\ ————————— — — — — —-_-- 2 
=——____—___ > -B> — S37” = a om 
Leis! / 
L / L : _ a aS fe 
d ae = 3 fi y } ma 
J a fone 
<= = ly, 7 der 
FF, / —s ‘ 
lf / fi / 7 oS _——" Scie 
iy LZ ad 
ap] 
t10) 
BADGER. (BOMBER: TUIO4 CIVIL IS SIMILAR) BISON. (BOMBER: SIMILAR CIVIL UNDER CONSTRUCTION ) BEAR (BOMBER: SIMILAR CIVIL UNDER CONSTRUCTION ens 
} in 1 
ing 
we 
his 
eo 
al res 








RIGHT: FIRST BREAK 











me 
DELTA FIGHTER ' “ihe 
tut 
Sci 
. . . . . Sar 
Fic. 5. Russian aircraft showing use of sweepback and boundary-layer fence formula. (Reproduced with permission from Flight 
ess 
way to account for the rapid growth of Russian air Decision and Scientific Criticism ‘i. 
power. When the first break was made, it was pre- rae 
: Decision is important everywhere and not less so it pr 
sumably because the formula was thought to be ex- Ae : : as 
<o , ae a aviation upon which so much hangs. Decision mus un 
hausted in that the aerodynamic demands arising from ; ; ‘sae a 
: : ‘ : : take account of a good balancing of time and technical {u 
increasing supersonic Mach Number had stretched the , : : Ay f 
gee ‘ risk; must recognize that a good deal of technologica or 
structural problem far enough to justify turning to the See : ‘ ; : of 
; . ¢ ; information vital to the success of a project cannot be 
stiffer delta plan form, but the tail-plane arrangement é' é : . i 
obtained until the project is started because, before | © 


was not changed, and there was thus some continuity of ye, ae aay ; Sf 
: j é ; that, the need for it is insufficiently pointed; and must 
experience brought into the new formula. 


att “ rae ‘ have regard for the vital part played by development be 
ro secure the benefits that arise in this way, certain as as ' 7 
‘ ; ; and the need to provide a basis from which it can begin 101 
penalties must, no doubt, be incurred. <A great burden on 
is placed on the machinery for making the decision be- Scientific criticism has become a considerable factor in th 
cause, if it is wrong, the consequences are severe. But decision in recent times, and quite rightly so, but there om 
this burden is less than might be supposed to this ex- are dangers that must be recognized. First, the sciet- | 4, 
tent: Provided sufficient work has been done on the tist, as he is at present educated, is endowed with a - 
fundamentals of the arrangement chosen to eliminate technique of analysis which is assumed to be logical th 
any possibility of a radical difficulty, then, given deter- but he is not necessarily endowed with initiative, in vi 
mination, the arrangement can be made to work. It spiration, and courage. There are some who argue - 
may not be ideal in some circumstances, but departure that a scientific education of the kind now common has on 
from ideal is presumably thought to be justified by the a positively adverse effect on these qualities; I doubt ce 


gains in other directions. if this is so, but it does not do much to encourage them 
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[fhe point of danger is this: Most scientists are edu- 


cated in a Way that gives them considerable faith in the 
logic of their discipline but little idea of its limitations 
ind how to overcome them. Belief in its accuracy is 
right provided the basis from which the deduction 
starts is full and complete, but when decision must be 
reached in a situation of any difficulty, it is rare that a 
complete basis for subsequent scientific deduction can 
be stated. Because there are always imponderables, 
scientific assessment can be misleading unless it is 
moderated by a reasonably full understanding of the 
processes Of application and what can occur in them. 
Most courses of study in science, and many in engineer- 
ing, leave this matter quite out of account. It has 
been the pose of many academic men whose disciplines 
are in the humanities that their students are not con 
cerned with any vocational aspect of their studies—yet 
what greater contribution can an educated mind make 
to life than by application through a vocation and how 
many do not, in fact, apply it this way? Equally, 
many teachers of science take the view that their stu- 
dents are concerned only with the purity of a line of 
scientific thought, untrammelled by considerations of 
application. It would be better if the study of a sec- 
tion of the history of science and of the history of some 
engineering developments were regarded as essential 
in the preparation of any student of science or engineer- 
ing and if some study of the methods of application 
were also seriously made. Anyone who studies the 
history of science will soon be convinced that its prog- 
ress, far from being highly logical, has been just the re- 
verse. The reason that it appears to be logical is only 
because it is taught that way; it is convenient to re- 
member only the work that was later seen to fit to the 
logical pattern. But an examination of the fine struc- 
ture of the history of science very quickly shows that 
scientific evaluation of the next move does not neces- 
sarily get the only, or even the ultimately right, answer. 

Scientific evaluation is an essential part of the proc- 
ess leading to decision, but scientific people concerned 
in it too often take less account of the other parts of the 
process than they could or should; they are sometimes 
underequipped in their feeling for the ultimate conse- 
quences of the influences they have felt themselves 
forced to neglect. When progress in the acquisition 
of knowledge is rapid, it is easy to advocate that de- 
cision be delayed because something new is about to be 
found or that a project be stopped because something 
It is also easy to state a case 
These 


better is now possible. 
for a large number of variants on any theme. 
are serious difficulties which undoubtedly contribute to 
the length of time taken to bring scientific possibilities 
to practical use. There is no escaping the fact that, un- 
less enormous resources are to be deployed, this position 
must be firmly controlled by decision that has regard for 
the great investment in technical effort that each indi- 
vidual project represents and for the necessity, in times 
when technical superiority matters, for reasonably 
continuous development of hardware as well as con- 


cepts. If scientists use their formidable techniques of 
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analysis to do no more than delay decision by argu- 
ments on marginal points and fail to study and take 
account of the other issues or if they advocate indis 
criminately too many variants on the theme, they do 


less service than they could and should. 


The System Concept 


I think that one of the factors that has had an adverse 
effect on the time taken to achieve aeronautical de 
velopments is the weapon system concept. Of course 
no reasonable person can think that a weapon should be 
other than a system designed to produce the end result 
desired; this is just as true in civil applications. But, 
under the cloak of a term such as ‘“‘weapon system, 
project planning can go adrift, and I think has done so 
in more than one case. There are many elements in 
for example, in a fighter, the air frame, 


any weapon 
and they must 


the engine, the radar, the armament 
operate together harmoniously. But if a high level of 
technical risk is taken in the choice of each of these ele 
ments, then the burden of the technical work needed to 
bring them all along to the same end point is very great 
indeed; if all the babies were born on one day there 
would be a shortage of midwives. Yet this is what has 
sometimes been advocated as necessary to a weapon 
The practical difficulties soon become almost 
the wing wants to 


system. 
overwhelming. The designer of 
know the aerodynamics and loads of the weapons he 
must attach, but these are not yet fully designed; later 
the weapon must be changed because of development 
snags, but this means changes in the wing when the jigs 
are completed and metal committed to them. This 
situation is multiplied many times over in a complicated 
system. If a new project is to come through reason 
ably quickly, it should be so conceived initially that the 
various elements that make it up are not all at the same 
state in their evolution and, in particular, are not all 
at the start of their evolution. By way of example, 
consider a fighter and suppose that the air frame repre 
sents a radical advance. I think the best progress, in 
time and in quality, will come from starting the new 
air frame around engines, radar, and armament, or at 
any rate two of them, which have reached or are ap 
proaching Mark I standard and which are now pro 
ceeding in their evolution by development. When 
the air frame reaches the Mark I phase, its de 
velopment should start and, during this process, radical 
departures made in one or more of the other items 
This means that at no time are all the elements that 
make up the system at their most advanced technical 
stage. By this staggering of the evolutionary stages 
of the elements in the system, I believe that a great 
deal of time can be saved. Those who argue that this 
will produce a less good article than could be had on 
the basis of the available knowledge take too little ac 
count of the time saved, with the consequent better 
matching to the military threat. They also underesti 
mate what can be gained in the development time of 
each element when there is reasonable certainty about 
the state of the other parts with which it must mate. 
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Of course, one cannot generalize too much; each case 
is one for consideration on its merits, but, unless the 
points mentioned above are kept in mind, the ideal 
weapon system can be so far off in time that, when it 


arrives, the crown of idealism has gone elsewhere. 


Production 

The production phase of aircraft is one of consider 
able difficulty in the following respect: If machines 
with good interchangeability of assemblies 
with reliability of 


leading to 


good maintainability—and per- 
formance are to be supplied, complex tooling and _ jig- 
ging is highly desirable. However, this impedes the 
rapid attainment of the best possible technical position 
because changes which would be advantageous cannot 
be introduced without serious impact on the tooling, 
unless this is not constructed until the prototypes have 
done full tests; in that case, the time for supply is ex- 
tended most seriously. Having got the tooling there is 
inevitably a tendency for its existence 
to hold back develop- 


in terms of the 
money and time it has cost 
ment. There cannot be any doubt that, to produce a 
sound product, extensive jigging and tooling is neces- 
sary, but I am inclined to wonder whether it is some- 
times taken too far. At a time when technical superi- 
ority may well be all important, the production engi- 
neer is presented with an opportunity to make what can 
be a vital contribution by achieving the maximum 
flexibility in his position. The methods adopted in the 
West and in the U.S.S.R. are in contrast; the Russians 
appear to use far less jigging and tooling than the 
West, and they no doubt pay some penalty in the 
maintainability and performance of the product and 
the manpower used to make it, but they probably gain 
in time and in ability to make change to an extent that 
justifies this position in their eyes. Research and de- 
velopment is as necessary in production methods as it is 
on the technical side of the aircraft; the creation of a 
position as receptive as possible to the most up-to-date 
technical requirement must take a position of impor- 
tance alongside the achievement of rapid, accurate, 
and cheap output. 


(II) THE CURRENT TECHNICAL POSITION 


In this part of the paper, some matters that affect 
the development of aviation for civil purposes are dis- 
cussed and serve to some extent to illustrate points 
made in the previous part. Though a discussion of ad- 
vanced military aircraft might be thought to provide 
better illustration, the operational efficiency of such 
machines depends on factors which cannot be discussed 
in an open lecture. I will therefore say little on the fu- 
ture possibilities of such types, but where digression 
from the civil theme is reasonable, the discussion is not 
confined to passenger machines. 

Three technical matters from the many that affect 
civil aviation are of particular influence: the direct 
operating economy attainable, the integrity and useful 


life of the structure, and its noise nuisance. All three 
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are interrelated, and all three place limits on what it ; 
attempt. The 


ticularly on these three issues but again without resist. } 


wise to discussion 1s centered par 


ing digression when it seems justified by interest 


THE RANGE EQUATION 


In Fig. 6 is plotted the relationship between weight 


speed, engine performance, and range for aircraft whic! 


spend a large part of their flight cruising in the strat 
sphere at constant lift/drag ratio. (These data deriy 
from the energy equation; the connection is outline 
in the Appendix.) The range needed Xx, the speed pr 
posed |’, and the specific fuel consumption c fix Re 

on the left-hand plot. This value is then related t j 
the line marked with the lift/drag ratio expected 


horizontal transfer to the right-hand plot then pr 
duces the relationship between the all-up weight 
take-off IV’), the basic weight II’,, and the weight of th 
pay load together with items essential to its carriag 
We. (The basic weight Il’, is the sum of the weight : 
the structure, the power plant, and the systems; se 
the Appendix for details of this item and of IV. 

The main feature of this type of presentation is the 
rapidly varying slope of the curves on the right. This 


slope is an indication of the sensitivity of the design 
where it is high, the aircraft is ‘“‘easy”’ to achieve sine 


— 


failure to produce the lift /drag ratio or the basic weight 
hoped for has relatively mild results on the pay load 
Where the slope is low, the design is very sensitive ti 
small changes in the aerodynamic or weight positions 
and the airplane is ‘‘difficult”’ in the sense that the con 
sequences of relatively small technical failures are great 
The regions of “‘easy,’’ “possible,” and ‘“‘hard’’ desig 
are indicated on the plot; they are not, of course, hard 
and fast divisions. 

The economics of the airplane can be derived roughh 


from the weight ratios emerging from Fig. 6; an aj 


proximate equation for this relationship is mentioned 
in the Appendix. “difficult 
category, its economics move adversely since the frac 


As the airplane enters the 
tional pay load is becoming small. The next sections 
contain a discussion of values to be expected for the 
parameters on the range equation which govern the 
position of the airplane on the plot of Fig. 6 and which 
therefore control the economics it can offer. 


THE TRENDS IN THE PARAMETERS IN 
THE RANGE EQUATION 


Lift/Drag Ratio 


At subsonic Mach Numbers approaching 0.5), 
lift/drag ratio of up to about 16 is achievable with 
conventional arrangements; thereafter there is a rapid 
fall with increasing Mach Number until, in the me- 
dium-supersonic regime (Mach Number 1.5 to 3), 
ratio of 5.5 to 6 appears likely for a variety of practica- 
illustrated im | 


ble airplane configurations. This is 
Fig. 7. (Subsonic values higher than 16 are, of course, 
obtainable with special arrangements, but this figure is 
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typical of good practice in conventional civil aircraft. 
Anticipating for the moment the discussion on trends 
in specific fuel consumption, a typical consumption 
figure for jet engines at high-subsonic Mach Number 
is about 1.0 and, at Mach Number 2.5, is likely to be 
bout 1.5. Using these figures for L,D and c, points 
for a high-subsonic aircraft of 5,000-mile range and a 
Mach Number 2.5 aircraft of 3,500-mile range are 
shown on Fig. 6. If, as is reasonable to assume, a 
basic weight of 45 per cent or better is achieved, both 
uircraft are at about the same position in the “‘not too 
difficult’’ class. The possibilities of medium-range 
operation at medium-supersonic speeds are clear. Un- 
less, however, a substantial improvement can be made 
in the lift drag ratio or in the specific fuel consumption, 
very long ranges at medium-supersonic speeds will not 
be easily achieved using the present type of chemical 
fuels; if the 3,500-mile cruise on which the supersonic 
point is based is changed to 5,000 miles for direct com- 
parison with the subsonic point, it is evident that a de- 
sign emerges which has a very low pay-load percentage 
ind which is extremely sensitive to small changes in 
L D and basic-weight ratio. The situation would re- 
vert to the ease of the 3,500-mile-range case if L, D were 
to rise to about S.5 or if the specific consumption were to 
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drop to 0.9—the latter being very improbable without 
departure from present-type chemical fuels. (In con 
sidering these range figures with civil applications in 
mind, it should be remembered that a still-air range 
exceeding the stage length by about 60 per cent is 
needed to provide for weather and safety allowances. 
The drop in lift /drag ratio in the supersonic regime 1s 
due to the onset of wave drag. I would like to make a 
brief review of the advances made in recent times in 
controlling it. 
appearance at subsonic forward speeds, but consider 


A few years ago, wave drag made its 


able progress was made in delaying its onset by keeping 
the local effective flow over the lifting surfaces wholly 
subcritical (below the speed of sound) by the use of 
sweepback and thinness. When this is done, supersonic 
flow patterns, with consequent wave drag, arise at the 
junction between the wing and the fuselage and at the 
wing tips. Dr. Kuchemann at the Royal Aircraft 
Establishment pointed out several years ago that if the 
body is suitably ‘“‘waisted”’ it can be arranged that, in 
transonic flight, the pressure field due to the combina 
tion of wing and body approximates to that of the 
A sub 


sonic flow pattern is thus retained at the junction, and 


wing alone at points remote from the body. 


the onset of wave drag does not occur earlier there than 
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We = OPERATIONAL FIXED WEIGHT EQUIPMENT. 
Wa = BASIC WEIGHT = Ws + Wp + Wy. 

Ws = STRUCTURE WEIGHT. 

Wp= POWER PLANT WEIGHT. 

Wy = SYSTEMS WEIGHT. 


The range relation for cruise in the stratosphere. 
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on the main part of the wing. 
area rule, which originated concurrently in the United 
States, was a more general form of the same principle; 
this can be stated as follows: A reasonably smooth 
airplane flying at near sonic speed has the same drag as 
a body of revolution with the same cross-sectional area 
in planes normal to the flight direction. Linear super- 
sonic aerodynamic theory was shown to be consistent 
with the transonic area rule as the Mach Number ap- 
proached unity, and this led to the investigation of 
whether linear theory could be used to develop an 
analogous rule applicable to wave drag on supersonic 
airplanes. The supersonic area rule emerged, expressing 
the drag in terms of the cross-sectional areas in planes 
tangential to the characteristic Mach cone; this area 
rule, when applied to general shapes, gives only an ap- 
proximation to the result provided by the correct linear 
theory. These ideas have led to attempts to design air- 
craft having a cross-sectional area at each section con- 
forming to that of a body of revolution giving the mini- 
mum drag for a given thickness/length ratio, such as the 
Sears-Haack body. 

This remains a potentially fruitful field of research 
directed to answering at least three questions: 

(1) Is further progress possible in getting the best 
compromise between the shapes which emerge for sub- 
sonic, transonic, and supersonic speeds? This question 
is relevant to the performance of medium-supersonic 
long-range aircraft in the climb and through-the-speed- 
of sound phases; in such aircraft the installed thrust, 
which balances total drag in the medium-supersonic 
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Whitcomb’s transonic 


SCIENCES MARCH, 19: 
regime, inevitably exceeds by little the drag in the trap 
sonic condition through which the aircraft must fly, 


(2) How far in terms of forward supersonic spee 
can the onset of serious wave-making drag be delayed 
carrying on the process that has already made so big ay 
impression on transonic drag? In the low-supersonj 
regime, methods of the area rule type may have appli 
cation to this end. The area rule approach does not 
specify the cross-sectional shape of the fuselage but onh 
the area, and it has been usual to waist the fuselage wit} 
circular cross sections; the Kuchemann  approacl 
specifies an indentation in the plane of the wing. It 
would seem that an extension of both approaches woul 
specify not only the area but the cross-sectional shape 
Such a method has been suggested by Dr. W. F. Hil 
ton.* 


not only at the surface but also in the whole flow field 


In his approach, the streamlines are determine 
surrounding the profile. Hilton claims that for wings 
with subcritical flow the shape produced is substantially 
independent of design Mach Number and is thus effec 
tive at subsonic, transonic, and supersonic speeds. This 
is an interesting field for research which may produce 
in the low-supersonic regime, substantially better lift 
drag ratios than are now current. 


9 


(3) Are there configurations which can show a reduc 
tion of drag at medium-supersonic speeds, thus in- 
creasing the lift/drag ratio to a value greater than the 
I think 
that this may prove to be quite difficult in practice 


figure of roughly 6 which is now attainable? 


but there are theoretical shapes which offer possibilities 
and these should be investigated vigorously because the 
I think it is likely that 


wing and body shapes which are successful in this way 


reward for success is very high. 


will not have good low-speed aerodynamic qualities 
so that they may involve the use of special arrangements 
for take-off and landing. But my impression is that no 
easy way has yet emerged for obtaining medium-super 
sonic lift/drag ratios greatly exceeding 6, and the vir- 
tues of new shapes may be found more in their struc- 
ture weight than in their aerodynamics 
Summarizing this discussion on lift/drag ratio, | 
think it 
regime, one can attain something approaching 6 with 


comes to this: In the medium-supersonic 
practical configurations, and this is enough to make 
medium-supersonic flight at ranges of 3,500—4,000 miles 
possible and attractive for some purposes. We should 
look hard for higher L/D, but we will not find it easy to 
obtain; if we find it, it may involve configurations 
which require special provisions for landing and take- 
off (this is not, of course, without its own advantages 
New configurations that can be built at a low structure 
weight are also of great interest, even if their lift, drag 
ratios are not high. In the low-supersonic regime, we 
may well attain a lift/drag ratio substantially better 
than 6, of the 
rule” type may prove to be fruitful. 


and research “extensions to the area 


* Tests of a Fairing to Reduce the Drag of a Supersonic Swept 
Wing Root, Journal of the Aeronautical Sciences, Vol. 22, No. 3, 


pp. 178-178, 188, March, 1955 


-- 
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that makes long-range flight at medium-supersonic 


: : € . as When Mach Numbers near 2.5 are reached, the two 
speeds reasonably possible. For supersonic flight the 


factors are no longer in serious conflict; the propulsive 


efficiency has risen to a high value, and the materials of 


I think 


practice 
combination of air frame and turbojet power plant is at 


its best at Mach Number of 2.5 to 3.0, provided that 
the problems of kinetic heating which will be encoun- 


tered at these speeds do not cause an increase of the 


sibilities, 


ause the ; ; 
” which the turbine is constructed can be worked to the 


ely that ‘ é ; 
ee peak of their properties, with benefit to the thermal 


‘his way } ficiency and without so adverse an effect on the overall 
. , . . . elhciency < ) é ‘rse ¢ Tect o overs 
jualities structure weight sufficient to reverse this favorable 





gements trend. 


that no It is interesting to compare the compression and 

n-super- expansion cycles for jet engines operating at subsonic PROPULSIVE EFFICIENCY 
the vir- ind supersonic speeds; such a comparison is shown in 

ir struc- | Fig. 10. (This Figure is plotted logarithmically, so 60 fF 


that the elements due to intake and compressor, for 


ie example, multiply to give the overall ratio.) The 
ratio, 1 . or : - 
change between the two conditions is clear and ac- 30 ——— 
eTSOMIC ° . ‘ —— . , ° 
counts, of course, for the higher efficiencies achieved in 
the supersonic case. THERMAL EFFICIENCY 
The simple-jet engine is at an ideal state at Mach 


6 with 
9 make 
10 miles . ahs ay , 
Numbers near 2.5 in the sense that it is then making the 


should Se ; 
best use of its materials, at least as we know them now. 
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efficiency as is shown in Fig. 11 for a lower Mach 
Number. 

If this interaction between propulsive efficiency and 
thermal efficiency is to be avoided in the lower Mach 
Number ranges, the intercoupling which is implicit in 
the simple-jet engine must be broken in a way that per- 
mits greater independence of the propulsive efficiency 
from the thermal efliciency. This is done in the pro- 
peller-turbine engine, which therefore permits much 
higher propulsive efficiencies in the lower subsonic re- 
gime than does the simple-jet engine. The types of 
engine known as the by-pass and the ducted fan also 
effect such a decoupling; the extent of the decoupling is 
greatest in the ducted fan of high by-pass ratio. 1 
doubt whether there are theoretical differences in per 
formance between simple-jet engines and by-pass en 
gines of low by-pass ratio that exceed the variations in 
either type likely to be found owing to degree of excel 
lence of design. There is, however, a substantial possi- 
bility that a ducted fan working at a by-pass ratio of 
about 1.5 would give a better product of propulsive and 
thermal efficiency. The curve on Fig. 8 shows the ex- 
tent, it seems to me, to which this possibility could be 
realized. At high-subsonic Mach Numbers, the higher 
efficiency available from such an engine, as compared 
with a simple jet or a low by-pass ratio jet, is sufficient 
to give a significant improvement to long-range air- 
craft, after taking account of the higher weight of the 
ducted fan. This advantage disappears at early super- 
sonic Mach Numbers; the ducted fan is an engine which 
offers advantage only over a restricted range of Mach 
Number, but within this range is the speed that is of 
great interest to a most important section of civil avia 
tion; we will return to this point later. 


Basic Weight: The Influence of Configuration 


For Mach Numbers in the region of 2, the straight 
wing, the swept wing, and the delta wing are fairly 
equal contenders from the viewpoint of supersonic drag. 
Undoubtedly all three can be made to work. We have 
had interesting experience in England recently with the 
tailless delta configuration at supersonic speeds, par- 
ticularly in the Fairey F.2. which raised the world’s 
speed record from 822 to 1,132 m.p.h. This configura- 
tion has proved itself capable of providing good control 
and performance characteristics at supersonic speeds. 
Incidentally, the breaking of the speed record is a task 
that requires more than speed alone; the aircraft must 
be capable of very delicate control in pitch to meet the 
altitude requirements, and the fuel capacity must be 
such that the take-off and climb, the first run, the turn, 
the second run, and the return to base can be accom- 
plished without the possibility of fuel lack limiting the 
accelerating phase, particularly that before the second 
run. 

Above Mach Number 2, the case for the straight wing 
steadily mounts, and there is much to be said for this 
configuration for aircraft starting their lives with a 
performance in the region of Mach Number 2.5. 

There will, no doubt, be the usual marginal gains in 
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basic weight arising from development. The specif 


weights of engines will continue to fall, though probabh 


slowly. The design of aircraft structures for Mach 
Numbers exceeding 2 must take account of the effects of 


kinetic heating, and materials new to aviation such as 
steel and titanium alloy must therefore be used, but it is 
improbable that structure weights will be 
affected. 
pect any substantial fall in structure weight unless 


seriously 


However, it seems unlikely that we can ex 


configurations different from those which have become 
conventional are adopted. 

The factor that can and probably will bring radical 
changes in aircraft configuration is the use of engine 
generated lift to relieve the wings of any duty at take 
off and landing. The development of vertical take-off 
can be viewed in several ways: It may be regarded as 
desirable in its own right (this could certainly apply to 
some Classes of military aircraft); it may allow the us¢ 
of a configuration with superior aerodynamic cruising 
characteristics, no longer compromised by the need for 
good low-speed properties; or, again because the low- 
speed characteristics are no longer important, it may 
allow the use of a configuration which can be built with 
a decidely lower structure weight than conventional 
layouts. 

Engine lift can be generated by a rotor (as in the 
helicopter), by a ducted fan, or by a gas jet. These 
methods vary downwards in their propulsive efficiency 
and therefore upwards in their fuel consumption per 
unit hovering time. The rotor is suitable only for use 
on aircraft of relatively low forward speeds. Fans re 
quire a wing of fair plan area and aspect ratio to permit 
their installation, and they are therefore most suitable 
for aircraft of low to intermediate (say Mach 1) speed 
range. The jet engine in conventional form is most 
suited to a plan form such as the narrow delta, which 
can provide the depth needed for installation while 
having characteristics suited to supersonic speeds. | 
expect that even more radical integrations of the power 
plant and structure will appear. 

The narrow delta vertical-lift aircraft, first suggested 
by Dr. A. A. Griffiths of Rolls-Royce, Limited, is an 
interesting study in exposing the problems. Work in 
England has suggested that the drag of narrow delta 
configurations at supersonic speeds can be predicted 
with fair accuracy as the sum of the turbulent flat- 
plate skin friction plus a margin, the theoretical wing 
wave drag as predicted by slender-body theory, the 
drag of the fins alone, and base drag (if any). The 
reason why the friction component seems to be greater 
than the flat-plate value is not fully understood and ts 
an interesting subject for further research. On_ the 
present basis and assuming a structure weight that 
is suggested by fairly detailed investigation, it emerges 
that a jet lift narrow delta configuration cruising at 
Mach Number 2.5 would have a range rather less than 
that of a conventional aircraft of similar weight, but 
not so much so as to make the concept of no interest. 
The dominant factor in this assessment is not that any 
marked improvement in lift/drag ratio emerges but 
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Ce th th adic: ’ differe . igurati . eS ; 
le specific that the 1 udic lly different configuration produces i 5 
structure weight distinctly lower than conventional 


probably °5 ns ; ; aad 
‘or Mac} ' wrangements, thus offsetting the weight of the lifting 
ach ‘ ae we ; 
effects of engines and their fuel. There is a long way to gO be Paonia pict : 
n such ; tween an assessment such as this and the engineering ann ; "3 , 
as z : $ RB prop ARN 
|, but it j whievement embodying the principle. However, there (000 #1 oF n-red 
is suflicient promise to make the position interesting FOR DUCTS AN , “ 


seriously : 

e can ex- and to show that enough may be gained in structure 

st unten ) weight by departing from conventional layouts to make 

what is, at first sight, an uneconomical method of 
Comparison of turboprop and jet-flap aircraft for 

\ll aircraft have same range and pay load 


Fic. 12 


dium-range duty 


‘ become ; 
landing the vehicle a reasonable proposition. Even 
: eon more radical departures from conventional configura- 
g radical 
f engine tions can be conceived and may well show greater gains; : 
sien pared the direct operating costs of a jet-flap aircraft 
at take- 
take-off 


arded as 


certainly this is a field with room for imagination. But ‘ 
and a conventional propeller-turbine aircraft suitable 


~ 


hopes for city center operation for such aircraft should : i : 
. : : , for medium-range operations (at the same range and 
not be raised too high, for they will be very noisy - “s ws 


The runway 1s, and ts likely to remain, an excellent pay load) plotted as a function of the runway length 


specified. These results are at the mercy of the assump 


apply to 
tions made about weight penalties and the like, but 


method of getting into the air and returning to the 
ground. Arguments are mounted that runways are 
altogether too costly a luxury, but I think that there 


’ the use 
cruising “o 
ing variations to the extreme of optimism do not seem to 


need for me to alter the following general indications 1) City 


~ 


is little justification for such a view. The capital 


the low- ta 
i ; : : fe center operation by jet-flap aircraft is not attractive 
it may costs of base equipment for other transport systems : es, . 
2 re ae hi : since the runway length below which cost rises rapidly 
lilt with such as railways, which have a much more restricted ; ; sh ; 
is not short enough. (2) Though the jet-flap aircraft 


command of distance, 1s also high; the case against ; : 
is cheaper to operate than the conventional aircratt 


for field lengths below 2,000—2,500 ft., it is still sub 
stantially more costly than a conventional aircraft 
Therefore, 


entional 
aviation base costs can hardly long survive if the issue 
is put into perspective by considering the capital put 
d into railways a century ago. I think that the con- 
tinued use of long runways is economically justified. 


sin the 

7; dae operating from a field of 4,000 ft. or more. 

iciency : ; » . 

we. is j , : : ; it seems to me that, since 4,000 ft. is by no means a 

But this is not to say that improvement in runway per- apa : : - 
difficult length of runway to find, operators will not 


think the jet-flap type attractive. 


ion per 
for use formance may not be desirable for some classes of air- 
craft, and reduction—or at least a curtailing of the 


‘ans Te: However, I expect to see interesting results coming 
permit increase —in approach speeds will be helpful to safety ; / ~ S “4 
: a. from unconventional approaches to the design of both 
suitable in bad weather. ' a 
. ; P , : the jet and the structure, the aim being a higher lift 
) speed ) For marginal wuprovement in runway pertormance, i ‘ ‘ - sai ; 
. - ; ae : drag ratio, less structure weight, and high mobility in 
s most blowing over the flaps’ is probably the best answer ; ‘ ; ae ; ee 
: : the direction of application of the thrust, without 
which and, indeed, may be the only one (apart from conven- ; , ae 
: : . ; ; weight penalty —or, indeed, with a reduction of weight 
. while tional slots, slats, and flaps and reduction in wing P 
eds. | loading). The emission of the propulsive jet along 2 
aa : ; sis ; . APPLICATIONS TO CIVIL AVIATION 
power the trailing edge of the wing, an arrangement known as 


‘gested 


the jet flap, is a method suggested for improvement of a 
more radical kind but short of full vertical take-off. 
The principle has, I think, been put forward sepa- 


Long-range high-subsonic jet-propelled civil air 
liners (typified by the figures 480 knots speed, 5,000 
still-air nautical miles range, 25,000-30,000 Ibs. pay 


, is an 
ork in rately from what is now known as the jet flap, but in load) will show distinctly lower costs of operation than 
delta more radical forms, and I think such approaches may the aircraft they will replace; this is the result of the 
dicted produce remarkable results in the end. But the follow- larger annual earning capacity with which their speed 
t flat- | ing discussion is centered on the jet flap as at present endows them. At ranges much above 5,000 nautical 
| wing publicized, and a rough comparison follows between miles the simple-jet high-subsonic aircraft is moving 
y, the aircraft so equipped and conventional aircraft for the into the ‘“‘difficult’’ class on the range equation diagram; 
The same duty. The jet-flap aircraft pays some penalty in if substantially longer nonstop ranges are required, 
reater weight for the distribution of the propulsive jet, and, they can best be got by the use of a ducted-fan engine 
and is since the lift-—as well as the thrust--depends on engine without sacrifice of speed or, for even longer distances, 
n the : power, there is an adverse effect from the ‘‘engine cut”’ by the use of a propeller as the driving element, with a 
that case as compared with conventional aircraft. How- drop of speed of 70 or SO knots. But the 4,500—5,000- 
lerges ever, the wing of the jet-flap aircraft is smaller than mile range airplane is undoubtedly capable of satisfying 
ng at that of the conventional aircraft with the same air- a very substantial part of the need, and in this class 
than field requirement, and the gains in structure weight the jet-propelled machine is the right answer; | think 
_ but can make the jet-flap machine lighter if the wing load- it will also be the right answer for the medium to long- 
erest. ing of the conventional aircraft with which it is being range bracket. There are two particular ways in which 
t any compared is low enough——that is, if the runway length greater economy can be brought to this class of ma- 
Its size can be increased, and it can use a 


; but 


specified is below a certain value. In Fig. 12 are com- 


chine: 
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Fic. 13. Comparison of seat-mile costs 


ducted-fan engine. The pursuit of such developments 
depends on whether the high-subsonic regime is likely 
to remain a field of high commercial interest for a very 
long time. I think it undoubtedly is, and I will return 
to the subject after attempting to justify this state- 
ment.* 
Present-day 
direct operating cost on British standards of about 24 
In practice, there is 


long-range aircraft have a formula 
pence per long-ton nautical mile. 
an addition to the formula cost which varies with the 
route, but in round figures it is about 10 pence. The 
indirect cost is about 22 pence. For comparative pur- 
poses, I will use the formula direct operating cost as a 
measuring stick. The high-subsonic jet aircraft suit- 
able for nonstop trans-Atlantic operation will bring 
the formula direct cost down to about 11 pence per 
long-ton nautical mile. To achieve a cost within 10 
per cent of this, a low-supersonic airplane, suitable for 
the same flight, would have to achieve a lift ‘drag ratio 
of about 14. Even if understanding of how to control 
wave drag in the early supersonic regime extends 
further, it seems imprudent to expect anything better 
than, or even as much as, 14, having in mind that 16 is 
a very good figure to achieve in a conventional sub- 
sonic civil airplane. Therefore, long-range civil avia- 
tion in the low-supersonic regime will be more expensive 
than subsonic—and probably by a substantial margin. 
Looking next at Mach 


capable of nonstop-Atlantic operation is in the ‘‘very 


Number 2.5, an airplane 


difficult” class; it is extremely sensitive to small de- 
partures from ideal in its design. Taking figures for 
lift ‘drag ratio and specific fuel consumption that are 
reasonable on present knowledge, the direct operating 
cost per long-ton nautical mile cannot be less than 100 
pence. This position could improve as our knowledge 
becomes more exact, but it would be imprudent to 
suppose that such an ideal state could be reached 
quickly without substantial experience in this regime of 
flight. 

Although I have quoted figures in the previous para- 
graph for Mach Number 2.5, which is where the com- 
bination of air frame and engines is in sight of its best 


* Some hold that greater economy can also be had by the use of 
boundary-layer control to maintain laminar flow; I do not pro- 
pose to discuss this here, having done so in another paper (J. Inst 


of Transport, Brancker Memorial Lecture, 1955) 
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performance, I think that, in considering civil aviation 
it would be more prudent to think of Mach Numbers 
in the region of 1.8 or 2 in order to avoid the more 
serious consequences of kinetic heating, thereby easing 
not only the structural problem but also that of tem 
perature control and air conditioning in the passenger 
cabin. Since overall efficiency is rising as Mach Num. 
ber increases, figures quoted above for Mach Number 
2.5 are optimistic for the more prudent choice of Mach 
Number 2. 

The place where supersonic civil aviation may be 
come more immediately interesting is clearly at shorter 
ranges, so that the aircraft is moved out of the “ver 
difficult’”’ into the ‘‘possible’’ class. If one thinks oj 
still-air ranges of 3,000 miles and less (stage lengths a 
little under 2,000 miles 
large proportion of travellers), direct operating costs of 
about 30 pence per long-ton mile at Mach Number 2 


and these are of interest to a 


can be found on prudent assumptions of lift drag and 
specific fuel consumption, and if optimistic figures are 
taken, a price of 12 pence comes in sight. This, inci- 
dentally, shows how sensitive aircraft of this type are 
However, I do not believe that, in practice, complete 
refinement can be achieved without development, and 
I would therefore rather base assessment of the im- 
mediate possibility on figures in the region of 50 pence 
with the knowledge that, when experience is gained in 
this regime of flight, there will be room for develop- 
ment that can be favorable to economy. 

To put the above figures in perspective, Fig. 13 
contains a plot of total operating costs—.e., the direct 
costs quoted above plus an element for indirect costs 
The scale is arbitrary, and no more should be read into 
the plot than that it gives a very rough indication of 
the overall position. Cost is not the only criterion ol 
service in aviation, but it does control the section of the 
market in which the greater part of the total demand is 
likely to be found. What these figures suggest is that, 
on present established prudent 
assumptions about what we may learn in the future, 


knowledge and on 


there is no great incentive in cheapness to impel a 
movement of long-range civil aviation into the super- 
sonic regime as there has been to bring about the present 
But the 
figures also show that costs may well be achieved which 
are sufficiently good to attract the smaller section of 
the market which is concerned more with time and 


increase of cruising speeds towards 500 knots. 


comfort. 

Against this background, I return to the question | 
left earlier, which was whether the high-subsonic air- 
plane merited long-term attention to reduction in cost. 
The answer is “‘yes’’ since the other regimes now opening 
up do not offer the high-subsonic aircraft such competi- 
tion that a large part of the market can be expected to 
turn to them, anyway, for a long time. The jet-pro- 
pelled high-subsonic airplane is an ideal tool of medium 
to long-range civil aviation—or will become so when the 
ancillary services, such as traffic control, match it as 
they undoubtedly can. I would like to discuss for a 
moment the two ways mentioned previously by which 
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creater economy could be sought in this class of airc raft. 
The influence of size is not apparent from the range 
equation because it comes from a greater economy of 
weight in the design, but there is little doubt that an 
airplane flying at Mach Number 0.87 and weighing 
150,000 Ibs., as compared with the next generation's 
size of nearly 300,000 Ibs., would show an effective 
reduction in the cost of the seat-mile. But there is a 
clash between ultimate economy and the service the 
passenger expects. It is roughly the case that settled 
transport systems provide a frequency of service com- 
parable with the traveling time——a not unnatural situa- 
tion since it would be odd to have to wait an hour to 
make a 5-min. journey or, equally, to have daily de- 
partures for a sea journey of 4 weeks. In Fig. l4isa 
diagram showing the way in which transport systems 
of different types fall into the pattern of ‘‘time between 
departures roughly The 
subsonic airplane, by greatly decreasing the journey 
time, will generate a demand for higher frequency, and 


equal to journey time.” 


this, on a given size of market, must delay the need for 
greater size in aircraft with the economy that it brings. 
Nevertheless, because the market for civil aviation will 
increase so much in coming years, still bigger jet-pro- 
pelled aircraft are by no means unlikely to be developed, 
but I think this will not happen until frequency has 
first been increased over present figures. There is 
another factor 


possibility of gaining further economy by increase of 


which may influence adversely the 
size, and that is the noise generated by the aircraft. 
Noise is, of course, related to the power installed in the 
machine, and a bigger machine will also be noisier; 
this point is discussed in a later section. 

The high by-pass ratio ducted fan seems in many 
ways the ideal power plant for these high-subsonie air- 


craft if its promise can be achieved in practice. Fig. 15 
shows the relationship between all-up weight and pay 


load for a 500-knot trans-Atlantic airplane powered 
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Fic. 14. General transport problem. Frequency of de- 
parture-journey time. (Due to Sir Arthur Gouge and reproduced 
with permission from Journal of The Royal Aeronautical Society.) 
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Fic. 15. Comparison of ducted-fan and simple-jet propulsion 
for high-subsonic air liners Phe aircraft all have a still-air 
range of 5,500 nautical miles at .1/ 0.83 and are designed for 


similar runway performance Phe runway performance used is 


severe and accounts for such differences as there may be between 
these figures and published information about particular air 
craft 


This 
takes account of the higher specific weight of the fan. 
I think the reason that a ducted fan is not available is 


not because these arguments are invalid but because 


with simple-jet engines and with ducted fans. 


the main sweep of jet-engine development has hitherto 
been pressed forward by the needs of fighters, where 
cruising fuel economy is by no means the sole criterion 
for success. But in addition to this, as mentioned in 
the discussion of specific fuel consumption, the ducted 
fan is an engine that has a narrow field of dominance; 
below Mach Number 0.8 it is defeated by the propeller 
turbine, and above Mach Number | it is defeated by 
the simple-jet engine. It happens that the region ot 
the fan’s dominance is that at which long-range civil 
aviation finds its best economy and at which it is likely 
to play a large part for many years in providing travelers 
with service. The arguments of the past which con 
demned the ducted fan on grounds of lack of flexibility 
for application over a wide speed range were valid as 
they were put but, perhaps, did not take account of 
the fact that an important section of aviation will 
operate in conditions in which the fan has much to offer. 


THE LIFE OF AIRCRAFT STRUCTURES 


The only structural topic I want to discuss, apart 
from the generalizations contained in previous sections, 
is the question of the strength of aircraft, such as those 
that must have a very long working life for civil use 
or military aircraft which are especially subjected to 
fluctuating loads when used, for example, for ground 
attack. I think two matters have become particularly 
controversial: (1) Should structural design be based 
on the philosophy of guaranteed life within which 
structural failure will not occur (the ‘‘safe-life’’ policy) 
or on the philosophy of if it fails, it does not cause im 
mediate catastrophe (the ‘‘fail-safe’’ policy)? (2) Is 
the use of a fatigue test on the complete aircraft justi 
fied ? 
Aircraft Establishment used a complete aircraft fatigue 
test for a particular accident investigation, such a test 
has become mandatory in British civil airworthiness 


Since I and my then colleagues at the Royal 


requirements; this has been a subject of much discussion 
as has also the question of whether the particular form 
of test we used (the tank test) is the best way to do the 
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work, assuming it has to be done. I am grateful to 
my late colleagues, Dr. P. B. Walker and Mr. Atkinson, 
of the Royal Aircraft Establishment, for discussions 
which have helped me to clear my mind on these points. 

The safe-life approach may be defined as the de 
termination of the length of service that an aircraft 
structure, or part of it, can endure without failure in 
When the safe life of any part is reached, the 
If this 


fatigue. 
part concerned is to be replaced by a new one. 
is impossible, the aircraft has reached the end of its 
useful life and must be withdrawn from service. 

The fail-safe approach permits failures in the primary 
structure provided the failure is certain to be found by 
inspection before the strength is reduced below specific 
standards. The implicit assumptions in this approach 
are that the fail-safe characteristic can be established 
and that reduced static strength of the structure can be 
accepted. The suggestion that a reduced standard of 
strength can be accepted for the period between the 
occurrence of a failure and its detection and rectification 
raises the question of how the reduced standards are 
to be arrived at. matter which has been 
under serious official consideration recently in the 
United States and in England, and I will not intrude 
But whatever arguments 


This is a 


on the discussion in detail. 
on probability of failure during the vital period are 
raised, I doubt if one can come below a strength such 
that the pilot can safely apply at any time a maneuver 
of reasonable severity since, otherwise, the confidence 
of the crew and passengers in the ruggedness of the air- 
craft must be in doubt, anyway, throughout its later 
life. 
difficulty, and this is a point on which close agreement 


The definition of ‘‘reasonable severity’ is the 
between various authorities is very unlikely to emerge 
quickly. 

If a reduced standard of strength is acceptable later 
in the life of the aircraft, were the original standards 
too high? The standards specified in the past have, 
on the whole, resulted in satisfactory service experience, 
and it would, I think, need substantial arguments to 
reduce them. These arguments might be provided as 
greater attention is paid to antifatigue design and 
fatigue strength since an element in the specification of 
static strength is an endeavor to cover fatigue require- 
ments in the broad. At present, design tensile stresses 
in civil aircraft are perhaps more governed by considera- 
tions of fatigue than of static strength; as antifatigue 
design improves or as materials with lower static 
strength but better fatigue properties are used (as some 
designers advocate), the whole question of static 
strength requirements will need further review. 

The use to which the aircraft is to be put must de- 
termine which design approach to adopt; this matter is 
not to be seen as a fight to the death between safe-life 
and fail-safe, and it is not a good service to aviation to 
seek to make it so. For example, for smaller aircraft 
for use in outback areas with limited maintenance and 
inspection facilities, there is good reason for using the 
safe-life approach; the operator is put in the straight- 
forward position that he changes parts at a stated time, 


FCAL SCIENCES MARCH, 
and the single-spar construction may well give him thy 
more economical aircraft. For aircraft taking a heayy 
beating in continuous low-altitude use, the overall 
economy throughout the aircraft's working life may 
best be met by safe-life, the cost of the new components 
being offset by the inspection and repair costs int 
which the fail-safe design may run in its later life 
am assuming, of course, that, in designing under thy 
safe-life approach, the designer will make adequat 
arrangements for ease of replacing components. 

For larger aircraft used for the most part in the high 
altitude cruise condition and worked by operators with 
excellent inspection and maintenance facilities, bott 
approaches are perfectly reasonable, and a combina 
ation of safe-life and fail-safe may well be appropriate 
For economy in inspection and maintenance, a reason 
ably long crack-free service must be attained in a fail 
safe design, even if the crack is not catastrophic, and 
it is then essential that the cracks can be easily found 
and the part repaired. In some parts of the structure 
there may well be difficulties of inspection and main 
tenance, and these parts must be put entirely beyond 
suspicion by a safe-life approach. I am only too con 
scious that, when it comes to deciding in a particular 
case, these principles can become distinctly troublesome 
as a guide; nonetheless, both approaches have their 
virtues in particular circumstances, and they are by 
no means mutually exclusive. 

In the present state of knowledge there must, | 
think, be a test to confirm that the trouble-free life is not 
too short, to determine where cracks are to be expected 
and their significance, to determine safe lives where 
necessary, and to develop inspection techniques and 
procedures that can be relied upon to find the cracks 
Recent test work in England has 
that it is not where 


before it is too late. 
shown what one would expect 
fatigue trouble has been anticipated and designed 
against that fatigue failures are found, but in local 
details and particularly in those that cause a rapid 
diffusion of load; there have been somewhat higher 
stress concentrations than were estimated. In nearly 
all our major tests, cracks and failures have been pro- 
duced at places other than main joints, and they can 
usually be explained by considering the way in which 
the load is diffusing in the area. Concentration on the 
joints alone is not sufficient. 

The complete aircraft fatigue test provides knowledge 
which it is difficult to get with certainty in any other 
way. A longer operating life can be allowed with 
safety if test results are available to substantiate the 
design. The method includes comprehensive strain 
gaging, and it is usually possible to extend the fatigue 
life by minor modifications at points of high stress 
revealed in the tests. In the application of fail-safe 
principles, the tests produce the cracks and establish 
their rate of growth; they prove that the alternative 
stress routes actually exist and that they operate suc- 
cessfully after the crack has formed. 

Increases in the life of civil aircraft to 20 years at 
full utilization is an aim which, if achieved, would make 
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4 further impact on the already favorable trend of 
aviation economies; such achievement will be largely 
governed by success in design against fatigue and is 
therefore dependent on the acquisition of the detailed 
knowledge which only testing can provide. The more 
fatigue testing is done, the more it seems to be justified. 

We have not found anything better than the water 
tank test. 
tests done under water and those done in air. 


No difference has been detected between 
There 
ire, of course, uncertainties in the correspondence be- 
tween a tank test and service use, but I doubt if they 
are serious, and further experience will shed more light 
is difficult the 


on them It to represent 


maintenance and of such factors as corrosion that are 


time dependent. Some commentators have been con- 
cerned as to whether corrosion taking place in the tank 
may not produce a serious adverse effect on the fatigue 
life of the specimen; this is not so, for less corrosion 
takes place in the tank than occurs to the aircraft in 
service as a result of flying or standing in the rain. 
The question is whether the tank overestimates the life, 
rather than underestimates it, on this score. 

The broad indication from work in England is that 
from the fatigue standpoint the cabin is an easier de 
sign problem than the wing. On the other hand, the 
cabin 1s vulnerable to danger from external 
sources, and therefore the form of construction adopted 


more 


and the stress levels used should have this very im 
portant consideration in view. As for tail planes, our 
experience has been good; the structure is usually rela 
tively simple and straightforward, and it may also be 
that designers tend to be more conservative with the 
tail. But which require 
careful attention. 
to have a low value of the reaction coefficient, and this 


undercarriages are items 


The tendency is to design the oleo 
means that, if the undercarriage is designed under 
considerations of static strength, then a relatively large 
fraction of the ultimate load is applied on every landing. 
The use of large light alloy forgings gives this particular 
point; unless the design and heat-treatment of the 
forging is suitable, large locked-in stresses can be caused 
and lead, in conjunction with the landing load, to early 
fatigue failure. Our experience has led us to the 
following principles in handling large forgings: these 
have helped to achieve success in avoiding serious 
locked-in stresses: The forging should be quenched in 
boiling water, and the alloy used should be such that it 
suffers only minor loss of mechanical properties as a 
result of this treatment. The forging should be heat- 

the fully 
Where the 
in a boiling water quench is not good enough, a very 


treated in machined condition whenever 


possible. dimensional control obtainable 
small amount of post-heat-treated machining is allow- 
able, provided it is symmetrically disposed. 

Before leaving the subject of fatigue it is interesting 
to note that the problem of the rate of growth of cracks, 
which is vital to the fail-safe approach to design, is an 
example of an almost forgotten research. In 1921, 
Dr. A. A. Griffiths put forward an analysis of the 
propagation of cracks in brittle materials and advanced 
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the hypothesis, which he was able to support by ex 
perimental evidence, that rapid propagation of the 
crack takes place if the rate of work done by the release 
of strain energy in the material exceeds the rate ol 
This 


research was apparently forgotten, except perhaps by 


energy demanded to produce crack extension. 


specialists, and was not connected by engineers with the 
problem of crack propagation in ductile materials used 
for structures; indeed, it is probably fair to say that, 
until an aircraft accident brought the question of crack 
propagation into strong prominence, most engineers 
had not studied the question seriously at all, though 
some steel structural designers had encountered and 
considered the problem of brittle fracture particularly 
in such structures as oil storage tanks. 

Recent work has shown, as would be expected, that 
the original Griffiths theory does not apply to ductile 
materials probably because the energy absorbed in the 
plastic phase of deformation of a ductile material must 
be taken into account. However, although attempts to 
evaluate the rate of work required to create new cracked 
surface, using present theories, give disappointing 
results, a great deal of empirical data is now accumulat 
ing. It appears from present work that the relationship 
between stress and the length of crack which 1s critical 
in the sense that it will run rapidly is affected by the 
width of the plate and to some extent by its thickness, 
but not by its length. The relation is unaffected by 
evelic loading or prestraining but is greatly influenced 
by the presence of reinforcing 

In Fig. 16 are shown some typical results on the in 
fluence of crack-stopping reinforcement, for which | 
indebted to The Bristol Aeroplane Company, 
An artificial crack was made in the plate to 


am 
Limited. 
which stress was then applied. The graph shows how 
the crack grew under additional stress. Rapid growth 
starts in each case at about the threshold for an unre 
inforced plate; this suggests that the condition for the 
start of self-propagation is local in character and is un 
influenced by the presence of the stoppers. There 
after, rapid growth took place in the specimen using 
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a riveting stopper until it was arrested at the inner 
rivet line; the specimen then took a higher stress until 
it failed. 
adhesive (Redux), the growth was slower; when the 


When the stoppers were attached by an 


stopper was machined integrally with the plate, the 
growth was also slower, but failure took place earlier 
than in the reduxed specimen, occurring when the crack 
reached the fillet on the stopper. 

With the accumulation of empirical data of this kind, 
it is likely that before long the mechanism of cracking 
and tear resistance will be much more understood, and 
fail-safe design will become firmly based. 


NOIsE FROM AVIATION 


The noise on the ground from aircraft in the air is 
not unnaturally becoming a considerable issue in public 
affairs. It is sharpened as an issue by the fact that, at 
present anyway, the people on whom the noise descends 
are, for the most part, not the people who enjoy the 
benefits of aviation. I think that there is no doubt 
that those who do enjoy its benefits must be willing to 
pay a little more in order to reduce the impact on those 
who do not. In this section I want to discuss, in the 
broad only, what has caused noise to emerge more im 
portantly and what the cost to the traveler will be to 
keep it in bounds. I am indebted to Mr. Greatrex of 
the Rolls-Royce company for discussions on this matter 
that have cleared my own mind, and I have quoted 
figures from a recent paper on the subject written by 
Mr. Greatrex. 

First, it is important to realize that the size and 
cruising speed of an aircraft controls the noise it makes 
because this is related to the power of its engines. A 
large part of the present-day worry about noise is 
caused by the fact that aircraft are larger than they 
used to be. There will, of course, be detailed differ 
ences in design between aircraft which will affect com- 
parison of the noise they make, but in the broad it is 
inevitable that an increase of weight from 30,000 Ibs. 
to 300,000 Ibs. will involve an increase of noise of 
about 10 db. This is a factor which has some bearing 
on whether aircraft will get bigger and, therefore, more 
economical, as discussed in a previous section, though 
it is an Open question whether a louder noise less often 
would be preferable to a lesser noise at higher repetition 
rate. It willalso have a bearing on whether jet-elevated 
vertical take-off aircraft could be used for city center 
operation; they have an installed thrust which exceeds 
their weight and must be very noisy. 

Then there is the influence of the introduction of jet 
propulsion; the character of the noise from a jet engine 
is different from that due to a propeller engine. In a 
propeller-driven machine, the propeller itself is a major 
source of noise; this is added to by the engine exhaust 
to an extent that depends on the type of engine and on 
its installation. 

Making allowance for increase of size and cruising 
speed, on a pessimistic basis future jet-engined air- 
craft of the 300,000-Ib. class, without provision for 
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noise suppression, will be about 18 db. noisier thay 
quiet present-day pistoned types and 8 db. noisier 
than noisy present-day types. We are now confident 
in England that we can more than bridge the gap be. 
tween the large jet transport and present-day noisy 
types. The use of a by-pass engine, which gives q 
lower jet velocity than the simple jet for the same 
3.5 db. In passing, the 
Next, 


means have been found by the use of a development oj 


thrust, drops noise by about 
use of a ducted fan would drop it much further. 


a corrugated nozzle to produce a further reduction oj 
about 7.5 db., giving a combined reduction of 11 db 
This will make the future large jet quieter than present 
noisy piston types, though still rather louder than 
present quiet piston types. To produce this effect the 
nozzle must be very corrugated indeed; relatively minor 
corrugations have a much smaller effect. 

Jet engines need thrust reversers, and it is therefore 
desirable that the jet pipe should incorporate a noise 
suppressing nozzle and a thrust reverser. Successful 
designs have been evolved incorporating both, but 
these are not, of course, without their cost. It is 
difficult to separate penalties resulting from each of 
these two items, but dealing with the noise suppression 
alone and taking together the effects of loss of cruising 
thrust, the drag of the installation, and the weight 
penalty, it emerges that the overall result on a large 
long-range jet aircraft is equivalent to an increase in 
cruising specific fuel consumption of about 2 per cent 
If take-off length is critical, the effect of loss of take- 
off thrust might be equivalent to another | per cent 
increase in specific fuel consumption. The 
in the cost of the seat-mile will depend on the flight 


increase 


plan, but in very round figures the passenger will pay 
about 6 per cent more on internal routes and 9 per cent 
more on trans-Atlantic routes than the fare the aircraft 
could otherwise offer in order to relieve the people on 
the ground. I think this is a fair bargain. 

The field remains a fruitful one for research. We 
have by no means a full understanding yet of the way 
their effect, 

As under- 


in which corrugated nozzles produce 

though there are some theories in outline. 
standing is gained, it may be possible to make further 
improvements. It that the 


which it will be socially acceptable to take the technical 


may well be extent to 
steps which can give greater economy or greater flext- 
bility of operation to civil aviation will be governed by 
the extent to which we succeed in gaining these im 
provements. 
CONCLUSION 

I have not set out in this lecture to reach, or to 

advocate, any particular conclusion. But, in writing 


about aeronautical research and development and 
some of the technical factors which will influence the 
next phase of aviation, it is difficult to resist the tempta- 
tion to think in broader terms of what aviation has 
done and will do. The world has needed aviation, and 
has set out to get it, to an extent that has caused it to 


be above most other things the mainspring under which 
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the methods of organized research have developed; 
these methods are now rapidly spreading throughout 
other technologies. But the reason aviation has been 
sy much sought after--and will continue to be sos 
that it has exerted a greater influence in a shorter time 
ver a wider area of politics and sociology than most 
other developments, and its powers in this respect are 


still nearer their beginning than their end. 


Mobility of people and of heavy goods has had as great 
an effect on the broad trend of world affairs as most 
other factors. Political integration has been greatly 
influenced by it because it is the capability of moving 
men and material with reasonable dispatch to the limits 
of a political area that gives point to the rule of what- 
ever law applies within it. Internal transport and 1n- 
ternal communication systems have set a limit to what 
it was reasonable to attempt to govern, and it has 
generally been the case that, when that limit was ex 
ceeded, the outlying districts did for long—~in 
historical terms—stay within the orbit. 
, means of movement at much greater speed, with by 
capacity for quantity and with 
than 


not 
By providing 


10 means negligible a 
sreater indifference to geographical features 
most other systems, aviation has already had, and 
will in future even more exert, a decisive influence on 
world affairs. 

In recent years, civil aviation has matured and 1s in 
sight of having tools for its job which are excellently 
suited to it; that this is so will become quite clear when 
the ancillary services on the ground have been adapted, 
as they can and must be, to meet the needs of the new 
aircraft. Civil aviation has before it an enormous ex- 
pansion; this will come from a pursuit of economy in the 
cost of the seat-mile, which will succeed and which will 
open up air travel to a far wider cross section of the 
world’s population than use it now. Civil aviation has 
expanded in recent years at a rate which exceeds that 
achieved by most other industrial enterprises at a com 
parable stage in their development; on air lines operat- 
the 2 million passengers were 
carried in million in 1953. Air 
a considerable expansion; 


western world, 


1937 


ing 1n 
and over 50 
freighting has also shown 
there tenfold 
carried by western world scheduled operators between 


1946 and 1955, and this has been done largely with 


was a increase in freight ton-miles 


adapted aircraft rather than with machines particularly 
designed for the purpose. Freight will follow the pas- 
senger trend——it will go faster and will be turbine pro- 
pelled and it is likely that the market will respond to 
the increasing availability of machines particularly 
suited to its purpose. 


In the military field, aviation’s first major influence 
was to restore to land battles the mobility that they 
had lost and so to alter radically the character of land 
war. Land battles had become static in World War 
I because it took a long time to concentrate for an 
attempt at a break-through, and, during this time, the 
other side, becoming aware of the build-up of strength 
before it, was able to group to meet the attack when it 
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came. Tactical aviation changed this; improvement 
on ground transport permitted more rapid concentra- 
tion, while, if air superiority has been attained, inter 
ference from the air frustrated countergrouping. This 
was the major effect of aviation on World War II. 
But, throughout the development of long-range avia 
tion, the concept of dominance by the bomber being 
decisive in itself had been canvassed. This possibility 
was not, in fact, in sight until the development of 
atomic explosive which, by providing a vast increase of 
hitting power at a moderate weight, permitted the 
design of aircraft for its delivery having far longer 
ranges than those previously used for bombing. This 
has made possible the policy of the nuclear deterrent; 
it has also made a radical change in the distribution of 
the air threat. The political reaction of countries to 
world situations is inevitably 
threat to their homeland, and one of the big influences 
which aviation will have on world affairs in the next 
decade will arise from the fact that the United States 
has joined the other major powers in being under a 


governed by the direct 


heavy direct aerial threat. 

The extent to which the political effects arising from 
this will be moderated by technical advance in air 
defense remains to be seen. The nuclear bomb, by 
multiplying something like a million times the explosive 
power that can be delivered by a single airplane, has re 
lieved the offensive of restraint in range, in end-point 
accuracy, and in sensitivity to losses, and defense has 
therefore a large burden to overcome if it is to be of 
account; this burden is made the heavier by the fact 
that defense depends vitally on electronic techniques 
which are, to a greater or lesser degree, vulnerable to 
electronic interference. Improvement of a few per 
cent, or even hundreds of per cent, makes little impact 
on a millionfold advantage. Much has been made in 
public debate of the high-explosive 
guided weapon as an antidote to this position, but it 1s 
difficult to feel happy with the argument. There are 
available to such weapons, in principle, no devices of a 


ground-to-air 


special kind for ensuring end-point accuracy or lethality 
that are not available to fighter aircraft, and the choice 
of whether or not to replace the pilot by automatic 
devices should be seen as one of many technical de 
cisions involved in creating a weapon rather than, as 
seems to have been the case with some commentators, 
an issue in itself. As an issue it is a relatively minor 
one from the technical point of view. The overriding 
factor is the millionfold increase in offensive power. 
This will, no doubt, eventually be offset by the use of 
nuclear explosive for air defense, and such a develop 
ment may well make a more powerful impact on the 
dominance of the bomber and ballistic missile, and, 
hence, on the political results that stem from that 
While these 


consequences of developments in long-range aviation 


dominance, than any other single factor. 


are going on, the initial impact of military aviation on 
the land battle will not have diminished in importance; 
indeed, the more successful the policy of the global 
deterrent becomes, the less likely it is to be employed 
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and the more likely it becomes that means will be 


sought to ‘get around the flank.” 

In a world which must increasingly recognize the 
interdependence of countries’ vital interests and realize 
that advance in the standard of living for all is best 
achieved by ensuring that the contribution that each 
can offer, whether in terms of special goods or of special 
facilities, is available to others, transport aviation for 
people and for goods will have a decisive contribution 
to make. 
guard for stability; the civil side an important means of 


The military side of aviation is an essential 
exploiting this stability for the common good. 
APPENDIX 


THE RANGE AND ECONOMY EQUATIONS 


The general performance of an aircraft that spends a 
large part of its flight in a cruising condition in the 
stratosphere is contained in the equation 


R = (L/D)(V/c) In (Wo/ Wr) 
where 
R = gross cruising still-air range 
(LD) = lift-to-drag ratio in the cruise 

I” = cruising speed 

c = specific fuel consumption of the engines 
IV) = total weight at the start of the cruise 

Wp» = weight of the aircraft less its fuel 


This equation applies to any type of propulsion, if 
the power-plant performance is converted to an equiva- 
lent specific fuel consumption, and to any configuration 
of aircraft provided it flies at constant lift drag ratio. 

Defining the operational fixed-weight equipment as 
the sum of the earning pay load, the crew, and the 
equipment of fixed weight necessary for the conduct of 
the mission, the weight less fuel IV, of the aircraft is 
the sum of the weights of the operational fixed-weight 
equipment and the weight of those items that are 
dependent on aircraft size, such as structure, power 
plant, and systems. The range formula thus relates 
these weights by the expression 


(We Wo) + WV Wo) = e 


Discussion of 


Hugh L. Dryden, Director, NACA: 
indebted to Sir Arnold for his broad review of the re- 


We are greatly 


search and development process by which new and use- 
ful aircraft are produced and for his assessment of the 
current technical situation. The view from the other 
side of the ocean is remarkably similar to that from this 
side, although there are some differences, especially 
when each of us is examining the other's handiwork. 

It is natural that I should be especially interested in 
Sir Arnold’s comments on the nature of research and its 
administration. I note that Sir Arnold lays special 


emphasis on what he calls ‘“‘applicable research’”’ 
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where 
IW = weight of the operational fixe 
weight equipment 
Ws, = basic weight = Ws + W, IV 


Is, Wp, Wy = structure weight, power-plant 
weight, 


tems weight, respectively 


IV Wo can be loosely described as the fractional pa 
load. 
weight of the crew and fittings for the mission. 

The data contained in this equation can be plott. 
in the form shown in Fig. 6, where Rc VV’, L/D, WV, W. 
and VW, The 
speed proposed I’, and the specific fuel consumption 


WV) are related. range needed 


available fix Rc’ WV, which is then related to the lin 


marked with the 1 D expected, on the left-hand plot 
horizontal transfer to the right-hand plot then pr 
duces the relation between II’, I!) and We WW 
A rough approximation to the direct operating cost 
a long-range civil aircraft is given by 
P = (A/R)(Wr/W.) + (B/V)(Wa/W. 


where 


P = direct cost in money per ton-mile 
A, B= 


IV, = earning pay load, which is II¢ less an alloy 


constants 


ance for crew and nonearning fittings 
IV, = basic weight 


The first term covers fuel costs, while the second covers 


such items as depreciation, maintenance, and crey 
salaries, which are roughly related to the basic weight 
of the aircraft. In most cases of interest, to the a 
curacy needed for exploratory work, II”, is not greath 
different from IV¢, and 


P = (Wo/We) [((A/R) (We Wo) + (B/V) (We/W 


The weight ratios are obtainable from Fig. 6. Suc 
approximate costing is not, of course, suitable for de 
tailed work, but it is good enough to gain an initial im 
pression. (This cost formula is due to Mr. Thorne 


of the Royal Aircraft Establishment.) 


the Lecture 


i.e., that directed toward a ‘‘broad search for detailed 
knowledge, from which the possibility of a radical en 
gineering change emerges.’ One of the current issues 
here in the support of research in science and technol 
ogy and in its administration is the degree to which re 
search should be guided by practical considerations ol 
NACA we 


advantages in the use of research airplanes or concepts 


applicability. In the have found great 


of new power plants to stimulate applied research t 
advance aeronautical technology. Such procedures en 
able us to pioneer later generations of practically useful 


devices. Nevertheless, we have attempted to do re- 
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SOME COMMENTS 


search of the undirected type over a broad front also 


Qur present aim is to apply about 15 per cent of our re 
sources to this longer range outlook. 

Sir Arnold's views that the Russian aircraft, until 
recently, constituted one aerodynamic type with vari- 
‘tions and that this concentration on sweptback wings 
with fences led to economy in research are very inter- 
esting. However, our own experience has been that 
seemingly minor differences in the air-frame shape often 
have surprising effects on stability and control or even 
m drag. In fact we do not see how completely suc 
cessful supersonic airplanes can be built in any reason 
ible time without access to facilities capable of repro 
ducing the full-scale Mach Number at high Reynolds 
We are gratified to see that our British col- 
I agree fully 


Number 
leagues have now obtained such tools. 
with the absolute necessity of individual brilliance and 
managerial skill in aircraft research and development, 
but some jobs simply cannot be done with primitive 
tools. 

Trained manpower of suitable quality seem to be in 
short supply on both sides of the Atlantic. It is en- 
couraging to learn that in the past agitation of this 
question produced results; for we are doing much talk 
ing and taking some actions in this country which we 
hope will be fruitful. 

Sir Arnold emphasizes the technical risks that must 
be taken to maintain sufficiently rapid progress to stay 
ahead of our opponents and rightly calls attention to 
the undesirable effect of compounding such risks simul- 
taneously in all components of a weapons system. 

I shall in the main leave comments on the assessment 
of the current technical position to others. We in the 
NACA along with others in this country have made 
studies of the Breguet range equations and of the effi- 
ciencies of jet engines similar to those described in the 
lecture. The results depend on the assumed values ot 
various parameters representing the present and po- 
tential state of the art. As illustrated in the lecture, 
such studies are a useful guide to needed research. 
They also demonstrate the way in which the design is 
sensitive to the mission requirements stated by the 
user. 

I was interested to hear the conservative position of 
the lecturer on vertical take-off and landing of air- 
craft and on aircraft noise. 

Sir Arnold has added another to the long list of 
iruitful exchanges of technical information and opin- 
ions between members of the Royal Aeronautical So- 
ciety and the Institute of the Aeronautical Sciences, 
which have been so beneficial to the progress of aero- 
nautics in our two countries and to the security of free- 


dom-loving people everywhere. 


Clifford C. Furnas, U.S. Assistant Secretary of De- 
lense for Research and Development: The author is 
to be complimented most highly on his meaty and 
cogent dissertation. He has favored us with a truly 
significant analysis of many of the most important 


phases of modern aviation. He has been signally suc- 
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rorPrics ISS 


cessful in avoiding excessive attention to minutiae, 
on the one hand, and overgeneralization on the other 
He has given us what we so often need but so seldom 
receive—a balanced picture, pleasingly highlighted in 
the right places and distinguished by both depth and 
breadth. 

with his mention, at 


I was particularly pleased 


several points, of the inherent virtues —only partially 
neutralized by its limitations—of the ducted-fan engine. 
This presentation should be particularly noted in this 
country because we have treated the ducted fan as an 
undesirable stepchild for far too many years. But it 
has persisted in not disappearing completely. I trust 
that Sir Arnold's observations will be instrumental in 
bringing this development up to its proper position of 
respect and dignity. 

Although he might not welcome the suggestion, I do 
wish that we might have the benefit of a series of lec 
tures by the same author rather than the single sample 
And I 


would like to suggest some topics that I am sure he 


we have had the privilege of hearing today. 


would handle in most admirable fashion. 

An intensive exploration of the optimum equipment 
and methods for the future of the air-freight activity 
would be of interest, and probably fruitful. Though 
speed is a variable of value in air freight, it is not as 
precious as in the case of passenger transport; nor 1s 
comfort. As a segregated operation, air freight would 
have lower fixed and overhead charges. Frequency of 
schedule is probably of lesser importance than in the 
passenger business. Would the optimum solution in 
volving these and the other parameters indicate the 
desirability of a different size and type of aircraft for 
purely freight service than for the passenger case? 
Or are some of anticipated passenger or military trans 
port types sufficiently close to the optimum to enable air 
freight to achieve its deserved importance in the near 
future? I am sure the author's thoughts would give 
valuable guidance. 

In another area, it would be most interesting to in- 
vestigate the implications of air-borne passenger trans- 
port at truly high altitudes—say, above 70,000 ft. 
Such an exploration would, of course, point up the new 
developments that would have to be carried to suc- 
cessful conclusion—in power plants, air frames, ete. 
before the idea would have the ring of practicality 
But given success in component development, would it 
be economically feasible to exploit the parameter of 
altitude to the full for passenger transport for medium 
to long distances? I am not sanguine of the economic 
feasibility of passenger transport at supersonic speeds 
unless we are able to take substantially more advantage 
of the third dimension than we have done to date. A 
careful analytical, projective study of the possibilities 
would seem to me to be most timely. 

In his discourse, the author recognized the potential- 
ities of vertical take-off, fixed-wing aircraft for certain 
special purposes but expressed serious doubts about 
their adaptation to regular civilian use. Granted the 
drawbacks of decreased overall performance and in- 
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creased noise, I should like to see the problem reana- 
lyzed for the short- and medium-haul cases with full 
recognition of the debilitating effect of ground conges- 
tion on human transport and the importance of office-to- 
office time. <A trans-Atlantic passenger does not ob 
ject particularly to an hour's elapsed time from airport 
to office, but the Boston to New York traveler takes 
quite a different view of the matter. He might even 
be willing to pay a little more to go a little slower in 
the air if his take-off and landing locations were more to 
his convenience. A substantial portion of potential 
short to medium haul passenger traffic will never be 
realized until we eliminate the great time loss of ground 
transportation. Intuitively I feel that fixed-wing air- 
craft with vertical or very short take-off characteristics 
will offer a better solution than rotary wing aircraft. 
Very low approach speeds may also be an important 
factor in maintaining adequate safety in the increas- 
ingly congested air space. If such a development is to 
proceed, the noise problem must, of course, be alle- 
This will call forth a most difficult research 
task. We need much more fundamental knowledge of 
the gas dynamics of noise generation. What are the 
By largely empirical methods we 
with 


viated. 


kinetics involved ? 


have made some progress in noise alleviation, 


modest penalties on performance. How much better 
could we do if we really understood the phenomenon? 
We need to find out. 

Another area in which our author's analytical acu- 
men would be most beneficial would be in the trouble- 
some problems of air traffic control. In foreseeing the 
injections of substantial numbers of jet aircraft into the 
civilian transport structure he expressed great confi- 
dence that the traffic control problems would be solved. 
I hope that this faith is fully justified, but the evidence 
in hand at present does not lead one to be very hopeful 
about an adequate solution in the near future. Some 
of our airports are now averaging about 1,000 take-offs 
and landings per day, a large proportion being under 
instrument conditions, and with present ground equip- 
ment and methods this appears to be about the satura- 
tion point. Most of us are now familiar with the spate 
of cancelled flights and delays of hours which are becom- 
ing quite common on the eastern seaboard when a 
widespread layer of only moderately bad weather comes 
in. These delays and cancellations are evidence that 
the traffic handling system is supersaturated. They 
are in the interest of safety-—and that is an exchange I 
am very happy to go along with—but as air traffic in- 
creases, we must be able to handle more flights and 
retain the safety. Despite the fact that the safety 
record is remarkable, we are always conscious of the well- 
advertised ‘four near misses per day’’ in the vicinity 
We wonder how many more 
We do not readily forget the 


of American airports. 
there may have been. 
recent collision of two air liners over the Grand Canyon. 
But these deficiencies are not due to the pilots or the 
traffic controllers on the Considering the 
communication and control equipment and systems 


ground. 


available, the personnel involved, in my opinion, do a 
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phenomenally fine job. The system actually works 
better than it has any right to 


and dedication of the pilots and ground controllers 


only because of the skjJ 


The most hopetul part of the picture is that a great 
number of individuals, corporations, and agencies ar 
constructively worried about the problem.  Researe! 
and development activities for control system com 
ponents are manifold. Many systems analyses, opera 
tional investigations, and generalized studies are under 


way. But the problems are technically so compli 
cated and are so intertwined with so many nontechnical 
aspects that no clean-cut solution is in sight. What 
new researches and new analyses are needed? What 


will be the impact of inadequate traffic control on the 
expected rate of growth of air travel? I would like t 
see our author's well-demonstrated objectivity focussed 
on this broad problem. 

Finally, I should like particularly to compliment the 
author on his observations on the characteristics, and 
sometimes foibles, of scientists and engineers in his 
section on “‘The Research and Development Process, 
One can tell that he has been there. His keen analysis 
of the process, the desirable and undesirable attitudes 
the pitfalls, the necessary relations between different 
kinds of research organizations, the hazards of creeping 
decay or dry rot, and the elements of the occasional 
triumphant success is a most admirable evaluation 
Would that many important executives would read and 
I trust that they will. This will contribute 


much toward assuring the continuing constructive and 


absorb it. 
rapid progress of aviation in both of our countries. 
William Littlewood, Vice-President, Equipment Re 


Sir Arnold Hall’s lec- 
restrainedh 


search, American Airlines, Inc.: 


is analytical, philosophical, and 


It is difficult to find much to disagree with 


ture 
prophetic. 
in his comprehensive coverage of a widely diversified 
a few 


list of aviation subjects. There are, however, 


implications—-or more properly emphases——which seem 
mildly controversial or which the discusser’s prejudices 
cause him to view somewhat differently. 

Sir Arnold makes quite a point of the lack of new 
discovery or invention in aviation history. It is be 
lieved that this critical evaluation does not give due 
credit to the thousands of successful researches and 
inventions peculiar and subordinate to each of the major 
areas of improvement over the years, each of which has 
made a small but essential contribution to the phenom- 
enal overall advance in aviation performances. Our 
first 50 years of flight have, it is true, been significant 
for development and refinement rather than for basic 
invention, but, as in all other areas of scientific advance, 
there are few inventors but many researchers, develop- 
ers, and engineers. 

The great technical needs of air defense are apparent 
to all, as Sir Arnold points out, but it seems significant 
to note that, in technical air defense equipment and 
operational development, the solution of the essential 
performance, control, and operating problems will be 
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ef direct and incalculable assistance in the solution of 
the essentially similar problems of civil aviation. 

It seems desirable to emphasize even more strongly 
than does Sir Arnold the basic need for solutions to the 
ancillary problems of aviation. It has become increas- 
ingly clear that the bottlenecks to aviation growth and 
progress during the next decade lie in the functional 
and capacity restrictions of operation both in the air 
and on the ground. This is particularly true with re- 
spect to civil aviation, but since both civil and military 
operations from here out will tend more and more to 
occupy the same air space, the problems are mutual. 
Airport and airplane operational capacity restrictions 
need no further emphasis here. 

Sir Arnold predicts that civil aviation costs will come 
down. We could wish to accept this conclusion but 
must point out that the continuing inflationary rises 
in all ‘‘costs of doing business’ keep on opposing the 
hoped-for economies. It also seems significant to point 
out that, over the years and for the indefinite future, our 
aircraft equipment and services continue to improve 
in performance, comfort, and safety, always attempting 
to give more at the same price. How long this can be 
maintained in the face of rising cost elements seems 
uncertain. 

It is believed that Sir Arnold has not properly empha- 
sized the problems and importance of the growth of air 
freight. To look for the automatic development of this 
potential business giant by improvements in equipment 
alone seems asking too much. Somewhat better freight 
airplanes can be developed, but their contribution to 
reduced costs will be small in comparison with the 
business building attributes of sales promotion based 
on adequate schedules, guaranteed deliveries, and prac- 
tical business solutions of the pickup and delivery prob- 
lems. 

It would appear that Sir Arnold has omitted some of 
the most important considerations essential to the 
achievement of desirable radical technical departures 
and improvements. It has often been said that neces- 
sity is the mother of invention, and urgency of one form 
or another is probably the greatest stimulus, if not to 
invention itself, at least to the applied research, de- 
velopment, and engineering which results in radical 
technical departures. The leniency which inventors 
and researchers frequently seek under the kindly rule 
that results are hoped for but not required does not pro- 
vide the incentive urgency to achieve reasonable results 
With the 


possible exception of pure invention and basic research, 


at reasonable dates or at reasonable costs. 


which are necessarily somewhat vague in definition, 
applied research and development as well as the more 
obvious engineering should be established on bases of 
clearly defined tasks and dated requirements. We are 
endeavoring to compete with a political system which 
does just that, much more forcefully than we would 
consider democratic or reasonable. Our system does 
attempt to substitute substantial remuneration and 
other forms of less tangible honor and reward to supply 


the drive necessary for reasonable accomplishment, but 
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man being what he is, there is no equivalent to a clearly 
defined and firmly stated requirement as to what is 
expected, even if the requirements are of necessity 
modified from time to time as the desired results are 
not achieved. Programs are frequently slowed down 
by stubborn resistance to change, pride of authorship, 
and reluctance to admit error. Engineers and scientists 
are very human, and those associated with aviation are 
no different than the rest. 

Sir Arnold refers to the valuable contributions of test 
models in expediting development and engineering 
This evaluation is acceptable but is no substitute for 
the economies and achievements of imaginative and 
thorough prior thought. Building ideas into models 
too soon proves mistakes early, but it also frequently 
accepts and perpetuates an answer which is less than 
the best. Mental inertia and economics tend to per- 
petuate the ‘“‘good enough.” 

I cannot too strongly agree with Sir Arnold on the 
essentiality of close contact and understanding between 
the research, development and design activities, and the 
users of the equipment. Research needs the practical 
touch of close association with the problem environ- 
It is believed that Sir Arnold overemphasizes 
Invention and 


ments. 
the importance of strength of resources. 
originality are to a large degree independent of physical 
resources. The very significant and substantial con 
tributions of British aeronautical research and develop- 
ment, in the face of restricted resources as compared 
with those currently available in America, provide an 
outstanding example of the case in point. These out- 
standing British aviation 
particularly during and since World War II, must be 


contributions to progress, 
attributed to the existence of scientific and engineering 
originality and aggressiveness prompted by great 
urgency and to some extent by the limitations of phys- 
ical resources and—it is believed—even more to the 
existence and implementation of a broad national avia 
tion program. This program has recognized a sensible 
and proper balance, including the needs of both military 
and civil aviation. Such a coordinated program has 
been entirely lacking in the American picture, and we 
have inevitably paid the price of our rugged individual 
ism which, lacking a national emergency, would not be 
subordinated to the good of truly national aviation. 
Sir Arnold refers to the shortages of technical man- 
power—a situation which does not appear too different 
in Britain and America, both leaders in individual self- 
determination. With no immediate solution in sight 
to the problem of creating substantial additional num- 
bers of scientific and engineering manpower, it would 
seem that our great obligation is to make far more 
effective use of the manpower and technical abilities 


The 


ineffectual use and frequent waste of our scientific and 


that we have and can have in the near future. 


engineering capacities by poor planning and lack of 
coordination should cause us the greatest concern. It 
might even appear that overemphasis on specialization 
in technical education tends to deprive our graduates 
of the broad-base scientific and engineering trainings 
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so essential to evaluating and accomplishing the un- 
known. Another very important factor in sensible and 
effective accomplishment is the elimination of too hasty 
and ill-considered programs. These are inevitably 
expensive in all costs, including manpower, and all too 
More 


early and thorough study with reasonable perspective 


often result in a complete waste of capacity. 


and adequate breadth of vision will frequently result 
in substantial economies and provide earlier and more 
effective accomplishment. 

It is easy to agree with Sir Arnold’s emphasis on the 
dangers of overplaying the system concept. In an ex- 
tremely complex and developing field like aviation, it is 
certainly desirable to make progress by slow stages 
with a relative minimum of all new essential elements. 
The American tendency to overtool has also proved a 
serious deterrent to early accomplishment and has fre- 
quently demonstrated the problems and costs of in 
flexibility. 

We are indebted to Sir Arnold for some new slants 
on economic and performance formulas to help in pre- 
liminary evaluations of proposed new aircraft. We 
could certainly concur with his emphasis on the prob- 
ability of civil aircraft performances being restricted to 
high-subsonie values for a substantial period of time. 
Such performances will for some time meet most trans- 
portation needs and minimize the economic and opera- 
tional problems which threaten to be increasingly seri- 
ous in the near future. It is agreed, however, that con- 
tinued research in the areas of higher performances is 
urgently needed. We hope for further verification of 
Sir Arnold's conclusions on the promise of suitable by- 
pass or ducted-fan power plants for high-subsonic civil 
aircraft. 
been devoted to this and other possibly beneficial types 


It is agreed that not enough attention has 
of power plant for civil use. In general, we can certainly 
agree with the inherent benefits of decoupling—that is, 
the separation of the propulsive and thermal elements 
of power-plant design. 

Sir Arnold’s discussions of STOL designs, jet noise 
problems, and the philosophy and design of aircraft 
Each of these 
With re- 
spect to fatigue, a determination of reasonable load 


for long fatigue life are most interesting. 
subjects would merit a full paper on its own. 


assumptions seems one of the most difficult problems of 
the safe-life concept. It is agreed that both safe-life 
and fail-safe principles must be incorporated in modern 
aircraft designs. It is interesting to note that, whether 
we admit it or not, the fail-safe philosophy which has 
successfully dealt with innumerable minor and more 
serious fatigue failures in the past and has permitted 
us to achieve the excellent record of dependability now 
being experienced with modern aircraft is an inherent 
and undeniable element of the picture. To the extent 
that we intelligently improve and provide these fail-safe 
fatigue characteristics and to the extent that we supple- 
ment them with some reasonable examinations of safe- 
life based on sensible load assumptions and dependable 
determinations of laboratory fatigue life, we will con- 
tinue to improve our fatigue safety and economy. It 
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is not believed that British and American engincering 
are very far apart on a true understanding of the signif. 
cance of these apparently controversial fatigue matters 

Sir Arnold Hall, in his 1956 Wright Brothers Lecture 
has made an invaluable contribution to thought and 
knowledge in the ever growing field of aviation and 
has added a very worthy document to the library oj 
Wright Brothers Lectures. 


T. P. Wright, Vice-President for Research, Cornell 
University: We have been privileged to hear another 
brilliant Wright Brothers Lecture. Sir Arnold Hal 
has continued the traditional excellence of these pres. 
entations. I am more than ever convinced of the 
soundness of the tradition of alternating the lectures 
between speakers from the United Kingdom and the 
United States. 

In Part (I), Sir Arnold discusses research and de 
velopment in aeronautics, indicating that research may 
be organized or unorganized, but, I hasten to say, by 
this he does not mean disorganized. This appears t 
mean a combination of evolutionary progress with oc- 
casional revolutionary advances or discoveries. He 
gives an excellent analysis of applied research and its im- 
portance, although here I would have expected a little 
fundamental research. I somewhat 


more stress on 


question the extent to which, and I quote, “Broad 
science, on which aerodynamics is based, was to a large 
extent completed before the Wrights flew the first 
airplane.” 

The discussion of the development stages in the prog- 
ress of a design is also well worth careful study. The 
need for risk taking is appropriately noted. One is re- 
minded of the recent quotation relevant to current jet 
re-equipment reported to have emanated from an air- 
line president, ‘‘We are buying airplanes that haven't 
been fully designed, with millions of dollars we don't 
have, and are going to operate them off airports that 
are too small, in an air traffic control system that is too 
slow; and we must fill them with more passengers than 
we have ever carried before.’ Well, it appears from 
this lecture that there will have been justification for 
the chances which taken, from both the 
standpoint of the efficiency of the jet aircraft to be de- 
veloped and the potential of the market to ensure pas- 


have been 


sengers to fill the seats provided. 

Emphasis is laid in the paper on the time factor in 
development and the very long period that transpires 
between initial concept and production. It appears 
that the Russians have reduced this time more than 
has been the case in the western countries. This he 
attributes partly to the fact that they have followed a 
similar configuration in all their designs. Nevertheless, 
I think our greater willingness to explore many design 
potentials will pay off in the long run, and we are well 
advised to have allowed greater freedom for the exer- 
cise of the ingenuity of our company designers in intro- 
ducing innovations. 

In dealing with education, Sir Arnold shows by means 


of rather disturbing graphs that the increase in tech- 
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SOME COMMENTS 


scientific education for Russia has ex- 


This, of 


ind 


nological 
ceeded that of our countries. course, iS a 
problem that we in education face and have only been 
able to answer so far by stressing the quality rather 
than the quantity of our graduates in science and en 
gineering. However, remembering the all-important 
need for war prevention and, as well, the need for bal- 
in the educational qualifications of our graduates 


ance 
for citizenship, I think we should not forget the impor- 


tance of liberal arts studies, including disciplines dealing 
with human behavior, the social sciences, government, 
and others which may supply us with people who can 
contribute to the solution of the problem of achieving 
a peaceful world through better understanding of inter 
national relations. What we is more 
technical grounding of our graduates in the liberal arts 
ind more knowledge of liberal arts subjects by our 


need, I guess, 


scientists and engineers. I am reminded of a quotation 
from Dutch Kindelberger: 


the best policy to adopt the course that is best tech- 


“It may not always be 


nically, but those responsible for policy can never form 
a right judgment without knowledge of what is right 
technically 

I think that those involved in aircraft development in 
this country, both in government and in industry, 
should weigh carefully the comments Sir Arnold makes 
concerning the weapon system concept and the implica 
tions that, in practice, this can be carried too far to 
our detriment in time of development and possibly in 
excellence of results. 

The discussion on production interested me particu- 
larly. This dealt with the degree of tooling which is 
justified in an art which takes so long in the develop- 
ment of new models and in which production is rela- 
tively so small. We had a good example of this during 
World War II when the automobile industry was intro- 
duced into the aircraft production picture and, in 
many instances, proceeded to overtool. This caused 
considerable delay in getting models into production 
in their shops and a great deal of trouble as regards 
necessary changes in production lines due to required 
modifications. I think we must keep in mind this 
flexibility to which Sir Arnold refers but of which, in 
general, I feel our aircraft industry is aware. 

In Part (II) of his lecture, Sir Arnold first discussed 
the Breguet range formula and its extension to bring in 
the factor he denotes as ‘‘earning load.'’ He shows the 
great importance of L/D and reduction of wave drag 
when considering aircraft design for supersonic flight. 
The difficulties of improvement in these factors and 
in reduction of specific fuel consumption as well, are 
dwelt upon, but he gives some hope of progress in both 
areas so as to make it possible to get into his “‘very 
hard” area of Fig. 6. In this section of the lecture, 
reference is made to the possibilities of vertical take-off, 
and, indeed, it appears that, for an aircraft flying at 
Mach 
for such utilization by some means of direct propulsive 
lift. Dr. von Karmén has pictured the future as the 
age of “‘the motorist as distinguished from the aero- 


Number 2.5, adequate power will be available 


ON CURRENT 


AVIATION TOPICS LS7 


dynamicist.”’ Relevant to vertical take-off, Sir Arnold 
mentions the Rolls Royce concept using clusters of lift- 
ing jets, which seems to have a great deal of merit when 


considering aircraft in the Mach Number 2.5 


range 
But air transportation at such speeds will be costly to 
passengers. 

For high-subsonic speeds, I was pleased to note his 
emphasis on the benefits to be derived from the use of 
the by-pass engine, and we will look forward with interest 
to flight performance when comparing pure jet-equipped 
DC-8’s with those powered with Rolls Royce Conways, 
to be placed in service in Trans-Canada Airlines. 

Sir Arnold referred to the speed record of the Fairey 
F.2, and I would like to pause to pay tribute to Sir 
Richard Fairey, who was my very good friend since the 
early 1920's. He contributed greatly to the advance- 
ment of aviation and, as well, to the furtherance of close 
relations between this country and Britain during World 
War II and at many other times. 

Structural fatigue is a matter dealt with at considera- 
ble length in the lecture, and a clear discussion is pre 
sented on the qualities defined as “‘safe-life’’ and “‘fail- 
safe.’’ It is brought out that either concept, as a guide 
in design, may be appropriately followed, depending 
on the specific aircraft mission, and that, also, some 
type of overall fatigue life test such as that developed 
at Farnborough should be conducted. It is probable, 
it appears to me, that it will be desirable for the de 
signers to give careful consideration to a combination 
of the safe-life and fail-safe philosophies. 

As to the effect of noise, Sir Arnold seems to be 
optimistic that sufficient reduction, although hard to 
obtain, can be attained in our coming jets to achieve 
acceptable good-neighbor relations with residents near 
airports. He cites the by-pass engine as being less 
noisy than jets, and he emphasizes the definite need 
for thrust reversal. 

In the Conclusion, some excellent points are made of 
the growing influence of the air in military affairs. 
Whereas there was but small impact in World War I, 
by the time of World War II aviation was the deter- 
mining factor in the outcome. As to any future major 
wars, I think we may infer from Sir Arnold's lecture 
that, with the atomic weapon some millionfold more 
powerful than previously developed explosives and with 
present supersonic delivery capabilities—a melding of 


we have arrived at a point 


the jet and atomic ages 
where the age-old conflict for supremacy between of- 
fense and defense has been settled in favor of the former. 
This condition will probably persist as we have passed 
a certain threshold of offensive potential that lies be- 
yond defensive This all that 
World War III must never occur! 

Similarly, the tremendous influence of civil air trans- 


capabilities. means 


port on world affairs is now apparent and is certain to 
increase because of its speed and obstacle avoidance 
capabilities. This may well be decisive in shaping 
world affairs. When reading there 
flashed across my mind Mr. Dulles’ recurring hops across 


this passage, 


(Continued on page 202 








Investigation of Local Laminar Heat Transfer 
on a Hemisphere for Supersonic Mach 
Numbers at Low Rates of Heat Flux’ 


I. KOROBKIN* anv K. H. GRUENEWALD** 
U.S. Naval Ordnance Laboratory 


SUMMARY 


Pressure and local heat-transfer measurements were made on 
the hemispherical nose of a 2-in.-diameter hemisphere-cylinder 
with laminar boundary-layer flow. The tests were conducted in 
the NOL 40- by 40-cm. Aeroballistics Tunnel No. 1 with atmos- 
pheric supply conditions. Pressure distributions were measured 
for the Mach Number range of 0.26 to 4.87, and heat-transfer 
measurements were made for the Mach Number range of 1.90 to 
1.87. The heat-transfer measurements were made with iso- 
thermal temperature distributions over the hemisphere and 
indicate that local heat-transfer coefficients are independent of 
surface temperature in the range of +50°C. However, the 
heat-transfer coefficients must be defined in terms of ‘“‘effective”’ 
temperatures rather than adiabatic wall temperatures. These 
“effective’’ temperatures, though lower, probably do not differ 
by more than 5 per cent from adiabatic wall temperatures. 
From local velocities outside the boundary layer as determined 
from static pressure measurements, it was found that the velocity 
could be made nondimensional so that a single curve represented 
the velocity distribution along the hemisphere for the entire 
Mach Number range. Using this velocity distribution, a single 
theoretical curve has been derived and used to correlate the local 
supersonic heat-transfer data in the Mach Number range 1.90 
to 4.87. 


SYMBOLS 


a = sonic velocity (m./sec.) 

( = specific heat of model material (k-cal./kg., °C.) 

Cp = specific heat of air at constant pressure (k-cal./kg., °C.) 

D = body diameter (m.) 

d = thickness of thin-walled model (m.) 

g = gravitational constant (9.81 m./sec.?) 

h = local convective heat-transfer coefficient (k-cal./hour, 
m.%, °C.) 

k = thermal conductivity of air (k-cal./hour, m., °C.) 

l = length (m.) 

m = model segment weight (kg.) 

q = time rate of heat transfer per unit area (k-cal./hour, 
m.?) 

r = temperature recovery factor (JT, — 7,)/(T) — 7)) 

S = surface area (m.?) 

T = temperature (°C., °K.) 

T.;; = effective temperature (°C., °K.) 

t = time (sec., hour) 

U velocity (m./sec., m./hour) 

w = specific weight of model material (kg./m.*) 


Received January 6, 1956. Revised and received August 30, 
1956. 

+t Based on Laminar Heat Transfer Characteristics of a Hemi- 
sphere for the Mach Number Range 1.9 to 4.9, presented by I. 
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tute, June, 1955. 

* Chief, Thermodynamics Branch, Aerophysics Division. 

** Formerly, Chief, Thermodynamics Branch. Now, with 
Fasertechnisches Institut, Frankfurt-Hoechst, West Germany. 


x = model contour length measured from stagnation point 
(1m. ) 

a = angular position on hemisphere measured from stag- 
nation point 

B = (dU,/dx)\,-) = gradient of local velocity evaluated at 
stagnation point (1/sec.) 

mm = absolute viscosity of air (kg., sec./m.*) 

v = u/p = kinematic viscosity (m.*?/sec 

p = density of air (kg., sec.?/m.*) 

Mo = U/a = Mach Number 

Nu = hl/k = Nusselt Number 

Pr = pc,g/k = Prandtl Number 

Re = Ul/v = Reynolds Number 

St = h/pg Uc, = Nu/Re Pr = Stanton Number 


Subscripts 


a = conditions for zero heat transfer 

D = diameter used as significant length 

e = adiabatic wall conditions 

1 = local conditions in the flow at the outer edge of boundary 
layer 

() = supply or stagnation conditions 

s = modelsurface conditions 

x = contour length measured from hemisphere stagnation 


point used as significant length 
© = free-stream conditions at infinity 


INTRODUCTION 


H™ TRANSFER to blunt-nosed bodies has recen 
attracted considerable interest due to the use. 
radomes on high-speed vehicles. Furthermore, the use 
of blunt-nosed sections tends to alleviate the aerody- 
namic heating problem, since the heat flux per unit area 
to blunt bodies is far lower than that to sharply pointed 
bodies. In addition, because of geometrical consider- 
ations, the reduction in heating rate for the blunt body 
is accompanied by an increase in heat capacity, the 
measure of a body's ability to absorb thermal energy. 
The presence of a region of subsonic flow in front of a 
hemispherical nose at all speeds led to the considera- 
tion of the incompressible case as a first approach to 
the heat-transfer problem, at least for the region near 
the hemisphere stagnation point. Several theoretical 
analyses and experimental investigations have been 
performed for laminar heat transfer on spheres at sub- 
sonic speeds,'~* and in recent years the research has 
been extended to laminar flow at supersonic speeds.** 
Stalder and Nielsen‘ determined average laminar heat- 
transfer coefficients on a hemisphere-cylinder at Mach 
Numbers ranging from 0.1 to 5.0. Stine and Wanlass° 
developed an approximate theory based on _heat- 
transfer solutions of the Falkner-Skan’ incompressible 
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wedge flow model and applied it to compressible axi- 
symmetrical flow incorporating the existing pressure 
gradient of the flow around the body at the outer 
edge of the boundary layer. They also determined 
local heat-transfer coefficients experimentally for an 
isothermal temperature distribution on a hemisphere- 
evlinder at a Mach Number of 1.97. 

The purpose of this paper is to present the results of 
both experimental and theoretical research on a 2-in.- 
diameter hemisphere at Mach Numbers up to 4.87. 
The wind-tunnel tests consisted of local heat-transfer 
measurements on a segmented copper model in the 
Mach Number range of 1.90 to 4.87 and of pressure 
distribution measurements in the Mach Number range 
In addition, heat-transfer check meas- 
nickel 


of 0.26 to 4.87. 
urements were made employing thin-walled 
shells at a Mach 3.21. The supersonic 
measurements were in the Reynolds Number (Rep = 
UD v.) range of 1.7 X 10° to 7.8 & 10°. An approxi- 
mate theoretical method developed by Drake® for axi- 


Number of 


symmetrical incompressible flow and based on an equiv- 
alent wedge flow model was used to derive a single theo- 
retical curve describing supersonic heat transfer to the 
hemisphere. The Drake’ differs from the 
Falkner-Skan’ wedge model used by Stine and Wan- 
lass’ by considering the actual growth of the boundary 


model 


layer along the body. 


THEORETICAL CONSIDERATIONS 


Newton's law of cooling is normally used to define 
the convective heat-transfer coefficient for energy ex- 


change between a fluid and the surface of a body. The 
defining equation is as follows: 
q = h(T, — T;) (1) 

iere 

q = heat transferred per unit time per unit area 

h = heat-transfer coefficient 

7, = surface temperature 

7, = temperature for zero heat transfer 


For compressible flow along a body with zero pressure 
gradient, such as a flat plate or cone, 7, is identified as 
the adiabatic wall temperature, 7,, the temperature 
attained on the surface of a perfect thermal insulator. 
When considering the laminar flow of gas along a body 
with zero pressure gradient, an isothermal surface 
temperature distribution, and small rates of heat 
transfer, it is known that the heat-transfer coefficient 
is approximately independent of Mach Number up to 
\f = 4. This means that the compressible data can be 
correlated with expressions derived from incompressible 
laminar theory. 

The adiabatic wall temperature for a body with no 
pressure gradient is constant along the body contour, 
whereas the adiabatic temperature varies along a body 
with a streamwise pressure gradient. When consider- 
ing heat transfer to such a body, it is possible to de- 
termine points on the body where g is zero locally for 
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temperatures other than the adiabatic wall temper 
atures. This effect is analogous to the case of a flat 
plate with a nonisothermal temperature distribution. 
This means that there is no unique way to designate 
the temperature for zero heat transfer in Newton’s law 
of cooling. However, if the problem is restricted to the 
case of heat transfer to bodies with pressure gradient 
but isothermal surface temperatures, it is possible to 
specify a temperature for zero heat transfer, known as 
“effective’’ temperature, which will always keep the 
heat-transfer coefficient The ‘‘effective’ tem- 
perature, 7\,,, at some point on an isothermal surface 


finite. 


is determined by choosing a temperature level such that 
the local heat-transfer rate is zero at the point under 
consideration. This means no heat is transferred at 
the point in question, whereas heat transfer may occur 
at other points along the body. This is in contrast 
with the adiabatic wall temperature, which requires 
zero heat transfer everywhere on the body. 

As mentioned previously, the selection of the zero 
heat-transfer temperature for the zero pressure gradient 
case makes it possible to correlate compressible data 
with incompressible laminar theory for small heat- 
transfer rates and Mach Numbers up to 4+. The choice 
of 7.,, for the body with pressure gradient also allows 
the use of incompressible theory for correlation. A 
more detailed description of 7\,, and the reasons for 
its use can be found in references 9 and 10. It should 
be pointed out that the “effective” temperatures for 
the case of the hemisphere are not greatly different 
from adiabatic wall temperatures (the exact difference 
depends on the pressure gradient and Mach Number), 
and the distinction need be considered only when 
(1, — T,) is small, perhaps of the order of 25 per cent 
of (T, — T. In the calculation of the heat 
transfer to an actual high-speed missile, the tempera- 
substantial, 


) or less. 


ture differences under consideration are 
and the designer need only determine the adiabatic 
wall temperature which may be readily computed. 


EQUIPMENT, TEST PROCEDURE, AND DATA REDUCTION 
METHOD 


The principal test facility used in the investigation 
was the NOL 40- by 40-cm. Aeroballistics Tunnel No. 
1, which is described in reference 11. The tunnel oper- 
ates intermittently on dried atmospheric air with a 
10-60-sec. blowing time and exhausts to a vacuum 
sphere. All models used for the tests were 2-in.- 
diameter hemisphere-cylinders oriented in the tunnel 
at zero angle of attack. The distribution 
model had eleven static pressure taps distributed on 


These taps were displaced at 18° in- 


pressure 


the hemisphere. 
tervals from one another, the angles being formed by 
the intersection of a hemisphere radius through the 
tap and the longitudinal axis of the model. Pressure 
measurements were made along the hemisphere sur- 
face for the purpose of determining local flow conditions 
at the outer edge of the boundary layer for free-stream 
Mach Numbers of 0.26, 0.49, 0.69, 0.83, 1.22, 1.90, 
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2.15, 2.48, 2.80, 3.26, 4.20, and 4.87. The pressure 
data were recorded on photographs taken of a special 
mercury manometer bank designed at NOL and de- 
scribed in reference 12. Using stagnation conditions 
behind a normal shock, these pressure data were reduced 
to local flow conditions by means of isentropic flow 
tables.!* 14 

The bulk of the heat-transfer testing was performed 
with a segmented copper model, and check measure 
ments were made with three thin-walled nickel models. 
The copper model (Fig. 1) consisted of a shell, 2 in. in 
outer diameter and 1/4 in. in wall thickness segmented 
into rings, the adjacent rings being connected by thin 
copper ribs and the remaining space between rings being 
filled with a thermally insulating plastic. A detailed 
description of this model, including the method of 
fabrication, can be found in reference 5. 

Model wall temperature measurements were made 
for 17. = 1.90, 2.80, 3.26, 4.20, and 4.87 at six different 
stations along the hemisphere of the copper model using 
36 A.W.G. copper-constantan thermocouples connected 
to General Electric photoelectric recording potentiom- 
eters and using the temperature of melting ice as refer- 
ence. The model was initially preheated or precooled 
to various uniform temperature levels, and a temper- 
ature-time record was obtained during the blow for 
individual model segments. It should be noted that 
negligible radial temperature gradients existed in the 
model because of the high thermal conductivity of 
copper. 

The second type of heat-transfer model consisted 
of thin continuous nickel shells (Fig. 2), which would 
not be subject to the uncertainties of segmentation and 
the presence of the plastic of the previous model. 
Three nickel models were tested in order to determine 
the effects of wall thickness, thermocouple wire diam- 
eter, and thermocouple material on the test results 
(see Table 1). 

The nickel models were tested only at a Mach Number 
of 3.21 using the same test procedure as in the case of 
the copper model. Details of the model and _ the 
method of its fabrication are described in reference 15. 

If the only energy exchange present is between the 
air stream and a model segment, the heat-transfer rate 
can be determined from the following: 


q = (m/S) c(dT, /dt) (2) 
where 

q = heat-transfer rate per unit time per unit 
surface area 

mS = ratio of weight to surface area of the model 
segments 

dT,/dt = instantaneous change of surface temper- 
ature with time 

c = specific heat of model material evaluated 


at temperature 7’, 


(In the case of the nickel model, the value of mS at a 
point on the hemisphere is equal to w-d, where w is the 
specific weight of nickel and d is the thickness of the 
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TABLE | 


Thermocouple Wit 


Approximate Diameter 


Configuration Wall Thickness in Material 
a 0.01 0.0050 Copper-Constantar 
b 0.01 0.0018 Chromel-Alumel 
c 0.08 0.0018 Chromel-Alumel 
nickel shell.) Since each segment changed temper- 


ature at a different rate, an isothermal temperature 
distribution existed only at the start of the blow. This 
is illustrated by Fig. 3 (the lowest Mach Number case 
was selected for illustration since the temperature 
changes were greatest). A nonuniform temperature dis- 
tribution would result in internal heat conduction 
effects invalidating the assumptions for Eq. (2). Fur- 
thermore, the convective heat-transfer rates for non- 
uniform surface temperatures are different from thos¢ 
for isothermal surfaces. Therefore, the only time that 
Eq. (2) may be applied directly is at the start of the 
blow, which means that the initial slope of the temper- 
ature-time record has to be determined. Because of 
the problems of recorder inertia and starting friction 
the direct measurement of initial slope on the poten- 
tiometer record was considered unreliable. Therefore 
slopes were measured at various time intervals after 
the start of the blow, and then the measured slopes 
were plotted against time and extrapolated to zero 
time. 

As pointed out in the theoretical considerations, 
“effective’’ temperatures and adiabatic wall tempera- 
tures are not greatly different, and the distinction be- 
tween the two temperatures need be considered only 
when (7, — 7,) is small. Such a situation actually 
existed for the present tests, since the maximum and 
minimum limits of surface temperature were +50°C. 
for a supply air temperature 7) of about 15°C. and 
resulted in considerable scatter of heat-transfer values 
when adiabatic wall temperatures were used in the 
analysis of the data. Consequently, a method was 
adopted for determining local heat-transfer coefficients 
and “‘effective’’ temperatures similar to the method 
used by Slack'® for determining adiabatic wall tem- 
peratures on a flat plate. The NOL procedure is out- 
lined in reference 17. Basically, the method consists 
of plotting heat-transfer rates versus body-surface 
temperature, the intersection of the resulting data curve 
and the line of zero heat-transfer rate giving the adi- 
abatic wall temperature or “effective’’ temperature 
whichever is pertinent. The heat-transfer coefficient 
is then calculated from its definition in terms of temper- 
ature difference and heat-transfer rate. The same 
method was employed by Stalder and Nielsen‘ for de- 
termining average values of heat-transfer coefficients 
on a hemisphere. 

The heat-transfer rates shown in Fig. + were deter- 
mined by applying Eq. (2) to the zero time slope ol 
the temperature-time record of the model segment, and 
these heat-transfer rates were plotted versus the corre- 
sponding surface temperatures, each point on the curve 
in Fig. 4 being evaluated from a completely independ- 
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ent test. Since all such curves seem to be straight 


lines, it was necessary only to determine the slopes of 


the lines in order to evaluate the heat-transfer coeffi- 
cients for each angular position and to determine the 
intercepts with the abscissa in order to obtain the value 
of the ‘effective’ temperatures. The straight line sig- 
nifies that the local heat-transfer coefficient is inde- 
pendent of temperature in the temperature range in- 
vestigated (+50°C.). 


RESULTS AND DISCUSSION 
Sibulkin'*® has shown that the heat transfer near the 
stagnation point of a body of revolution in incom- 
pressible laminar flow with constant fluid properties 
is given by 
Nup — 0.763 (BD? yp)" Pr.4 (3) 
where 
GB is defined as (dl dx) ,-0 


For the case of incompressible potential flow about a 
sphere 


B= 3U./D (4) 


This gradient of local velocity, 8, evaluated at the stag- 
nation point is one of the significant parameters which 
can be determined from pressure measurements in the 
study of heat transfer to a sphere. 

Since it was believed that the supersonic data would 
also satisfy Eq. (3), it was deemed desirable to find the 
influence of Mach Number on 8. It was a simple mat- 
ter to evaluate 6 directly by measuring the slope of lL 
versus x curves determined from pressure measure- 
ments for both supersonic and subsonic Mach Num- 
bers, since U; varies linearly with x. near the stagnation 
point. The values of 6 for all cases were combined 
with the body diameter and the stream velocity at 
infinity into a dimensionless parameter (8D), U 
which is plotted against undisturbed free-stream 
Mach Number, J/,, in Fig. 5.* It should be noted 
that for Mach Numbers greater than 4 it is possible to 
calculate 8 with a fair degree of accuracy from the prin- 
ciples of Newtonian flow.'’ Fig. 6 represents the dis- 
tribution of the nondimensional velocity, U\/@8D, at 
the outer edge of the boundary layer along the hemi- 
sphere for supersonic flow in the undisturbed stream 
plotted as a function of angular position. The points 
on the figure are not actual data values but are deter- 
mined from curves drawn through the pressure data 


* It is interesting to note the influence of the cylindrical after- 
body on the subsonic data for the hemisphere-cylinder. The 
data seem to agree more with potential flow around a Blasius 
half-body (solution for a point source in a uniform stream) than 
with flow around a sphere. The half-body was selected for com- 
parison at Mach Number zero since the problem of the hemi- 
sphere-cylinder in potential flow has never been solved mathe- 
matically, to the knowledge of the authors, because of the dis- 
continuity in curvature at the junction of hemisphere and cylin- 
der. It should be pointed out that the diameter used for the 
correlation of the half-body value was twice the radius of curva- 
ture at the stagnation point. 
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Fic. 7. Nondimensional local heat-transfer parameter ys 
angular position 


and are presented merely as a means of distinguishing 
the various Mach Numbers. It is significant to note 
that incorporating 6 (evaluated only at the stagnation 
point) into the nondimensional velocity results in a 
single curve for the velocity distribution in the free- 
stream Mach Number range of 1.22 to 4.87. 

All the heat-transfer data for the copper model are 
plotted in Fig. 7. The coordinates for the plot are 
Nup (8D? v)'* vs. angular position measured from the 
stagnation point, all fluid properties being based on 
local conditions at the outer edge of the boundary 
layer. The outer edge of the boundary layer is a con- 
venient location for the evaluation of fluid properties 
because the parameter Nu,» (8D? v)"? is approxi- 
mately independent of temperature for small rates of 
heat transfer and moderate Mach Numbers (the local 
Mach Number on the hemispheres never exceeded a 
value of about 2.5). 

As mentioned in the section titled Theoretical Con- 
siderations, compressible laminar heat-transfer data 
may be compared with incompressible theory for small 
rates of heat transfer and Mach Numbers up to 4, and 
these requirements were met by the current measure- 
ments on the hemisphere. The dotted line in Fig. 7 1s 
an approximate calculation of laminar heat transfer 
to spheres in incompressible flow with constant fluid 
properties and Pr = 0.7, using a method outlined by 
Drake‘ and a velocity distribution determined from the 
pressure measurements for incompressible flow by 
Fage.” An examination of the velocity curves of 
Fig. 6 indicates that the velocity distributions are 
quite different for compressible and incompressible 
cases. Since the theoretical calculation is sensitive to 
the nature of the velocity distribution at the outer 
edge of the boundary layer (1 U’., as a function of a), 
it was decided to recalculate the incompressible curve 
using a velocity distribution which better approximates 
the supersonic case. Consequently, an equation was 
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Fic. 8. Comparison of thin-walled model data with theory. 


fitted to the nondimensional, supersonic velocity curve 
in Fig. 6, and a hypothetical velocity distribution was 
determined. The heavy solid line of Fig. 7 is the heat- 
transfer curve for the new velocity distribution. All 
the supersonic data in the entire Mach Number range 
investigated have been compared with this single 
theoretical curve. 

The thin solid line is a cross-plot of the Stine and 
Wanlass® theoretical curve for a Mach Number of 1.97. 
It is presented here for comparison since it is the only 
other prediction of local heat transfer to a hemisphere 
in supersonic flow known to the authors. Most of the 
NOL data and the NOL theoretical calculation (heavy 
line) lie above the Stine and Wanlass® prediction and 
data. As mentioned in the introduction, the Stine 
and Wanlass® calculation is based on the Falkner- 
Skan’ wedge flow model. Eckert*!' explains that the 
Falkner-Skan’ model gives results which become more 
inaccurate as the distance from the body stagnation 
point increases. This is demonstrated by Fig. 11 of 
reference 10 which shows that the heat-transfer co- 
efficient predicted by the Falkner-Skan’ method for a 
cylinder in incompressible flow is about 15 per cent 
lower than experiment at the 75° position. It is im- 
probable that the calculation for heat transfer to the 
hemisphere in supersonic flow is not subject to the same 
kind of inaccuracy. 

The data obtained with the copper model were sub- 
ject to a scatter of up to 30 per cent, and many of the 
data were higher than the theoretical prediction. For 
this reason it was decided to make check measurements 
with the nickel models, which were considerably differ- 
ent in construction from the copper model and were not 
subject to the uncertainties of segmentation. The 
check was made at a Mach Number of 3.21 at which 
the earlier measurements deviated greatly from theory, 
and the results are presented in Fig. 8. No appreciable 
differences due to wall thickness or thermocouple size 
and material were detected in the data obtained with 
the 3 thin-walled nickel models, and so only the average 
values need be presented.* It can be seen that the 

* It should be noted that the 49° and 82° data points for the 
nickel model are not average values but were obtained on one 


model only (configuration a). 
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newer measurements are in reasonably good agreement 
with the theoretical prediction. 

Fig. 9 is a plot of all the ‘effective’ temperatures 
determined by the method of Fig. 4. In addition to 
the data points, two curves are plotted for the adiabatic 
wall temperature, the solid line representing the case 
where the recovery factor 7 is equal to Pr’? based on 
wall conditions and the dotted line for the case where 
the recovery factor equals 0.95 Pr’. It has already 
been pointed out that adiabatic wall temperature and 
“effective” temperature are not greatly different, and 
this fact is verified by the experimental results. The 
limitations of the data reduction procedure (extrapola 
tion of slope readings to time zero), however, make it 
impossible to draw any accurate quantitative con 
clusions about the data. An error as small as 5 per 
cent could completely change the entire pattern pic- 
tured in Fig. 9. It appears from the data, however, 
that the ‘‘effective’’ temperature is somewhat lower 
than the adiabatic wall temperature, a result which 
agrees qualitatively with reference 10. More detailed 
descriptions of the NOL tests, procedures, and results 
may be found in Navord Reports 3841*° and 3980" by 


the authors of this paper. 


CONCLUSIONS 


Within the accuracy of the data, the following prin 
cipal conclusions may be drawn from the current in 
vestigation on local heat transfer to a hemisphere in 
supersonic laminar flow with small rates of heat transfer 
for the Mach Number range of 1.90 to 4.87, isothermal 
surface temperature range of +50°C., and atmospheric 
supply conditions: 

(1) Asingle curve derived from incompressible theory 
and the actual velocity distributions along the surface 
can be used to represent nondimensional heat-transfer 
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coeflicients. Local fluid properties at the outer edge 
of the boundary layer are used in calculating heat- 
transfer rates from the incompressible correlation. 

(2) The local velocity distribution along the surface 
of a sphere may be expressed in terms of a nondimen- 
sional parameter, l/8D, which eliminates Mach 
Number dependency in the investigated range. 

(3) The “effective” temperature does not differ by 
more than 5 per cent from the adiabatic wall temper 
ature, but it appears to be somewhat lower than 
adiabatic wall temperature. Because of the small 
difference, the designer of high-speed vehicles subject 
to high rates of heat transfer need not consider the 
distinction between the “effective” and adiabatic wall 
temperature. 

(4) The local heat-transfer coefficient is independent 


of surface temperature. 
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Travel Arrangements 


Institute of the Aeronautical Sctences and The Royal Aeronautical Society Conference, 
London—Folkestone, September 1-16, 1957 


Arrangements have been made with “‘The Travel Center of Manhasset’’ to accomplish all travel and hotel reservations 
for the delegates and their wives. This travel bureau handles all arrangements on a ‘‘personal”’ basis. The bureau has now 
obtained commitments from air lines, steamship companies, and hotels to assure the availability of accommodations. The 
air-line commitments include group arrangements on TWA and Pan American flights which will allow delegates to travel 


Special arrangements also have been made by the travel bureau with American Express, Thomas Cook & Son, and Ask 
Mr. Foster for preconvention and postconvention tours of Europe by the delegates—individually or in groups. 

A travel brochure for the conference is being prepared and will be circulated to all members early this year. In the 
meantime, full travel information can be obtained and reservations made by addressing inquiries to: 


IAS Conference Travel, 2 East 64th Street, New York 21, N.Y. 
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Inviscid Hypersonic Flow Over Blunt-Nosed 


Slender Bodies’ 


LESTER LEES* ann TOSHI KUBOTA** 
California Institute of Technology 


SUMMAR\ 


At personic speeds the drag/area of a blunt nose 1s much 
ger than the drag/area of a slender afterbodyv, and the energy 


contained in the flow field in a plane at right angles to the flight 


is nearly constant over a downstream distance many 


direction 
greater than the characteristic nose dimension Phe 

rse flow field exhibits certain similarity properties directly 
inalogous to the flow similarity behind an 
found by G. I. Taylor, S. C. Lin, and A. Sakurai 
experiments of Hammitt, Vas, and 


13 in helium shows 


intense blast wave 

A comparison 
Bogdonoff on i 
te with a blunt leading edge at 
iccurately by this 


it the shock-wave shape 1s predicted very 


imilarity analysis he predicted surface pressure distribution 


Experimental results on a hemi 
VW = 7.7 in the GALCIT 


ite that not only the shock-wave shape but also the 


is somewhat less satisfactory 


air tunnel 


cylinder obtained at 
surlace 


pressures for this body are given very closely by the similarity 


eory, except near the hemisphere-cylinder junctio1 
considerations combined with a detailed study of the 


Energ\ 


equations of motion show that flow similarity is also possible for 


class of bodies of the form 7, ~ x”, provided that n m <1, 
where mm’ = 2/3 for a two-dimensional body and m’ = 1/2 fora 
body of revolution When 7 m’ the shock shape is not 
imilar to the body shape, and the entire flow field some dis 


rom the nose must depend to some extent on the details 
nose geometry 
finds 


By again utilizing energy and drag considerations one 


that at hypersonic speeds the inviscid surface pressures gener 
ted by a blunt leading edge are larger than the pressures in- 


duced by boundary-laver growth on an insulated flat surface 


over a distance from the nose of order /, where 

l/d = 2(Reqg/M 
For an insulated blunt-nosed slender body of revolution the 
corresponding distance is given by 1/d = 2(Reqg/ AM Here 


Reynolds Number based on leading-edge 
thickness, or nose diameter In free flight these constants are 
1,700 and 20, 


action effects are important over the forward portion of a blunt 


Reg is free-stream 


replaced by respectively, so that viscous inter 


nosed slender body only for relatively low values of Re How 


ever, ‘‘far downstream” of the nose the inviscid over-pressure is 


} 


small and viscous interaction phenomena will have to be taken 


Into account 
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SYMBOLS 
Free-stream quantities are denoted by thi ibseript 
while the subseript * denotes quantities ev te the bor 
surtact 
= sound speed 
1 constant 
( Chapman-Rubesin factor in relation p/p ( 
Cp nose drag coefficient, D 1/2)¢ rd } 
( pressure coetlicient, 1/2 
nose diameter or leading-edge thickne 
D = drag 
E energy in transverse flow field 
I p(z)/p 
geometric index 
influence lengt 
body length 
” exponent, /p ~ A 
VW = Mach Number, u 
pressure 
0” iny physical quantit 
distance normal to bod xis or chord line Xi 
R distance of shock wave from x-axis 
Re = Reynolds Number, ; d/ 
= distance along body surface measured from forwar« 
stagnation point 
t time 
T = absolute temperature 
u, 7 = velocity components parallel and normal t 
t = shock velocity in direction normal to x-axi 
U.(dR/dx 
\ = distance along body axis or chord line, m« red 
from forward stagnation point 
= R 
a = exponent, (1 m)/n 
ratio of specific heats, ¢ ( 
6 = angle between tangent to the body surface and thi 
¥-aXIs 
6 shock angle with respect to x-axis 
Ll = absolute viscosity 
p = density 


stream function 
|) INTRODUCTION 
\ \ J HEN A FINITE AMOUNT of energy is suddenly re 
leased at some “‘point’’ in a gas initially at rest, 
the radius of the intense spherical blast wave generated 
by the explosion grows like 
R Fi(y) (E/p F 


The flow field in the wake 
of the shock wave exhibits a certain similitude, in the 


as shown by G. I. Taylor. 


sense that the pressure, density, and outflow velocity 
are described by relations of the form 
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Q(r)/Q(R) = f(r/R) 


This similarity holds only in the intermediate zone not 
too close to the origin of the explosion, where the theory 
predicts that 7 —~ © and p — 0, yet not so far away 
that the shock strength has decayed to a level where 
the strong shock approximations are no longer ap- 
plicable. Taylor's analysis was later extended to the 
case of a cylindrical blast wave by S. C. Lin,? who 
found that 
R = Fi(y) (E/p,)”4 t”? 


in this case. Lin also remarked that, according to 
Hayes’*® concept of hypersonic similitude, this relation 
for R(t) should describe the shape of the bow shock 
wave behind an unyawed, axially symmetric body 
traveling at a uniform hypersonic velocity. The axial 
flow velocity with respect to coordinates fixed in space 
is very small, provided that the shock angle 64, is not 
too large, and the flow in a transverse plane fixed in 
space behind the body resembles the flow generated 
by the explosion of a long, highly concentrated cylin- 
drical charge at the time ¢ = 0. Here t—x«/U,, and 
the energy F/ per unit length of charge is identified 
with the total drag of the body. 

In reference 4, one of the present authors pointed 
out that these considerations are equally applicable to 
the shock wave generated by a blunt nose of finite 
radius on an unyawed slender body.* At hyper- 
sonic speeds the drag of the nose per unit cross-sectional 
area is much larger than the drag/area of an afterbody 
with a uniformly small slope in the meridian plane. 
To be specific, the drag of a blunt nose of diameter d 
(or leading edge of thickness d, for a two-dimensional 


body) is given by 
Dy ~ (1/2)p,,U,? d**! 


while the drag of a conical (or wedge-like) afterbody 
of half-angle @ and length LZ, for example, is ~ 


(1/2) p., U2 (207) (@L)**! 


where k = 0 for a two-dimensional body and k = 1 for 
a body of revolution. The drag of the afterbody be- 
comes comparable with the nose drag only when 


L/d ~ 1/9&+%/4+) 


In other words, the shape of the bow shock wave, the 
inviscid flow field, and the surface pressure distribution 
on a slender body are dominated by the blunt nose or 
leading edge over a downstream distance many times 
greater than the characteristic nose dimension. The 
analogy with a constant-energy, nonsteady similar 
flow of the type investigated by Taylor and Lin is com- 
plete for the particular case of a blunt nose followed by 
a cylindrical afterbody (@ = 0). In this case the shock 
shape is described by 

R ~ [(x/8)p,, U..*d?Coy/p..]"* (x/U,)”? 
or 


* Since reference 4 received only limited distribution the 


main points of the discussion are repeated here. 
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R d = Ki(y) (Ca) its d)/? 


The analogy is readily extended to two-dimensional 
bodies,* where it is complete for the case of a flat su; 
face with a blunt leading edge at zero angle of attack. 
such as the case studied experimentally by Hammitt, 
Vas, and Bogdonoff.* © (Here the energy in the trans 
verse flow field over one surface corresponds to one-hali 
of the nose drag.) This particular analogy was also 
noticed independently by Cheng and Pallone.’ Dr. 
S. C. Lin** pointed out to the present authors that 
Sakurai* has extended Taylor's constant-energy analy 
sis to the plane blast wave problem, where the shock 
Fi3(y) (E/p t? 


In a second paper,’ Sakurai also treats not-so-strong 


wave shape is of the form R 


shocks (later times) by constructing asymptotic series 
solutions in inverse powers of (a@,/v,)*, or shock 
strength. This technique can be readily taken over 
to the problem of hypersonic steady flow over blunt 
nosed slender bodies [see Section (4b) ]. 

These considerations suggest that it would also be 
worth while to investigate the more general case in 
which the energy of the transverse flow is increasing 
with distance from the nose, but the shape of the body 
is such that flow similarity is preserved. The corre- 
sponding nonsteady flow problems are the expanding 
sphere, expanding cylinder, and motion of a piston in 


a long straight tube. 
(2) SmmILAR FLOWS: ENERGY AND DRAG 
CONSIDERATIONS 


Taylor's assumption of flow similarity in a fixed 
transverse plane is satisfied only for ‘‘strong’’ shocks, 


where 
v(R) = 22,/(y + 1) 
P(R)/p,, = 2 95?/(y + 1) 
and 
p(R)/p, = (y + 1)/(y — 1) 

The strong shock approximation in turn is applicable 
only when 

ly — 1)/2] (e.?/a,2) = My — 1)/2] M202 > 1 Tf 
In addition, 


v(R) ~ R-“ or R*(dR/dt) A(const) 


and 
R = (A/U,m)™x™ 
where a = (1 — m)/m 
Also, the boundary condition v(7) U.(dr,/dx) on 
the body requires that v(7)) = 22;, or 3 = 1%/R = 


** Now with Avco Research Laboratories 
t Previous experience with hypersonic similarity suggests that 
this approximation is useful when 


((y — 1)/2] M,20,2> 2-3 
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INVISCECID 
const, if flow similarity is to exist—1.e., 7, ~ x”——and 
the shock and body are similar. 


When these conditions are satisfied, we may write 


o(r, t A (z)/R* b(r, t)/p.. = A? F(z)/R™ 


r/R (see Fig. 1), and the energy E associated 


and pl’, t)/p 
where 2 
with the flow field in a transverse plane is expressed as 
follows :* 


ek 


E ey aa | p[c,T + (v?/2)] r* dr (1 
J ks 
or 
71 
Eats" s 2° a HF/(y¥ -— 1] + 
(1/2)p?}2* dz (la) 
where k k’ 0 for planar flow, k k’ 1 for 
axially symmetric flow, and k’ = 1 and k 2 for 
nonsteady spherical flow. An energy balance shows 
that 
dE/dt = U,(dE/dx) = 2* x*’ ry* po vp (2 


in other words, the rate of change of energy of the fluid 
motion is equal to the rate at which work is done by 
the pressure forces acting on the fluid along the body 
surface. Evidently, from Eq. (la), £ const. when 
2a 1+ k, orm = 2/(3 + k), and by Eq. (2) » 

Zp () everywhere, except right at the nose. For 
spherical flow a = 3/2, m = 2/5, and R ~ t*® (Tay- 
lor'); for axially symmetric flow a = 1, m 1/2, and 
R~x 1/2; m = 2/3, 
and R ~ ?#/3 (Sakurai’), or R ~ x/* (Cheng and 
Pallone’ and Lees‘). On the other hand, when 2a < 
1+ k, orm > 2/(3 + k), then dE /dt > 0, and w% > 0, 
z, > 0. For a positive body slope (or an expanding 
sphere, cylinder, or piston), similar solutions exist (1 


(Lin?); for planar flow a = 


; ae dian + na 
at all) only form > 2/(3 + k). The same conclu- 
sion is reached by considering the pressure drag. For 


these bodtes, 


+) >R 
D = 2 “| ry Py (dr,/ds) ds = const Re -ak 
or 
‘i 
D = const ; tai t 


the drag is finite only when m > 2/(3 + k). 


The quantity E has the dimensions of energy /area for planar 
flows, energy/length for axially symmetric (cylindrical) flows, 
ind the energy itself for nonsteady spherical flows 

+ Stewartson (private communication) also found the re 
striction m > 2/3 in his study of boundary-layer — shock-wave 
interaction over a planar body of shape ry ~ x”. 
the lateral velocity given by the inviscid solution at the outer 


2 


For m < 2/3 


edge of the boundary layer is negative, and he was unable to 
match it with the positive (outward) lateral velocity given by 
the sum of the body slope and the boundary-layer growth 
However, no explanation was offered for this behavior, and the 
special significance of the case m = 2/3 was not explored. 
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There is no smooth connection between the class of 


solutions for bodies of the form 7, ~ x” and the special 
constant-energy solution with m 2/(3 + k). This 
point will become clearer when the mathematical 
properties of these solutions are discussed (see Section 
3). The similar-flow solution yields a shock wave 
that lies farther and farther away from the body sur 
face as m — 2/(3 + k) and zs, — O, and the slender 
body approximation must become progressively poorer 
as this critical value of m is approached. 

The criterion m 2 2/(3 + k) is a necessary condi 
tion for the existence of a similar solution, but a study 
of the mathematical properties of the equations of mo- 
tion is required to determine whether this condition 
is also sufficient. 


(3) SIMILAR FLows: MATHEMATICAL CONSIDERATIONS 


Within the framework ot the hypersonic slender-body 
theory* the equations of motion for similar flows have 
the same form as the blast-wave equations’ * * ° 
except that a now takes on any value less than or equal 
to (1 + k)/2, rather than the special values 1/2, 1, and 
3/2. The transverse momentum equation becomes 


(b — 2)%’ — ab+ F’/y 0 (3 
the continuity equation yields the relation 
(@ — z)p’ + YP’ + ky b/s = 0 (4) 


and the condition of constant entropy following a fluid 
particle leads to the equation 


(@ — 2) [F’ —_ (yFy’ vy) | — Yak 0 (2) 


where the prime denotes differentiation with respect 


to zs. The boundary conditions are as follows: on 
the shock wave, F(1) ?(1) 2/(y + 1), and 
v(1) = (vy + 1)/(y — 1); on the body surface ®, 


2, when m > 2/(3 + 
z — 0 for the constant-energy solution with m 


2), ora < (1 + k)/2; ?>O as 


(3 + k). 
By multiplying Eq. (4) by 2a/¥(® — sz) (1 + &), 
dividing Eq. (5) by F(® — 2), and adding the two 


results one obtains the equation 


F'/F — {y — [2a/(1 + &)]} (W’/¥) 4 
[2a/(1 + k)] [(@ — 2)’ 
[2ak/(1 + k)] (1/z) = 0 


(db — sg 


) 
“ye 
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which admits of the integral* 
F(z aa @:p) (2 1+k)] a J 
gltak/ (+) const (6) 
pir Pa/+h)] 
From Eq. (6), we see that y > (2 — ©)"as#—>2—> %, 
where 
1 2a [y(1 + k _ 2a] 
2(1 — m)/{m[yQ1 + Rk) + 2] — 25 


Thus ” must be positive if the density is to be finite at 


the body surface, or 


In fact, p > 0 as s > &, exactly as in Taylor’s case, 
unless m takes on one of these extreme values, in which 
case p — const. But for any real gas the lower limit 
is always less than or equal to the value 2/(3 + k) 
imposed by energy and drag considerations (see Section 


2), so that actually m is confined to the range 


dt (yt — 2a) {((y +1) 0 —s)l2 + [(2 

d a(s 1) [i 

According to Eq. (7), = Ga/(y + 1) at the shock and 

t 2a/y at the body surface for k Q. Hence, in 
the s-t plane the body surface is represented by the 
point (s 1, t 2a/y), and the shock by the point 


(s 2/(y + 1), ¢ Ga/(y + 1)]. There are nine 


singular points of Eq. (8), but only four of these are 


relevant to the present discussion. These four singular 


points are as follows (see Fig. 2a): (1) s = 1,t = 2a/y, 
a node (body surface, point B); (2) s i 2a/y, a 
saddle point (point A); (3) s ae. (1 + 2a)/¥ 
(point C), a node when a > (y — 1)/2, or a saddle 
point when a < (y — 1)/2; (4) s f= 201 + @)3 
(y + 1) a saddle point when a > (y — 1)/2, or a node 
when a (y — 1)/2. Inthe regions S$ 1,t 2 2a/y 


the slope dt/ds is negative below the line ¢ s, while 
above this line dt/ds is positive, provided that ¢ S f,(s), 
where ¢;(s) is the upper branch of the curve defined by 


the relation 


(y¥t)HNd-sy@+ fv +3)a4+2j) 5 —- 
(3a + 2)]t — a(2a +1) 5 0 


Along the line / = s the slope dt/ds is infinite, except at 


the saddle point s = ¢t = 2a/y. Near s l,¢ = 
2a/y the integral curves are of the form t — 2a/y7 
C(l — s)”%~*”, and hence dt/ds ~ 0 as s—> 1 for 


0<a< ¥/2. 

So long as 6a/(y + 1) S 2/(y + 1),0r0 Sa £ 1/3 
and 3/4 S$ m S 1, thent < s. In this case, the in- 
tegral curves connecting the shock wave and the body 
surface need never cross the line ¢ s, and an analytic 
solution certainly exists (see Figs. 2a and 2b). When 
1/3 < a < 1/2, the integral curve must cross the line 


t = s in order to reach the point in the ¢ — s plane 


* A first integral of this form is also obtained by Sakurai® and 
Latter” for the special case of constant energy. 
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for a body with a positive slope. 
» 


Since the basic equations, Eqs. (3)—(5), are invariant 
under the transformation @ —~ a®, z ~ az, Ff — q?f 
y — y, these equations are reducible to a single non 
linear first-order equation in the new variables {¢ 

d@/dz and s P/s. 
first eliminating F’ between Eqs. (3) and (5) to obtain 


This reduction is accomplished by 


the relation 


The quantity y/F is then eliminated by differentiating 
Eq. (7), and utilizing Eqs. (4) and (5). By employing 
the identity 


2(d/dz | ies (d/ds 


the following equation is finally obtained for two 


dimensional flows :+ 


> > ‘ 
Va “7 ja) s — (Sa +? al(2a 1) sy 
‘\ 
A A - ff 
Y é 


corresponding to the shock wave. Clearly, all of the 
integral curves that lie above the dividing curve pass 
ing through the point ¢ s 2/(y + 1) are ‘“‘turned 
back”’ at the line ¢ 
Ya). This fact led the senior author t 


s, and never reach the shock point 
(see Fig. 
the erroneous conclusion’ that no similar solutions 
exist in the range 1/3 < a < 1/2 
ent analysis shows that every integral curve contained 


However, the pres- 


in the triangular region ABC (Fig. 2a) is bound to 


reach the body-surface point B. Since the shock point 


s 2/(y + 1), ¢ Ga/(y + 1) hes within this region 
so long as a 1/2, analytic solutions also exist in the 
s-t plane for 1/3 < a < 1/2, and the corresponding 
curves cross the line ¢ s below the point ¢ 5 

2/(y + 1) (see Figs. 2a and 2b). The fact that 


dt/ds > © at the point ¢ = s does not lead to any singu 


larity in the physical plane. 

To summarize: similar solutions exist in the entire 
range 0 < a < 1/2, or 2/3 < m € 1, for two-dimen- 
sional flows (see Fig. 2c), and in the entire range 0 < 
a <1, or 1/2 < m € 1 for axially symmetric flows 


From the present point of view the constant-energy 
solution m = 2/(3 + ) is a singular solution passing 

t The restriction m < 1 may also be derived by considering 
the entropy distribution imposed by the shock-wave shape 
Along the shock the quantity p/p? is proportional to p, or p/p’? ~ 
MK) so that p/p? ~ w?(" 


("™—) while vw ~ R*t! ~ x 
But the quantity p/p? is constant along streamlines, and p ~ 
warm /ymk+) pV throughout the flow field. In particular 
near the surface WY — 0 and p — O (or a constant) if m 
and p— ~ifm> 1. Actually, the domain m > 1 is also ruled 
out by the fact that the shock wave at the nose or leading edge 


l, 


is of zero strength, while the entire analysis is based on the 
strong-shock approximation. On physical grounds one would 
not expect the flow over concave bodies with m > 1 to exhibit 
similarity properties. 

t The analysis for axially symmetric flows is similar, but the 


algebra is more extensive. 
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INVISCID 
through the saddle-point (see Fig. 2b), with a slope 
ip/dz l/yats 1/¥ 


Included within the admissible range of values of m 


l 


spe Fy 2% 
(See Fig. aC). 


wre of course the wedge and cone (m 1), and also the 
or body of 


hypersonic optimum shape” r, ~ x 
revolution of minimum zero-lift drag for a given fineness 
ratio, as determined from Newtonian impact theory neg 
lecting centrifugal foree by Eggers, Dennis, and Resni- 
Cole 


koft By including centrifugal force, J. D. 
obtained the value m 2/3 for this optimum shape. 


For two-dimensional flow Cole obtains an optimum 
shapé with m 0.87; both of his cases also lie within 
the range 2/(8 +k) <<m S31. 


When these similar nonsingular solutions do exist one 


expects them to provide a good approximation to the 
pressuré¢ and velocity fields not too close to the blunt 
nose, provided also that m is not too close to the critical 
3 +k 


for example) is given by 


value 2 The surface pressure distribution 


b() m* _ yM 
F(2,) + 
p (e/E 
or’ 
( »m?F(%,) r7/2 v/s 
where r } a 
For these bodies, the results obtained by utilizing 


iny one of the purely “‘local’’ hypersonic approxima 
tions, such as tangent-wedge (or cone), or Newtonian 
plus centrifugal force, are similar in form, which gives 
one some confidence in these approximations. When 
+ k), however, we conclude that the shock 


>) 9 
- ” 


shape is not similar to the body shape, and the entire 


ml « 


flow field some distance from the nose must depend to 
some extent on the details of the nose geometry. It re- 
mains to be seen whether any simple local hypersonic 
approximation is applicable to a blunt-nosed slender 
body in these cases. 

t) COMPARISON BETWEEN THEORY AND EXPERIMENT 
FOR THE CONSTANT-ENERGY SOLUTIONS 


(a) Flat Plate with a Blunt Leading Edge 


A clear test of the analogy between hypersonic flow 
over a blunt-nosed slender body and the constant- 
energy Taylor-type flow is provided by the experimental 
investigation carried out by Hammitt, Vas, and Bog- 
donoff® © on a flat plate in the Princeton helium tunnel. 
The blunt leading edge is formed by taking a plane cut 
normal to the upper plate surface, which is parallel to 
the oncoming flow. The lower surface is inclined at 
10° to the flow, but does not influence the upper sur- 
face. In these tests the Mach Number ranges from 
11.4 to 13.8, and the shock angles are such that the 
assumptions of the strong shock theory are fully satis- 
fied. Shock-wave shapes were determined from inter- 
ferograms over a range of leading-edge thicknesses 

“Here F(z,) and zp» are functions of m and y; their values are 
now being determined for a few cases of interest. 
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OW 


0.17 X. 10° m. sd Ss 59 X 10-* m., or 120 Sv Ke, & 


70.6 X 10%. For Re, > 16 X 10° viscous effects are 
negligible (see Section 4c below), and the empirical fit 
to the data presented in reference 5 is R 1.36 d 

x?) Our numerical results give R 1.44 d’* x7'* in 


this case, which is quite close to the observed shock 
shape (see Fig. 3).7 
According to the similarity theory the surface pres 


sure 1s given by the expression 


p/p 0.40 AM v/d 
As shown by Cheng and Pallone (Fig. 1 of refer 
ence 7) the predicted surface pressures are certainly 
of the correct order numerically However, the best 


empirical fit to all of the data offered by the experi 


menters themselves,” ° is 


Ap p 0.0161 V/ v/d 


in the range 4 X 10° <= Re, = 70.6 X 10 A close 
examination of the Fig. 7 of 


shows that the surface pressures are strongly influenced 


data (see reference 5 
by the sudden expansion and recompression around 
the 90° corner at the leading edge 

had 
tributions for a similar geometry in the 11 X 11-in 
heated air tunnel at the NACA Langley Laboratory 
at 


inverse 2/3 power dependence on x/d in 


Previously, Bertram’* measured pressure dis 


1,960 show an 
the 


6.86, and his data for Ré 
range 


10 « S/d 
of over-pressures and leading-edge thicknesses is not 


70 (approximately However, the range 
wide enough to permit any definite conclusions to be 
drawn. Bertram’s'* more recent measurements of sur 
face pressures and shock wave shapes for blunt-tipped 
cones and a flat-nosed cylinder also show the strong in 
fluence of the local expansion around the shoulder. 

7.7 in Air 


(b) Hemisphere-Cylinder at M 


The experimental results discussed in Section (4a 
suggest that the analogy with a constant-energy Tay- 
lor-type flow might be more complete for a blunt-nosed 
slender body with continuous slope, such as a hemi- 
spherical nose followed by a cylindrical afterbody. 
An experimental study of the shock-wave shape and 
surface pressure distribution for such a body has been 
carried out by the junior author in the GALCIT air- 
operated, continuous flow hypersonic wind tunnel at 
M.. = 7.7 (Leg No. 2 
700°F., thus ensuring one-phase flow. 


For most of the tests p 315 
psi and 7» 
The model is 0.75 in. in diameter and 10 in. long, and 
is supported from the rear by a sting 0.5 in. in diam- 
eter. The model is mounted in a portion of the test 
rhombus in which the static-pressure variation is +4.0 
per cent along the centerline. The Reynolds Number 
based on the diameter of the cylinder is 1.80 X 10°. 


= 


+ For the constant energy flows [m = 2/(3 +k 
authors obtain solutions of the equations of motion in 
form. This property first for the 
(Taylor) case by Latter,’ but it holds also for axially symmetric 


the present 
closed 
spherical 


was discov ered 


and two-dimensional flows 
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FIG. 4 SHOCK SHAPE FOR HEMISPHERE-CYLINDER AT MACH NO 7.7 
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For this case the shock-wave shape given by the 
strong-shock (Taylor-Lin) theory is found to be R/d = 
0.78 (x/d)' *, and the second approximation, obtained 
by Sakurai’'s method,* gives the result 


211 


R/ad 0.78 (x/d) /2{1 + 1.62 [(x/d)/M.?]; 


a 


In Fig. 4 these predicted shock shapes are compared 
with the shock shape determined from schlieren pho- 
tographs. Evidently the form of the bow shock is 
closely reproduced by the blast-wave analogy; in fact 
the experimental and calculated shock-wave ordi- 
nates differ by an almost-constant amount of about 
0.3 d. This outward displacement of the shock rep- 
resents the effect of the cylindrical afterbody of finite 
diameter. 

According to the similarity theory, the surface pres- 
sure on the hemisphere-cylinder is given by 


p/p. = 0.0665 M,?/(x/d) 
in first approximation, and by 


p/p 0.0665 [MW .?/(x/d)]{1 + 6.08 [(x/d)/.7_° 
(0.0665 M_?/(x/d)} + 0.405 


in a second approximation obtained by Sakurai’s 
method.’ The slow decay of the pressure over the cy- 
lindrical afterbody is predicted very closely by the 
second approximation to the blast-wave theory (see 
Fig. 5), while over the hemispherical nose the pres- 
sures are given quite accurately by the modified New- 
tonian approximation,’ except near the nose-cylinder 
junction.* In the transition region between these two 
flow regimes no analytical treatment is yet available. 
An attempt is made to predict the surface pressures by 
a modification of the “‘local-expansion’’ method, in 
which the Newtonian pressure distribution C, 

C, sin* @ and the Prandtl-Meyer relation 


(1/p) (dp/dé) = + (yM?2/V M? — 1) 


are matched at the point on the surface where both the 
pressure and pressure gradient given by these two 
formulas are equal. For J, = 7.7 and y = 1.4, the 
matching point on the hemisphere-cylinder occurs at 
§ = 33.8°, and J = 1.37. The surface pressures are 
predicted quite accurately by this approximation up 
to the shoulder (@ = O°) (see Fig. 5). 

These experimental results have encouraged us to 
carry out an extensive experimental study of shock- 
wave shapes, surface-pressure distributions and flow 
fields for bodies of the form 7 ~ x”, as well as the 
hemisphere-cylinder. This investigation is now in 


pre gTeSS. 


(c) Viscous Effects 


Most of the theoretical investigations of hypersonic 


boundary-layer -external-flow interaction tacitly assume 
that the leading edge or nose of the body is infinitely 


“Note that there are no empirical constants in the blast- 


Wave theory, and only the constant Cp, is introduced in the 


modified Newtonian approximation. 
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sharp. But at hypersonic speeds even slender bodies 
must be blunt-nosed to some extent because of the 
critical heat-transfer problem. The question natu 
rally arises as to the relative importance of the inviscid 
pressure field associated with a blunt leading edge 
or blunt nose), and the self-induced pressure generated 
by boundary-layer growth. An estimate of the rela 
tive magnitude of these two effects can be obtained by 
considering the energy introduced into the transverse 
flow field by the blunt leading edge and by the pressure 
drag, D,, associated with the ‘effective’ body shape. 
For example, the quantity D,. for a flat plate is given by 


Dy = 'p,u,? (M2 C¥4/(Re,)*] [(3po0)/y] 


according to the strong-interaction theory.'® This drag 
becomes comparable with one-half of the nose drag 
D,,/2 ~ pu? d/2 when 


I/d = (y/6pod0)* (Rea/CM,2)* = a(Rey/M.2)3 (9) 


! 


For an insulated surface under wind-tunnel conditions 
at W. = 15, a = 1.0 for helium and a = 1.9 for air. 
Thus, at JJ, = 15, the viscous interaction effects 
should be important for Reg < 450, but when Re, in- 
creases beyond this value, the inviscid effects spread 
rapidly over the plate surface. Certainly for Re, > 
1,000 the inviscid pressure field is dominant over the 
forward portion of the body and determines the bound 
ary-layer development, skin friction, and heat transfer 
in thisregion. These rough estimates are in qualitative 
agreement with the experimental results of references 
Sand 15. 

For bodies of revolution the corresponding distance 
from the blunt nose at which the viscous interaction 
drag is comparable with the nose drag is given by the 
expression 


l/d = b (Reg/ M,? (10 


where b = (y/6) [1/(Cpodo?) ]. 
drical afterbody, b = 2 for air and b 


For an insulated cylin- 
1.7 for helium 
under wind-tunnel conditions at .\/ 15. Here the 
viscous interaction effects are important for Re, 
1,000 at 17. = 15. As Reg increases, the inviscid ef 
fects spread downstream much more slowly than in 
the two-dimensional case, because of the more rapid 
decay of the inviscid pressure field, relative to the 
viscous pressures. However, for Re, > 10% the in 
viscid pressure field is dominant over the forward por- 
tion of the body. 

The situation is quite different in free flight, where 
the surface temperature is much lower than the re 
covery temperature, and the boundary layer is con- 
sequently much thinner than it is on an insulated sur- 
face. For example, for air at JJ. = 15, with 4 

.29 and a surface temperature five times higher than 
free-stream temperature, the constant factor in Eq. (9 
is about 1,700, and the constant factor in Eq. (10) is 
about 20. Thus, viscous effects are important over the 
forward portion of a slender, blunt-nosed body only 


t The constants fp and 6» are defined in reference 15 








202 OF THE 


JOURNAL 


when Re, 
when Re, 
Of course 
inviscid over-pressure is small and viscous interaction 
phenomena will have to be taken into account. 
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the Atlantic! We look forward to an ever-in- 


creasing importance of these possibilities of civil avia 


may 


tion which, coupled with the military nuclear deterrent, 
The 
spective effectiveness of a United Nation's police force, 


should prevent a catastrophic world war. pro- 
which at long last has gotten a start from the Near 
East troubles, may well be the answer. 

I agree with Sir Arnold's forecast of an ever-increas- 
ing volume of civil air transport. I am convinced that 
this will achieve what may truly be called mass air 
transportation in the not distant future. Certainly, 
the advance has been by geometric progression during 
the past years. My own forecast indicates that present 
United States domestic air transport, in terms of bil- 
lions of passenger-miles, will triple by 1965, and that 
United States international and world aviation, in the 
same terms, will extend the 1956 achievement by three- 
to fourfold. 
Arnold’s comments on improved efficiency in terms of 


I am fortified in these estimates by Sir 


break-even load factor of the coming generation of jets 


187 


which points toward the practicability of fare reduc- 


tions, which, in turn, have a big effect on increasing air 


transportation by reaching more and more _ persons 


C 


sapable of paying the tariff. This, I think, will follow 


a cube rule which states that the potential market for 


air transportation varies inversely as the cube of the 


fare. 


He also paints a rather rosy picture for the in- 


crease in air freight. 


I 


There are many quotable passages in Sir Arnold's 


yaper. I particularly liked the concluding one: ‘‘The 


military side of aviation is an essential guard for sta- 


bility; the civil side an important means of exploiting 


t 


he stability for the common good.” 


I think Sir Arnold Hall has given us much food for 


thought in his excellent paper, and I am sure I express 


t 
c 


a 


he views of this audience in indicating our great appre- 
iation for the time and effort he has spent in its prep- 
ration and for his coming to our country to present it 


to us today. 
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On the Flutter of Swept Wings 


PETER F. JORDAN* 
The Glenn L. Martin Company 


SUMMARY 


Tt} concept of a sweep parameter «x which combines sweep 


ispect ratio 1s developed on the basis of the streamwise 


rip method Charts of reference speeds varying sweep param 
wing density, and Mach Number are shown. The theo 
| possibility of bending flutter and its effect on the flutter 


wings and of wing-body combinations are considered 


il flutter ind improved analytical methods are discussed 


INTRODUCTION 


I THE DAYS of simple, straight wings, low speeds, and 
large aspect ratios, the number of nondimensional 
parameters that determined the flutter speed was rela 
tively small. The result of a flutter calculation, based 
on aerodynamic strip method, could almost be predicted 
from charts of standard results and could be relied 
upon to agree fairly well with the experimental flutter 
speed (except for some control-surface configurations, 
but control-surface flutter no longer seems to be the 
problem that it was). Today we have to deal with a 
number of additional parameters, with more compli- 
cated structures, and with much more complicated 
aerodynamics. It would seem hopeless to prepare the 
comprehensive charts that would allow us to predict 
the flutter speed of any given wing. 

This leads to a situation that is quite unsatisfactory. 
For a given wing at a given Mach Number we still can 
solve accurately, in theory at least, both the structural 
and the linearized aerodynamic problem. In practice, 
we have to simplify at least the aerodynainic problem. 
We then calculate a flutter speed which may or may not 
agree with the experimental speed. If it doesn’t, what 
is wrong? Have we simplified too much? If so, where? 
Does linearized aerodynamic theory fail in this particu- 
lar case, or is there an error in our arithmetic? Perhaps 
the experiment is inadequate. Another possibility, 
not unusual, is that both analysis and experiment are 
all right, but we hit on one of those configurations where 
a slight change in one of the parameters produces a 
large change in the flutter speed. How are we going 
to find out? And if calculated and experimental speeds 
agree, what conclusions can be drawn? Is it a chance 
result, or can it be generalized? If so, to what extent? 

The only way to find out is to gain experience, to col- 
lect this experience and to coordinate it; indeed, much 
useful work is being done in this direction by a number 
What we need in order to be more economical 
It is 


not sufficient, obviously, and indeed quite dangerous, 


of teams. 
in this effort is a common basis for coordination. 


Session, Twenty-Fourth An- 
1956 


Presented at the Aeroelasticity 
nual Meeting, IAS, New York, January 23-26, 
* Head, Dynamics Research Staff 


to condense results in statements like, ““On a bomber, a 
variation of this parameter by one per cent increases 
the flutter speed by x per cent. 

An accepted and useful basis is the reference speed 
calculated by means of incompressible strip aerody 
namics. However, our basis, in order to be most eco 
nomical, should be as diversified as can readily be 
handled, though it must necessarily remain fairly sim 
ple. The first aim of this paper is to point out that the 
reference speed can be extended to include wing sweep 
and Mach Number and still be covered by quite com 
prehensive charts of flutter speeds, which can be pre 
pared once and for all with relatively little effort 

The second part of this paper deals with a few addi 
tional points: the effect of the theoretical possibility 
of pure bending flutter of swept wings on wing body 
flutter; a confirmation of strip aerodynamics with a 
somewhat complicated configuration—namely, a T 
tail with swept fin and swept stabilizer; checking if two 
or three structural modes are sufficient for a given flut 
ter problem; and a new three-dimensional aerodynamic 
method. 


AXIS-NORMAL AND STREAMWISE APPROACHES 


In calculating the aerodynamic forces that act on a 
swept wing of very large aspect ratio only the axis 
normal component I’ cos A of the air stream need be 
taken into account. This leads to a simple rule: if a 
large aspect ratio cantilever wing is swept by rotating 


it about its root (see Fig. 1), its critical speed I” varies 
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Swept wing and unswept reference wing 


Fic. 1. 
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Fic. 2. Chordwise pressure distributions one chord length 
from wing root (swept rolling oscillation ). real, - imagi 
nary. Left: streamwise method. Right: difference between 
methods-——‘‘streamwise’’ minus ‘‘axis normal.” 


approximately as sec A: 


V(A) =~ V°(A) = V(O)sec A (1) 


According to this zero-order approximation, sweep 
would always increase the flutter speed. 

A logical refinement of this approach arises from the 
consideration that the oscillating cantilever wing ex- 
periences spanwise deformations, and that, in conse- 
quence, the spanwise component I’ sin A of the air 
stream also produces downwash. The full downwash 
on the wing is the sum of the two downwash compo- 
nents. Several authors have proposed to start from this 
full downwash but to apply to it only the axis-normal 
the air stream, therefore using axis- 
This first-order approach we here call 
It has been developed 


component of 
normal strips. 
the “‘axis-normal approach.” 
in particular in the well-known paper by Barmby, 
Cunningham, and Garrick.' Further refinements are 
indicated in the same paper but these refinements tend 
to become cumbersome. 

An obvious alternative would be to consider the full 
air stream IV and to apply it to streamwise strips (see 
Objections to this approach are justified only 
In the aero- 


Fig. 1). 
if it is not followed through consistently. 
dynamic strip method each strip is considered as part 
of a uniform infinite wing. Logically, on a swept wing 
each streamwise strip has to be considered as part of an 
infinite swept wing. It follows that the aerodynamic 
forces to be used on a streamwise strip which is part of 
a swept wing are cos A times the forces that would be 
calculated for the same strip, and for its identical mo- 
tion, were this strip part of an unswept wing. This 
approach, which we call the “‘streamwise approach,”’ 
leads to the same limit for the case of very large aspect 
ratio as the axis-normal approach, namely to the zero- 
order result, Eq. (1). The details of the streamwise 
approach are developed in reference 2. 

As an illustration of the difference between the two 
approaches, consider a swept wing that performs swept 
rolling oscillations with the rolling axis normal to the 
wing axis. The downwash is constant along any axis- 
normal strip of this wing but is linear in x along any 
streamwise strip. Thus the axis-normal flow compo- 
nent |’ cos A meets a constant downwash while the full 
flow V meets a linear downwash. Obviously, the two 
results cannot be identical. Fig. 2 shows the difference 


between the two pressure distributions along the 


AERONAUTICAL 


SCIENCES MARCH, 1957 

streamwise strip at the distance of one chord lengt} 
The imaginary part of this dij 
On the other hand, the rela 


from the wing root. 
ference is not negligible. 
tive importance of this difference is inverse to the dis 
tance from the root. 

It would be somewhat futile to examine the differ 
ence between the two methods closely; they lead t 
practically identical flutter speeds if the aspect ratio js 
so large that a strip method is fully justified. Th, 
streamwise approach has a two-fold appeal; it js 
complete in that it does not neglect any part of the flow 
and, further, it does not involve the practical difficul 
ties that arise for the axis-normal approach at wing root 
and wing tip where some axis-normal strips end at the 
side-boundaries of the wing. A complication arises for 
the streamwise approach if the wing has rigid axis. 
normal ribs, as in this case camber of the streamwisy 
strip will arise. The same is true for the axis-normal 
strip if the wing has rigid streamwise ribs. Modem 
wings, on the other hand, can rarely be considered t 
have rigid ribs, so that this difference between the tw 
approaches is somewhat removed. 

In a few sample calculations the effect of the rib di 
rection on the flutter speed was investigated; in the 
cases considered, there was no effect worth mentioning 
in the practical range of frequency ratios. It was als 
established that this is only true if the comparison is 
consistent, in that camber deformations, where they 
arise, are properly taken into account. In a practical 
case leaving out the camber terms reduced the calecu- 
lated flutter speed by about 15 per cent. 


FLUTTER DETERMINANT OF SWEPT WINGS 


The zero-order effect of sweep on the flutter speed is 
We call first-order effect the 
In calcu- 


described by Eq. (1). 
effect that is found if strip theory is used. 
lating this effect we use the streamwise approach 
Further, we define the unswept reference wing to our 
given swept wing in a manner illustrated by Fig. 1. 

The wing is supposed to have a straight elastic axis, 
the root of which is held rigidly. In the transition from 
swept to unswept wing, each spanwise strip keeps its 
length and is always composed of the same mass par- 
ticles; while it rotates about its root point, this point 
slides along the root chord. The wing ribs are always 
parallel to the air stream. Then, if b, is the half chord 
at a given spanwise station of the swept wing, and ) 
the half chord of the straight reference wing; 


b, = by sec A (2 


Due to the skewness between elastic axis and mb, 
bending of the axis produces both vertical translation 
and pitch of the rib, the pitching motion being propor- 
tional to the local slope of the bending mode. Further- 
more, the moment of the aerodynamic forces acting on 
the rib produces both a torsional moment and a bending 
moment in the elastic axis. The details of the integro- 
differential equation that thus arises are worked out in 
Applying to it the Rayleigh-Ritz method, 
This determinant 


reference 2. 
we arrive at the flutter determinant. 
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we now describe for the case of a uniform swept wing factor depends, of course, on the particular choice of 
with two degrees of freedom, the flutter modes being flutter modes.) 
the fundamental bending and torsional modes of the Now let Eq. (4) be the flutter determinant of the 
uniform beam. unswept reference wing of our given swept wing. The 

As an illustration of the notation we write the flutter flutter determinant of the swept wing then is 
determinant of the two-dimensional problem (wing in ; " 

Aw + (Am + Ai)«k 4 0.96Ao + 1.07AnK 
parallel motion between walls) : 
r 1.16A qk" “ee = 0 
Aw — x Ay = 
, 0) (3 0.96A, + 1.07A 1K Ay -y 
Ay Aj = 2 ss 
(o 
re a4 Ss Se (w y- a: vy , : , g = 
Here A ih Ln, ind so = ' A@B/@)" By J If the terms with «x are disregarded, Eq. (5) looks ex 
or/w)” Mr. » is the non-dimenstona wing mass, 7 6 actly like Eq. (4 The A,, also are the same coef 
rotational — and ws and or ase the ancoupied ficients as the A,, in Eq. (4); this means A Mo 
natural frequencies. " 4 L,, and so on, with yu referred to the halfchord by of 
3 : ‘ ewent We ; ; é ; 

We turn now to the ore and pe wep cman , vel the unswept wing and thus invariant with respect to 
wing. Its flutter determinant is O tained i certain sweep. Further, wz and w7 (in x and y) in Eq. (5) are 
integrals alate bending ae and amano mode are still the natural frequencies of the unswept wing 
introduced - Eq. . ees ue pam by dividing However, the flutter speed is now connected to the 
through the factors of Aw and Ay, we obtain frequency parameter & (of Ls, etc.) by the relation 

Aw — x 0.96 Aa (0 (4 V" wh,/k (why sec A)/R (5a 
0.96 Ay Ay = as : b , 
Thus Eq. (5) contains the zero order approximation, 
The nearness of the factor 0.96 to 1.00 shows the close Eq. (1). 
relationship between the two-dimensional and_ the The difference between zero-order and _ first-order 
cantilever flutter problems. (The exact value of this approximations is determined by the terms with x. « 
E at 40c C.G. at .40c 
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is the swept wing parameter; it combines sweep angle 


A and aspect ratio of the unswept reference wing A Ry: 


K 2 tan A/AR» (5b) 


Terms with « arise from both the effects mentioned 
above—from the pitch of the rib that is coupled to the 
slope of the bending mode and from the bending 
moment in the elastic axis that arises from the rib 
moment. 
that both effects are governed by one and the same pa- 


The fact that x contains the aspect ratio A Ry 


The simplicity of Eq. (5) is due to the fact 


rameter k. 
is readily explained. The rib moment is measured in 
terms of the wing chord, while bending moments are 
measured in terms of the length of the elastic axis. 

The factor 2 in Eq. (6) has been introduced in order 
to let all the numerical factors in Eq. (5) be of the order 
one. Thus the numerical value of x is the order of 
magnitude of the second-order effect. An example: 
for A 15° and AR, 0.4. 


Eq. (5) refers to the aforesaid uniform wing and to a 


5 we find x 
particular choice of modes. The numerical factors de- 
pend upon these modes, and certain spanwise means 
of the derivatives take the place of the A,, if the wing 
uniform. 


is not However, even for the nonuniform 


wing, Eq. (5) remains a fair illustration. 


SAMPLE CHARTS OF RESULTS 


Figs. 3-5 are sample charts of results of Eq. (5); 


this means of a two degree of freedom calculation, 


Chart of reference speeds, c.g. in front of E-axis 


See Fig. 3 for legend 


based on aerodynamic strip method, and referring t 


the uniform swept wing. Over the frequency ati 
wp/wr 1s plotted the critical value of the reciprocal rei 


erence frequency 


| @7 { V cos A) wb V A V neni wrb Vv f ( 


Here .J/ is the Mach Number of the aerodynamic coel 
ficients and accordingly the Mach Number of the axis 
normal stream component. (See the definition of the 
strip method in the section on axis-normal and strean 
wise approaches. ) 

For an elastic axis at 40 per cent chord and for ont 
position of the weight axis each chart shows flutter 
speed and divergence speed for five Mach Numbers 
for three values of the sweep parameter «x, and for thre 
values of the mass coefficient u m/ ap be. 

Figs. 3-5 are samples of a set of 20 charts now 1 
preparation covering five positions of the /-axis with 
four positions of the c.g. axis. Further sets might be 
prepared to cover the effects of taper or of concentrated 
masses. Once completed, such sets will be useful for 
the coordination of experimental results, as guidance 
in the design stage and as indications of ranges where 
the flutter speed is sensitive to small variations of some 
parameter. 

There is not enough space here to discuss the sample 
charts in all details. Generally speaking, the swee{ 
parameter x tends to decrease 1/a@7. This tendency) 
can be traced to the existence of a solution for 1/7 


¢l|- 
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(, this means, to the existence of flutter of a swept wing show this, we write the upper left-hand term of Eq. 
with infinite torsional stiffness. This freak solution is (5) in the form 
discussed in the next section. The tendency, which is ‘ ’ ) 
roughly linear in «, in connection with the factor cos LEG Aa + ee FA) (el + Seale) 50 

\ in Eq. (6) has the effect that the minimum flutter The term inside the braces is almost identical with the 
speed does not usually occur at zero sweep but occurs pitching moment about an axis b/« forward of the quar 
ta moderate positive sweep angle A. This conclusion ter-chord axis. Glauert flutter occurs if both real and 
is well confirmed by experiments." * imaginary parts of this moment are positive. Bending 
As an example of the usefulness of such charts con flutter of the swept wing occurs if the term inside the 
sider the effect of increasing the Mach Number .\/ for bracket fulfills this condition. In both cases, a large 
siven value « and a practical frequency ratio, wg/w7 mass coefficient u is required to make the real part posi 
15. The value of 1/7 will remain roughly constant tive. 


ip to transonic speeds, decreasing somewhat if both x 
nd uw are large but increasing 1f both are small, and in WING-Bopy FLUTTER OF SWEPT WINGS 


a 
reasing at supersonic speeds. ‘a ; : : : 

rhe essence of the preceding considerations was de 
. . veloped during and shortly after the war There 
BENDING FLUTTER : ‘ites 
was an obvious conclusion to be drawn: pure bending 































































































\ solution for 1/@7 0 indicates that flutter occurs flutter, although a freak theoretical occurrence, had 
ven with infinite torsional stiffness. This means it been shown to affect the flutter of swept wings in a wide 
indicates pure bending flutter. In Figs. 3—5 this occurs parameter range; it might affect even more the flutter 
it transonie speed only. However, this limitation is of the combination of a swept wing with a heavy but 
' due to a limited range of uw. With very high values of free body. This possibility was investigated at the 
u, pure bending flutter occurs down to incompressible RAE in England. As suitable calculating facilities 
speeds (theoretically at least). In fact, the pure bend- were not then available, the experimental rig shown in 
J I 
ing flutter of swept wings is the same physical occur- Fig. 6 was designed. It is described in a paper by 
rence as the pure pitching flutter of straight wings hav- Jordan and Smith.6 The arrangement shown here is 
ng a forward pitching axis. (This occurrence—dis- for symmetrical flutter. A floating platform carrying 
covered by Glauert* in 1929, and sometimes called the vertical wing is arranged by means of three vertical 
negative damping’’—is being widely discussed). To legs with universal hinges at both ends. The platform 
4 | E at.40c C.6. at .475¢ 
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wr = Fic. 5. Chart of reference speeds, c.g. behind E-axis. See Fig. 3 for legend 
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Wind-tunnel arrangement for symmetrical wing-body 
flutter 


Fic. 6 


is held in its plane by two pairs of horizontal springs. 
Mechanically, these springs are negative masses. 
amount of these negative masses depends upon the 
flutter frequency which is not known a priori. There- 
fore, it is important that the springs are arranged sym- 
metrically with regard to the center of gravity of the 
oscillating unit thus leaving its position invariant. 
Through variations of spring stiffnesses, distances, addi- 
tional weights, etc., large parameter ranges can be 
covered. 

This test rig has proved very useful. 
with test programs done on it will be 

A copy of the rig has been built in 


A number of 
reports dealing 
released shortly. 
Sweden.’ 

A sample result is shown in Fig. 7. Here ‘“‘center 
position” refers to the case where the overall c.g. coin- 
cides with the aerodynamic neutral point. In the 
‘front position’ the c.g. is about 1/5 mean chord in 
front of the neutral point. 


(1) wing-body flutter may reduce the flutter speed con- 


Two main conclusions arise: 


siderably but occurs only if the body radius of inertia 
is small; and (2) the position of the overall c.g. is an im- 
portant parameter. 

[At the Second International Aeronautical Conference 
in New York in 1949, Broadbent* reported on the con- 
siderations that led to the test program (that was then 
in progress) and on the test program itself. He also 
showed some numerical calculations that he had made. 
In these calculations the body c.g. position was kept in- 
variant while the sweep angle was varied. Hence, the 
results of Broadbent’s numerical calculations have to 
be used with caution. ] 

T-TAIL 


RESULT CONCERNING 


FLUTTER 


AN EXPERIMENTAL 


Fig. S shows that the strip method can yield correct 
results even if the sweep parameter x is not small. 
Shown are results for a T-tail configuration. Fixed to 
the top of a vertical fin is a horizontal stabilizer which 
= 0.3. The 
main aerodynamic forces come from this stabilizer, 
while the fin spar is the elastic member of the configura- 


is swept and has the sweep parameter x 


tion. The fin root is rigidly held. 

The experimental points are taken from a test pro- 
gram that was undertaken in cooperation with the U. S. 
Navy Bureau of Aeronautics. Varied are fin sweep 


The 
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All test 


are referred to an invariant coupled fin torsional fy 


angle and stabilizer dihedral angle. 


quency. 

The main experimental result is that fin sweep is a 
vantageous from a flutter point of view. This conely 
sion remains true even if the results are referred not to 
constant fin frequency but to a constant fin stiffness 
Another experimental result is that stabilizer dihedral 
has an increasing effect at larger sweep angles. 

The theoretical curve that is also shown in Fig 
refers to zero dihedral. It was calculated using th 
streamwise strip method with fin bending and fin tor 
sion as degrees of freedom. There is agreement for 4 
A = 30 
this is readily explained by the observatio, 


O° and There is disagreement for 
A = 60 
that second-order bending took part in the flutter ny 
tion at this sweep angle. (Calculations including sec 


ond-order bending are being undertaken.) 


IMPROVEMENTS OF THE FLUTTER CALCULATION 


Even if the wing can be represented reasonably well 
by means of an elastic axis, the simplified calculations 
used for Figs. 3-5 often need to be improved by the 
introduction of additional modes and the use of three 
dimensional aerodynamic methods. 


On Checking the Modes 


In the present age of high-powered calculating ma 
chines it no longer seems to be a problem to take a fen 
more modes “‘just in case’ and thus to make fairly sure 
that the flutter motion of the actual wing is adequately 
represented. In practice, however, additional modes 
require additional preparations, and neither machines 
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ON TFHE FLT 


nor manpower nor time are always available. Thus 
it may uot be out of place to recall an old but little- 
known trick? which allows checking the adequacy of a 
given set of modes without increasing their number. 

The integro-differential equations of the flutter prob 
lem 
wise point of the wing. 
by reducing this set of (7 X ©) 


Usually this is done by the Rayleigh-Ritz proc 


represent ” conditions of equilibrium for each span 
We form a flutter determinant 
conditions to ” condi- 


thon 
OMS 


ess which, mathematically, means multiplication with 


virtual modes and subsequent partial integration. 
What we do in effect is to take a weighted mean of the 
infinite number of conditions. 

However, there is no basic reason why any other 
weighted mean could not be taken instead. In real 
valued problems, the argument is that the Rayleigh- 
Ritz mean yields an error in the unknown frequency w 
that 1s proportional to the square of the error in the 
modes, while any other mean yields an error in w that is 
linear with respect to the error in the modes. Cum 
grano salis, this position still applies in the complex 
flutter problem. However, considering the degree of 
accuracy that we require from a flutter calculation, the 
linear error might well be the smaller one. (See the 
example below. 

The trick, then, is to form two different means, using 
the same set of modes each time. If both determinants 
vield (nearly) the same values for both flutter speed 
and flutter frequency, then very likely we have a relli- 
able flutter result. If the two sets of results disagree, 
then something is wrong. 

An example, admittedly a freak one, is presented in 
Fig. 9 


[Two two degree 


The wing is unswept but strongly tapered. 


freedom results are shown; E 


and refers to the Rayleigh-Ritz 


of 


stands for “‘energy”’ 


method, D stands for and was obtained 


“deformation” 


by calculating the deformation of the wing tip. The 
Aw + (0.99 Aor + 1.07 Ayo)x + 1.21 
0.97 Ay + 1.04 Ain k 


The smallness of the differences between the two deter- 
minants indicates that two degrees of freedom are usu- 
ally sufficient for uniform wings. Larger differences 
arise with strongly tapered wings and in particular in 


the case of concentrated masses. 


A Three-Dimensional Aerodynamic Method 


An improvement of the streamwise strip approach 
has recently been developed by van de Vooren and 
Eckhaus* who take the spanwise slope of the lift dis- 
tribution into account. 

This, however, means a complication with a rela- 
tively small reward. In view of the modern trend of 
wing plan-form designs it might be more practical to 
develop a proper lifting surface method that works well 
and quickly in particular for small aspect ratios. If 
we know the two ends of the aspect ratio range, the 
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0 
15° »stabilizer dihedral, experimental result 
30 
calculated for O° stabilizer dihedral 
“correct’’ solution was obtained using five degrees of 


freedom, three in bending and two in torsion. 

The strong disagreement between the three solutions 
that Fig. 9 shows is not, in fact, due everywhere to a 
true unsuitability of the modes. The parameters of 
the wing are such that a slight change in the position of 
c.g. axis produces a large change in the flutter speed 
It is useful, though, to have a method that reveals the 
existence of such a sensitivity. At low values of the 
frequency ratio the second bending mode plays an addi 
tional role. At high values a branch involving second 
torsion appears in the correct solution. 

The D-determinant for the uniform swept wing that 
corresponds to the E-determinant Eq. (5) follows for 


comparison : 


0.97 Au + 


A «eile 


1.15 Ayn K 


medium range can, if necessary, be covered by inter- 
polation. 

The differential equation connecting downwash func 
tion w and pressure function p on the lifting surface 


can be written 


1 ( (pen 

W (Xx, ¥ x 
4rd Js ieee 
G [(é —x)k, |\y — n k]) dé dn (8) 
In incompressible unsteady flow 
- ins F eds 
G (u,v) =e (Sa 
J u/?(S* + ] 


The fairly formidable integral Eq. (Sa) becomes much 


simpler in the steady flow case, withw =v = 0. Then 
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G (0, 0) 


I(x — &)/V (x £)? ++ (y — n)?] (8b) 
A very successful method for this case has been de 
veloped by Multhopp." This method offers a num- 
ber of advantages. It is a continuous method, which 
means that all distributions, including the pressure dis- 
tribution, are represented by continuous functions; 
thus, no singularities arise at the borders between two 
lift boxes or similar units. It makes maximum use of 
the result of the two-dimensional theory and of the 
number of pivot points that it uses. The application is 
fast and convenient once the necessary tabulations of 
influence functions are available. 

The influence functions in question are suitable chord- 
wise integrals of the G-function. In Multhopp’s 
steady flow method, they are functions of two continu 
ous parameters, spanwise distance and chordwise dis 
tance. At first sight, extension of this method to un 
steady flow did not seem feasible; it seemed necessary 
to tabulate the influence functions for a third continu 
ous parameter as well—namely, the local frequency 
parameter k. It has been possible however to circum- 
vent this difficulty by choosing a fixed number of pivot 
points. These are shown in Fig. 10. The two pivot 
points that appear at each spanwise station play the 
same role with respect to lift plus moment as does the 
three-quarter chord point with respect to lift alone. 
However, while the three-quarter chord point loses its 
advantage soon if the frequency parameter k differs 
from zero, the present two points keep their advantage 
up to much higher values of k. 

Before applying the method, the given downwash is 
split into a downwash of standard form and a remainder 
which, in two-dimensional unsteady flow, would pro- 
duce neither lift nor moment. This remainder can be 
shown to have little three-dimensional inductive effect 
and is treated by strip method. The standardized 
downwash forms the loads in a system of eight linear 
complex equations. The coefficients of this system of 
equations are read from graphs. 

The method is simple and fast and can be applied 
without difficulty and in a routine fashion to any plan 
form and to any mode of deformation. It yields the 
complete pressure distribution and is thus suitable for 
The graphs required for its appli- 
0 and a number of 


flutter calculations. 
cation have been prepared for \/ 
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Pivot points in three-dimensional method ) 


values of the frequency parameter formed with thy 
wing span. Extension to compressible flow is planned 
Results for rectangular wings have been compare 
with the results obtained by Lawrence and Gerber. 
Good agreement was obtained up to AR = 4. Exten- } 
sion of the method to higher aspect ratios is feasible 
for instance, by increasing the number of spanwist 


pivot stations. 
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7 Lower Buckling Stress of Circular Cylindrical 
Shells Subjected to Torsion 


YOSHIMARU YOSHIMURA* ano JUN’ETSU NITSAWA** 


et 
Tokyo University and Nthon Unwersity, Respectively 
aoe 
} SUMMARY buckled surface far apart from the initial one. This 
Phe critical stress for torsional buckling of circular cylindrical is the reason plane plates do not reveal Durchschlag, 
ined from the linear theory with small deformation is but circular cylindrical shells do. The fact that the 
iderably greater than the experimental buckling stress buckled surface is applicabie to the circular cylindrical 
present paper clarifies the cause of this discrepancy between surface means that it is developable, because the latter 
I experiments on the assumption that the torsional . * oy 
ae oy ; is developable. In fact, the buckled surface of circular 
uckling well as the compressive buckling, is a Durchschlag ‘ : ‘ ; ; 
nomenon 3ased on the experimental fact that the buckled cylindrical shells under axial compression or torsion 1s 
face is approximately developable (the mean surface being approximately developable. In particular, in the case 
ped), the equilibrium state after buckling was calcu of torsional buckling, the developable surface is realized 
[by means of the slain premegee oF potential energy, because the mean surface of the buckled waves is not 
. for the case of buckling under constant load, the lower . s 
ae oye cylindrical, but somewhat drum-shaped, and resembles 
uckling stress was obtained, which is considerably smaller than . ‘ ; ? ‘ 
tl upper buckling stress and which explains the experimental a hyperboloid of one sheet of revolution. 
sults approximately In the case of buckling under constant According to the above reasoning, the torsional 
gle of torsion there exists no lower stress buckling of circular cylindrical shells was treated as a 
with th ‘ ; 
dee Durchschlag phenomenon by means of the energy 


» planned INTRODUCTION 
compared 


Gerber if IS WELL KNOWN that the theoretical buckling load 
ixtes. of circular cylindrical shells subjected to torsion, 


ee based on the small deformation theory, does not coin 

- feasible, - 

spanwise cide with the experimental results. For example, ac- 
cording to the Donnell’s calculation,' which is repre- 


sentative, the experimental buckling load is about 60 
to 70 per cent of the theoretical value. 
The authors believe that this discrepancy between 


fret theory and experiments in the torsional buckling can 
be ascribed to the same mechanism—.e., the Durch- 
raft Flut schlag phenomenon —as in the case of compressive 
ons 1012 | buckling, which was clarified by one of the authors’ 
Chapts alter references 3 and 4.  Durchschlag’s is the phe- 
i as RAI nomenon of passing from the stable state just before 
Wi buckling to that stable state after buckling, which is 
far apart from the former, jumping over an unstable 
ng Aer state between them. And whether the stable state 
after buckling with minimum potential energy exists 
7 ae depends upon whether an intensely buckled surface 
— approximately applicablef to the initial undeformed 
s. British surface exists. Further, the existence of this applica- 
ble surface depends on both the geometrical shape of 
rhalter the initial surface and the type of load applied. While 
ho in the case of buckling of plane plates there is no 
Problem | buckled surface applicable to the initial one, in the case 
rntaional of the circular cylinder there exists an applicable 
a Received November 15, 1955 
eal “Pr fessor of Applied Mechanics, Institute of Science and 
lechnology, Tokyo University 


* Lecturer, Faculty of Technology, Nihon University 


i Iwo surfaces that have the same first fundamental quad 
aiid tic form but have different second fundamental quadratic forms 
canal said to be applicable or isometric. See, for example, Eisen 
11 ‘a rt, L. P., Introduction to Differential Geometry with Use of the 

- Te ( lus, p. 147, p. 215, Princeton, 1947 
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method based on the finite deformation theory and 
results satisfactory for explaining the experimental re 


sults were obtained. 


Basic RELATIONS 


The fundamental equations for the finite deforma 
tion, which determine simultaneously the stress func 
tion x and the displacement w normal to the undeformed 
cylindrical surface, are given by the compatibility 
condition 


[(0?/Ox?) + (07/dy?) }*x E[(0°w/Oxdv 
(0?w/Ox") (O°w/dOy") — (1/a 


(O-w On , 
and the equilibrium equation 
[Et?/12(1 — v?)] [(0?/Ox?) + (07/Oy?) ]*2 
p + (0°x/Ov") (O?w/Ox*) — 2(0°x/OxXOY) X 
(O° /Oxdy) + (02x /Ox?) [(1/a) + (d°w/dy 2 


in the direction normal to the deformed middle surface, 
where x and y are measured in the axial and circumfer- 
ential direction in the middle surface of the cylindrical 
shell, a denotes the radius of the shell, ¢ the thickness, 
FE Young's modulus, v Poisson's ratio, and p the ex 
ternal force normal to the deformed middle surface 

It is obviously difficult to solve the simultaneous 
equations for w and x, so the minimum principle of the 
potential energy is applied as usual instead of Eq. (2 
The total potential energy with respect to the half 
length / of the cylinder is the sum of the extensional 


elastic energy 


Wy (f 
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the bending energy 


W it® ‘f f~ (= 4 4) 
, 24(1 — pv”) J0 J70 | Ox? oy’ r 
O*w O*w O-w \2 
v) | a ee ( ) | dx dy (4) 
Ox? Oy? Oxdy f s 


the 


2(1 + 


and the potential of external force (torsional 


moment) 


) 


W3= -t | (Try) x dy | y dx (5) 
« J0 


0 


where the stress components o,, 7;,, ¢, are obtained 
from the stress function x by the relations 
= Q*y/Ox? (6) 


a, = O*x/Oy", Try = —(0*x/Oxdy), a, 


and y denotes the angle of torsion and is given by 


y = 


y = (Ou/Oy) + (Ov/Ox) (7) 
where and v denote the axial and circumferential dis- 
placements, respectively. Therefore, from 
try = [E/2(1 + v)] [(Ou/dy) + 

(Ov/Ox) + (Ow/Ox) (Ow/Oy) | 


we have 


y = [(2(1 + »)/E]rz, — (Ow/Ox) (Ow/Oy) (8) 


THE ASSUMPTION FOR DEFORMATION 


To obtain a solution as exact as possible by the energy 
method, it is required to assume the deformation as 
close as possible to the actual one. According to the 
experimental results,’ the surface buckled under torsion 
is regarded to be expressible by the superposition of the 
helical wave pattern upon a drum-shaped surface, but 
That the 
average surface of the buckled waves is drum-shaped 
or one sheet hyperboloid-shaped is due to the fact that 
the wave amplitude at the middle part of the cylinder is 
Although the decrease 


not upon the initial cylindrical surface. 


larger than that near the ends. 
of the wave amplitude toward the ends seems to be due 
to the end conditions, that is not the case; the decrease 
is regarded as resulting from the fact that the buckled 
surface is strengthened to become developable or in- 
extensional. 

From the above, the displacement w can be assumed 


as 


w/a = fo[l + cos (x/l)x] + 
fi cos [(Ai1/a)x + (n/a)y] + 
fo cos [(A2/a)x + (n/a)y] (9) 


where afo, afi, and afz denote the amplitudes, 27a/\, and 
2ma/d2 the wave lengths in the axial direction, 2/ the 
length of the cylinder, » the circumferential wave num- 
ber, and x is measured from the middle of the cylinder 


in the axial direction. The first term in Eg. (9) rep- 
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resents the drum-shaped mean deformation, and th, 


terms following represent the wavy deformation super. | 


posed upon it. 

Eq. (9) satisfies neither the conditions of support, 
ends nor of fixed ends perfectly. However, since it has 
already been shown by L. H. Donnell! that the buckling 
stress does not depend much on the end conditions 
simply supported or fixed, we need not consider th 
difference between them. 
that Eq. (9) is not required to satisfy any condition 
other than 


w= 0 x + l If 


In order that the condition (10) may be fulfilled by | 


Eq. (9), 
f, cos [+(A,/a)/ + (n/a)y] + 
ff cos [+ (Az a)l + (n a)y] 0 (i 
or fi cos (A;/a)l + f2 cos (A2/a)l = 0) 
fy, sin (A; /a)/ + fo sin (As/a)l = OF 
Since f; and fe do not vanish for the buckled surface 
Eq. (11’) reduces to 


cos (A;/a)/ cos (A:/a)/ 


.; = 0 (12 
sin (A;/a)/ sin (As/a)/ 
or sin (Ay — Ae) (//a) = 0 (12’ 
so that (Ay — As)(i/a) = 1x7, += 1,2,3,... (8 
From Eqs. (11’) and (13), 
fo cos (A;/a)/ 
fi cos (Av/a)/ 
sin (A;/a)/ {—1 i: even 
sin (Ao/a)/ (+1 i: odd | 


Because the cylindrical shells buckle with the mini- 
mum value of A; — A2,®7 = 1 from Eq. (13), and, there- 
fore, from Eq. (14), 

h=h=f (15 


indicating f; and f/f. by f. By means of Eq. (15), Eq. 


(9) can be written as 


a 


w = afol[l + cos (x//)x] + af; cos[(A;/a)x + 
(n/a)y| + cos [(As/a)x + (n/a)y]{ (16 


each term of which satisfies the condition (10) indi 
vidually. 

Now it should be remarked that the first term in Eq 
(16) satisfies the condition of fixed ends Ow/dx = 0 
for x = +/, but not that of simply supported ends 
0*w/Ox? = Oforx = +/. The second and third terms 
satisfy neither of the conditions exactly, but are con 
sidered to satisfy either of them in the mean over the 
circumference. Therefore, (16) is regarded as 
representing the deformation approximating the case 
of fixed ends rather than that of simply supported ends. 


Eq. 


TOTAL POTENTIAL ENERGY 


Substituting Eq. (16) into Eq. (1), we have 


Consequently, we assuny | 
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92/dx2) + (07/Oy?) |? (x/£) = fo(r/l)? cos (x/l)x + if(Ay/a)? — (1/2)fof(n/a)* [x(a/l) 7; xX 
cos [(A,/a)x + (n/a)y] + }f(r2/a)? — (1/2)fof(n/a)? [x(a/l) 2} cos [(A2/a)x + (n/a)y 
(1/2) [fn(m/l) }?( cos | [Cr + Az) /a lx + 2(n/a)y} + cos (x//)x) 
(1/2) fof [n(a/1) ]? (cos } [(2A1 — Az) /alx + (n/a)y} + cos }[(—Ar + 2d2)/a}x + (n/a)y; 17) 


The particular s¢ ution of Eq. (17) is given by 


y/E A/(m/1)?| cos (r/l)x + Ba? cos [(A;/a)x + (n/a)y| + Ca* cos [(Ao/a)x + (n/a)y| — 
Fa? cos }[(A1 + A2)/a)]x + 2(n/a)y} — Ga® cos } [(2A; — A2)/a]x + (n/a)y; 
Ha? cos } [(—A1 + 2d2)/alx + (n/a)y) + (7/E)xy (18) 
where A = fy — (1/2) (fn)? 
B = [fd? — (1/2) (fofn?) (wa/l)?|/(x? + n?)? 
C = [fr2® — (1/2) (fofn?) (wa/1)?]/(A2? + n?)? ™ 


F = (1 2) ({n)° (3a f)?} [(Ay 4. he)? 4. in? ]? 
G = [(1/2) (fofm®) (wa/1)?]/[(2M — Az)? + n?}? 
H = [(1/2) (fof?) (wa/l)?]/[(—dr + 2d2)? + n°}? 


and 7 represents mean shear stress. The stress components in the middle surface are obtained, from Eqs. (6 and 


1S), as 
E — Bn? cos |(A;/a)x + (n/a)y] — Cn? cos [(A2/a)x + (n/a)y| + 
tFn? cos } [(A1 + A2)/a] x + 2(n/a)y} + Gn? cos ; [(2A1 — Ae) /a]x + (n/a)y;y 4 
Hn? cos } [(—Ai1 + 2A2)/alx + (n/a)y} 
1,/E —A cos (m/l)x — Bd,? cos [(Ax/a)x + (n/a)y] — Crs? cos [(A2/a)x + (n/a)y] 4 
F(A, + Az)? cos } [(Ar + Av) /a]x + 2(n/a)y} + G(2A, — Az)? cos } [(2A1 — Az) /alx + (n/a)y; 4 (20) 
H(—); + 2d2)? cos } [(—A, + 2Az)/alx + (n/a)y; 
tr, /E — Bnd, cos [(A;/a)x + (n/a)y] — Cndrz cos [(A2/a)x + (n/a)y|] + 
2Fn(A; + d2) cos } [(Ar + Az)/a]x + 2(n/a)y} + Gn(2d1 — dz) cos } [(2A1 — Az) /alx + (n/a)yy + 
Hn(—)d, + 2d2) cos } [(—Ay + 2Ar2)/alx + (n/a)y{ — (7/FE 
While 7,, is composed of the uniform stress 7 and the The potential energy can now be calculated by means 
periodic terms referred to x and y, o, and a, have only of the results for the displacement w and the stress 
the periodic terms. Since the cylinder ends are usually components o;, T,,, and g,. Introducing the non- 
reinforced by rigid rings in experiments, the displace- dimensional quantities 


ments u and v at both ends should have constant value 


: iat ae n = ¢?/(12(1 — v’)a?] 
independent of y. These boundary conditions regard- 


ne 
& = foVn g =J/V1 


ing w and v will be satisfied by the general solution of ‘ ‘ 
: me ¥ : a ae . ‘ 3 = Aiv U] e2 Ae U] 9 
Eq. (17), which is the sum of the particular solution (20) ‘ (21) 
: : - = ky ~~ ko = mra/l)n n 
and the complementary solution of Eq. (17). But, _ 
- d = n 


since the effects of these boundary conditions are neg- 
ligible and it is not the purpose of the present paper 
to discuss them, Eq. (20) is sufficient for the stress com- the extensional elastic energy Eq. (3) divided by 
(1/2) zaltEn can be obtained, by means of Eq. (20), as 


@ = (7/E)/Vn 
ponents in the following calculation. 


I, = W,/[(1/2) raltEn] 
= gy? + gl [Ry4/(Ri? + N22] + [Rot/(Ro? + N2)2]} — gog?N2(1 + wu?) [Ri?/(Ri? + N?)?] + (99) 
ho? /(Ro? + N?)?]}) + (1/4) g4*N4(1 + fut/[(Ri + Re)? + 4N2)2}) + (1/4) g02g?N4u4((1/(Ri? + N*)?] + < 
(1/(ko? + N2)2] + {1/[(2k, — ke)? + N2]2 + $1/[(—k + 2k)? + N2]°t) + 4(1 + ve? 





and the bending energy, by means of Eqs. (4) and (16), ; —— bhai :; , 
sae y dx in Eq. (5) becomes a periodic function of 
0 


as 
. ; y. This means that the circumferential displacement 
TT. = We/[(1/2) xaltEn] oo ’ : a ‘ 
= p2[(k,? + N2)? + (ke? + N2)2] + go2u! (23) at the cylinder ends is not constant but depends on y. 
all it ; : cai _ This result is due to the fact that the boundary condi- 


In calculating the potential of the external force tions for u and v are not fulfilled, as was mentioned 


Eq. (5)], it is obvious that the periodic terms in 1,, above. Consequently, it is reasonable to take the mean 
4 ° ° 2 a ‘ 
need not be taken into consideration, because they = _ . . 
=a , value (1/27a) dy y dx of the integral 
: 0 0 
have no effect on the integral f (Fislcmn 4g OP. 1 
- _ : ee ‘ [ y dx as the angular displacement of the end 
But if the periodic terms in y [Eq. (8)] are considered, Jo 
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This mean value is shown to be equal to the integral 


| 
| y dx for the case when only the nonperiodic terms 
J 0 


in y are considered. 
to Eq. (5), should be rewritten as 


7 nonperiodic terms of }[2(1 + v)/F]r,, — 


(Ow/Ox) (Ow/Oy); (8’ 


Hence, from Eqs. (16) and (20), 


7 -{201 + v)/E|r — (1/2)f?n(A, + r2) = (24) 


or y’ y/V 1 —2(1 + vd — 
(1/2)g2?N(ki + he) (24’ 
The potential of the external force is obtained as 


Il. W; {1 2) raltkn]| 


—8(1 + v)¢*? — 2¢°N(k; t ko) (25) 


by means of Eqs. (5), (20), and (24). 


STRESS-STRAIN RELATIONSHIP AFTER BUCKLING 


The state that may be actually realized after buck- 


ling must be determined by minimizing the energy, 
with respect not only to the magnitude of deflection 
fo, f but also to the parametric quantities such as \,, 


At first 


it is necessary to determine the stress-strain relation 


A», m, Which characterize the wave pattern. 


after buckling by minimizing the potential energy 
with respect to go and g, taking 1, ko, and N as param- 


eters, before obtaining the equilibrium state after 


buckling. 
From the minimum condition of the total potential 
energy with respect to go, 


(0/Ogo) (Ih + Ile + II3) = O (26) 
K, [kit/(Ri2 + N2)2] + [ot/(Re? - 
Ky N4y?([Ri?/(Ri? + N*)?] + [Re?/(Re? 






























Ky = (1/16) Nu4([1/(k:2 + N2)2] + [1/(ee? + N2)2] + 
/ 2N(ki + k 
1.0 eet 
Ei = SSSii) SS SS see: 
@ FH te | a 
+ ++ + + ++ + ++ + 
asf ttt | ir Ht 
} }} 
ie; e | 
| | | | 2004 (019) g . | | 
0.0 eee Sree oe | 
TT t19.008 (0.10)4—_# TTT OE RE ST | 
0.051 a~a $=03p'(46+ | 7.8+ 24) fixed ends, 1 
Donnell 


b~b g =0.3'(28+ 26+ - ) simply Supported 
ends. 





* Experiments by Donnell (steel ) 
| --°-- Lower buckling stress. 
| --— = Eg. (23) 
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Fic. 1 
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results 
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Therefore Eq. (8), to be applied 


MARGH, 1957 


we have 


{ | | (Ry? \ 
| | 
(ko? + N [(2k; —k \ 
. «a 
(—k 2k: \ os 
For thin cylindrical shells, of which t/a < 1 and a 


uw < 1 by Eq. (21). So, neglecting the small terms 
higher than the fourth order with respect to wu, Eq. (27 


reduces to 


go = (1/2)g?N? (1 + we [Rr?/(Ri? + N?*)?] 
[Ro?/(Ro® + N*)?]{) (27 


As the coefficient A in Eq. (19) is shown to have a 
positive value according to Eq. (27’), the average valu 
—A cos (/l)x of the stress component o, in Eq. (20 
is compressive in the central part —(//2) < x <(//2 
of the cylinder and tensile in the region —/ < x < (//2 
and (//2) <x </near the ends. This effect is due to 
the interaction of the drum-shaped deformation of the 
average surface and the periodic waves. 

The minimum condition of the total potential energy 


with respect to g 
(O/Og) (Il; + Hl, + IIs) 0 (28 
gives the stress-strain relation 
od (1/J) (Ko — Aog? + Ayg'*) (29 


where 


+ N22] + (ky? + NY)? + (bo? + NY)? 
+ N)2] + (1/2)m) [Ri?/(Ri? + N?)?] + [Ro?/(Ro® + N?)?])? | 
) 
(30 


(1/2) }u2/[(ki + Re)? + 4N2]2) — (1/2)? 


11/[(—ky + 2k)? + N2}}) 


(1/[(2k: — be)? + N2]2} + | 


The parameters in Eq. (29) are ki, ko, and NV, but by 
mM k, — k» [Eq. (21)| the stress-strain relation (29) is 
determined with k, and N as parameters, when yu, which 
gives the length of the cylinder, is prescribed. 

With g = Oin Eq. (29), the critical stress for buckling 
or the upper buckling stress in the small deformation 


theory 


re) Ko/J (31 


is obtained. In Fig. 1, the buckling stress given by 
Eq. (31) is compared with both the theoretical result 
and the experiments for steel cylinders by L. H. 
Donnell.! It is found that the result obtained by the 
authors is about 4 per cent smaller than that of Donnell 
for the fixed ends and about 6 per cent larger than that 
for the simply supported ends. The theoretical critical 
stress by both the present authors and L. H. Donnell 


coincides with the experimental buckling stress for 
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Stress @ r/(Ey») after buckling against wave 


unplitude g = v/s 


small u, but does not for large uw, and the experimental 
results scatter in the range of 60 to 70 per cent of the 
theoretical value. 

In order to explain the discrepancy between the 
linear theory and experiments, we must consider the 
For example, in the case of 
0.05), for which the dis 


post-buckling state. 
U 0.1 (a 
crepancy between theory and experiments is large as 


0.7 for \ n 


shown in Fig. 1, the relation between @ and g after 
buckling is calculated by Eq. (29) and is shown in Fig. 
2, with kp and NV as parameters. Practically, the quan- 
tity that is observed in experiments is not the wave 
amplitude g, but the angle of twist y. Thus, the rela- 


tion between ¢ and y’ is shown in Fig. 3. According 
to Figs. 2 and 3, the drop of stress after buckling is not 
so great as in the case of compressive buckling. 
AFTER BUCKLING 


THE STABLE EQUILIBRIUM STATE 


AND THE LOWER BUCKLING LOAD 


The equilibrium state realized after buckling can be 
determined by considering the condition of minimum 
potential energy with respect to the variables k;, ko, 
and .V, or to ky and N with prescribed yu, based on the 
stress-strain relation obtained in the last section. Fur- 
ther, the equilibrium states after buckling decided by 
the minimum condition of potential energy are differ- 
ent according to the conditions of loading—namely, 
the rigidity of testing machines. One of the extreme 
cases of loading condition is the case of buckling under 
constant stress; the other is that of buckling under 


constant angle of torsion. 


(a) Buckling Under Constant Stress 


The equilibrium state after buckling in this case is 
obtained by minimizing the total potential energy 
I] Il; + Il, + Tlz. In Fig. 4, the total potential 
energy II is calculated against the independent vari- 
able ¢, using Eqs. (22), (23), and (25) and the result 
of Fig. 2. The parabola O—O represents the total 
potential energy corresponding to the linear stress- 
Strain relation before buckling, and the curves inter- 
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TABLE | 


Comparison Between Upper and Lower Buckling Stresses 


Stress and 


Lower Buckling 


Upper Buckling Stress and 
Corresponding Values 


Corresponding Values 


of ky and of ky and 
faa) 0.453 0.375 
k 0.02 (1). O4 
\ 0.3 0.3 


secting the parabola O—O represent those corresponding 


to the stress-strain relation after buckling (see Fig. 3). 

Buckling is regarded to occur at such stress as the 
potential energy just before and after buckling is equal, 
jumping over the energy barrier between the two states. 
According to Fig. 4, the stress is given by each point 
which the stress-energy curve after buckling intersects 
with that before buckling O-—O and depends on the 
5, the stress corresponding to 


pl rtted 


parameter ky. In Fig. 
the intersecting point 


Buckling occurs at the minimum stress 


mentioned above is 


against ko. 
among the stresses represented by the curve in Fig. 5. 
This minimum stress, which is called the lower buckling 
0.375 for k: 0.04, with the 
prescribed values of the parameters given in the Figure. 
This result is compared in Table 1 with the upper 
buckling stress obtained from Eq. (31) for yu ek 
The lower buckling stress is about SO per cent of the 
upper buckling stress, and the qualitative explanation 


stress, is obtained as @ 


of the experimental results is considered to be sufficient. 

The drop of stress after buckling, as shown in Fig. 2 
for long cylinders (small uw), is very small, because the 
coefficient A, in Eq. (29) has the factor u Accord 
ingly, the lower buckling stress approximately coincides 
with the upper buckling stress due to Eq. (31) for small 
uw. Considering this fact, the curve representing the 


lower buckling stress is given in Fig. 1. 


(b) Buckling Under Constant Angle of Torsion 


In the case of prescribed angle of torsion, the equi 
librium state is determined by the minimum condition 
of the elastic energy I]; + IIy of the cylinders, not in- 


volving the potential of the external force. The elastic 
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energy given by the area under the stress-strain curves 
in Fig. 3 has different values according to the curves. 
Although the equilibrium state after buckling is given 
by the minimum value of this energy, the lower buck- 
ling stress does not exist in this case because the angle 
of torsion-energy curves after buckling does not inter- 
sect with that before buckling. Therefore, if a per- 
fectly rigid testing machine were used and the shell had 
a perfect cylindrical surface, the buckling would occur 
at the upper buckling stress. 


CONCLUSION 


The essential feature of the buckling of circular 
cylindrical shells under torsion is that it is a Durch- 
schlag phenomenon, as in the case of compressive 
buckling. And this phenomenon (Durchschlag) is due 
to the fact that the buckled surface is approximately 
developable and, therefore, that the equilibrium state 
after buckling, which has as low potential energy as 
that just before buckling, does exist. 

The equilibrium state after buckling was determined 
by minimizing the potential energy with respect to 
all the parameters which characterize the wave pattern 
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of buckling, based on the finite deformation theory 
The results of calculation show that the lower buckling 
load exists for the case of buckling under dead load, 
which is smaller than the critical stress by the linear 
theory and which explains experimental results approxi 
mately; but it does not exist in the case of buckling 
under constant angle of torsion —1.e., by perfectly rigid 
testing machines. 
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Experimental and Analytical Investigation of 


Secondary Flows i 


HOWARD Z. HERZIG* 


AND 


Lewts Flight Propulsion 


SYMBOLS 


mstants, coefficients of polynomials in powers of 


\ 


constant 
F. I Blasius function (Table 1 
1H function of y, J/i(n) = Piln F'(n 
constant, 7 = 0, 1, 2, mI 
function of n, J;(n) = Pdn)/F'(n 
kK constant of integration for stream equation in 
v, y-plane 
constant 
constant 
) = component of boundary-layer velocity normal 
to main-stream flow direction 
se 
function of n, w = > WP 
static pressure 
] Vl? 4+ W? 
component of boundary-layer velocity tangential 
to main-stream flow direction 
| main-stream velocity components in x- and 
directions, respectively 
inlet velocity 
boundary-layer velocity components in the 
y-, and directions, respectivel\ 
A, & components of external force field in x- and 
directions, respectively 
rectangular coordinates 
constant 
similarity variable, » = yV U°®/vx 


viscosity 


coeflicient of kinematic 
density 
Subscripts 
variable index, 0: 4.2; ) 
mstunt 
Superscripts 
variable exponent, = (), 1, 2, j 


constant 

differentiation with respect to 
differentiation with respect to 7 
denotes division by ( 


INTRODUCTION 


' ‘iE IMPORTANCE of secondary flows (here considered 
boundary-layer flows having 
to 


aeronautical research is attested by the great number 


three-dimensional 
components normal to the main-stream direction) 


of recent reports on the results of secondary-flow in- 
In aerodynamic configurations at oper- 
jet en- 


vestigations. 
ational conditions typical for turbomachines, 


Flow in Jet Engines Session, 


IAS, New York, 


the Internal 


Annual Meeting, 


Presented at 
lwenty-Fourth 


26, 1956 


January 23 


* Research Scientist 


in Ducts 


ARTHUR G. HANSEN* 
Laboratory, NACA 


gines, etc., the investigations should be concerned with 


three-dimensional turbulent boundary-layer behavior 


The mathematical complexity of the nonlinear partial 


differential equations which represent such boundary 


layer flows make direct analytical solutions almost 
impossible and certainly quite impractical at this time. 
Accordingly, the secondary-flow investigations have 
been either experimental in nature or have been 
analytical studies of laminar boundary layers 
Experimentally, a great deal of secondary-flow 
information has been provided by the use of flow 


use of smoke, dves, 


naking 


or lampblack traces, as 


visualization procedures, 


paint, well as pressure and 


flow-angle measurements. (Reference 3 in par 


good survey and bibliography of 
By these n 


main-stream 


ticular contains a 


secondary-flow research. eans, the second 


ary flows, which result from turning, 


were found to describe complex flow patterns. For 
typical secondary-flow patterns for annular 
Here the effects 
of the cross-channel (more than main-stream turning ol 
the boundary-layer on the 
and radially inward flows are shown to result in sizable 


fluid in the 


example, 
nozzle cascades are shown in Fig. 1. 
outer shrouds 


inner and 


local accumulations low-energy nozzle 
passages. 

Equally as important as the physical insight into 
secondary-flow behavior provided by these studies is 
the observation that laminar boundary-layer behavior 
can usefully provide qualitative information concern 
ing turbulent boundary-layer secondary-flow behavior 
This experimental observation has far-reaching theo 
retical consequences. The simplifications introduced 
by the 
flow to be laminar make it possible to obtain analytical 


(It should be noted that 


considering three-dimensional boundary-layer 


solutions to the equations. 
simplification of the mathematical problem can 


also 
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Fic. 2. Cross-channel flow in turning 


two-dimensional 90 
circular-are duct. 


be achieved by assuming nonviscous flow. However, 
the existence of a limiting cross-channel secondary-flow 
streamline in the direction of the resultant shear stress 
at the wall, observed in references 1 and 2, eliminates 
Refer- 


ences 6-9, for example, present theoretical analyses of 


consideration of the nonviscous-flow theories. ) 


the laminar boundary-layer cross-flow developments 
over a nearly flat surface having a leading edge. In 
references 7 and 8, exact solutions are obtained for cases 
of main streamlines consisting of parabolic-shaped 
translates (i.e., the entire streamline pattern can be 
obtained by translating a parabolic-shaped streamline 
in a direction parallel to the leading edge) with no re- 
and 9, 


the turning. In references 6 


obtained for 


strictions on 


approximate solutions are somewhat 
more general main streamline shapes, restricted to re- 
gions of small turning. 

The present report is directed toward obtaining a 
deeper insight into the secondary-flow phenomena in 
curved ducts, as a result of experimental and analytical 
studies. Experimentally, the laminar boundary-layer 
cross-channel flows and vortex roll-up,' together with 
their relations to the main-stream flows and the suction- 
and pressure-surface boundary layers, were investi- 
smoke-flow-visualization 


gated herein by means of 


procedures in several such ducts. Analytically, assum- 
ing laminar main streamline translates, this report 
extends the analyses of references 6-9 to obtain exact 
solutions for general streamline paths without restric- 
tions to regions of small turning. Tables of special 
functions of the Blasius similarity parameter required 
for computing such boundary-layer cross flows are 
presented. The streamline paths at various positions 
in the boundary layer are described for thick and thin 
boundary layers. Several interesting examples of cross- 
channel flows are computed. Certain 
flow effects are then observed experimentally and the 
results compared with the theoretical predictions 


cross-channel 


the only way to provide assurance of the physical 
significance of the theoretical solutions for laminar 
flows. 


AERONAUTICAL 


SCIENCES MARCH, 


EXPERIMENTAL SECONDARY FLOWS IN Ducts 


Flow-visualization procedures were employed in refer 
ences | and 2 to determine the secondary-flow patterns 
in various nozzle configurations. It was noted in these 
references that turning of the main stream in a channel 
or blade passage always results in overturning (more 
than main-stream turning) in the end-wall boundary 
This flow ' 
upon the suction surface typically gives rise to a vortex 


layer. cross-channel where it impinges 
roll-up, the ‘‘passage vortex”’ so-called, near or on the 
suction surface. In the present investigation, smoke 
was used to trace and compare the cross-channel flows 
in two-dimensional channels, the relations between the 
occurrence of their vortex roll-ups and the boundary- 
layer development on the suction surfaces, and _ the 
effects of these vortex roll-up boundary-layer accu 
mulations on the main stream. In addition, the 
smoke-flow-tracing method is used to depict the com 
plicated boundary-layer flow behavior in three-dimen 
sional ducts which accelerate and turn the flow inde 


pendently. 
Cross-Channel Flow and Vortex Roll-Up in 


Two-Dimensional Ducts 


Vortex Roll-Up—lIt 
that the formation 


was suggested in reference 2? 


and location of a passage-vortex 
roll-up in a curved passage depends upon the boundary- 
layer growth on the suction surface. This in turn isa 
function of the main-stream static-pressure distribution 
and velocity profile along the suction surface. For 
purposes of comparison, then, the secondary-flow 
patterns were investigated in a two-dimensional (rec 
tangular cross-section) 90°-turning circular-are duct 
(Fig. 2) and in a two-dimensional 90°-turning acceler- 


» 


ating duct (Fig. 3). The latter is a Lucite model of a 


rectangular elbow designed": '! for a prescribed ve- 


locity distribution that decelerates nowhere along the 


















pressure or suction surfaces in an effort to avoid 

boundary-layer separation. 
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The secondary-flow patterns in the circular-are duct 


) the that 


were quite commonplace (Fig. 2) in 
there were the boundary-layer overturning on the end 


sense 


walls and the passage vortex roll-up in the suction 
surface end-wall corner so frequently observed pre- 
viously! The 


secondary-flow patterns in the accelerating rectangular 


that it requires no repetition here. 


90° turning duct disclose some interesting differences. 
As can be seen in Fig. 3, a directly overhead picture, 
the end-wall boundary layer overturns and forms the 


In Fig. }, look- 


ing in through the top and pressure surfaces, one can 


passage vortex near the suction surface. 


see that the cross-channel flow does not roll up near 
(The lines forming the squares on the 
suction surface are reflections on the the 
stream and potential lines drawn on the end wall. 


the end wall. 
Lucite of 


Instead, it flows spanwise out onto the suction surface. 
The actual vortex roll-up is seen in Fig. 5, looking up 
stream toward the exit of the duct, to occur at about 
one-third blade-span height. 

The suction-surface boundary layer was seen to flow 
smoothly along, remaining quite thin far downstream 
from the inlet. At about 20 1), the 
suction-surface boundary layer begins to be deflected 
toward mid-span. When 
been reached, the spanwise deflection levels off at 


of turning (Fig. 


about 70° of turning has 


nearly one third the passage height. Because the same 
kind of spanwise deflection toward mid-span is occurring 
at the top of the duct (near the Lucite end wall through 
which these pictures were taken —seen in the pictures 
to be held in place by bolts), the entire suction-surface 
boundary layer at the inlet is crowded into the middle 
third of the surface by the action of the cross-channel 
flow. Even when thus thickened, few signs of sep 
aration appear in this portion of the boundary layer in 
the accelerating duct. 

The relation between the suction-surface boundary 
layer and the location of the passage-vortex roll-up 
may be instructive. In references | and 2 it was sug- 
gested that when nozzle blade suction-surface boundary 
layers remained thin, because of well-designed velocity 
profiles or removal of the accumulations by means of 
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radial flows, the passage-vortex formation might be 
delayed or even prevented. In the present case, cer 
tainly, the vortex formation did not occur where the 
cross-channel flow first encountered the suction sut 
face and its thin boundary layer but was delayed until 
the cross-channel flow had reached well up onto the 
surface spanwise. In the course of numerous visuali 
zation experiments on secondary flows by the present 
authors, the conviction grew that the passage vortex 
results primarily from deceleration and halting of the 
component of cross-channel flow along its path in the 
direction normal to the main-stream flow rather than 
from its encountering a wall or a thick suction-surface 
boundary layer or from suffering a sharp change in 
direction. This point of view appears to be supported 
rhe 


surface boundary-layer deflection indicates a negative 


by the flow patterns seen in Fig. 4. suction 
static-pressure gradient in the spanwise direction on the 
suction surface. The measured static pressures (Fig 


14, reference 10) prove it is so. In this region the 
cross-channel boundary-layer flows may be accelerating; 
at any rate, they are probably not decelerating. The 
the 


indicates 


leveling off of suction-surface boundary-layer 


region where the 


Where 
the cross-channel 


streamline arrival at a 


spanwise gradients are absent. this occurs, 


flow com 


there is deceleration of 
ponent in the direction normal to the main stream 
With no pressure gradient to urge it on normal to the 
main stream and having little momentum, the cross 
channel flow soon has nowhere to go but downstream 
in a direction parallel to its vorticity vectors, thus 
giving rise to the flow vortices. This purely mechani 
cal, continuity-preserving action provides a_ reason 
why passage-type vortices are sometimes seen to appear 
in mid-passage with no wall or thick boundary layer 
encountered. 


As discussed in references + and 7, the cross-channel 
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Fic. 5. Vortex roll-up at exit of accelerating duct 
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component actually is a “double” boundary-layer flow 
with two oppositely directed vorticities. Owing to 
viscous action, the vorticity near the wall is rapidly 
dissipated, the remaining outer half appearing as the 
passage vortex. 

Main-Stream Deflections ~The overhead views shown 
in Fig. 6 illustrate the deflection of the main stream 
around the region occupied by the passage vortex. For 
these photographs, the smoke was admitted near the 
inlet at one fourth the spanwise height. Deviations 
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Displacement of main-flow streamlines due to passage 
vortex formation 


from the calculated potential flow streamline shapes 
drawn on the base measure the secondary-flow effects 
on the main stream. A more dramatic illustration o 
the main-stream deviations caused by the secondary 
flow effects is shown in the sequence of pictures (Fig 
7) taken from a high-speed movie. For this movi 
the vortex shown results from smoke injected in the 
boundary layer by a probe at the inlet near the pres 
sure surface. For the main-stream flows, a spanwise 
traverse was made with the probe at the position shown 





in Fig. 6(a). The pictures were taken looking up- 
stream through the pressure surface as shown in Fig. 3 
As the traverse starts at mid-span and _ continues 
toward the mid-wall, the main stream can be seen to 
deflect around the passage-vortex region. It should 
be remembered that only one of the streamlines of the 
passage vortex appears in the picture. Because the 
vortex persists to the passage exit and beyond, the 
sizable displacement of the main stream likewise 


persists. 


Secondary-Flow Patterns in Three-Dimensional Ducts 


SMOKE TRACE OF MEAN ‘ 


STREAMLINE OF MAIN FLOW b.In actual turbomachines, the end walls of the ducts 
and passages vary radially, converging or diverging, thus 


causing radial turning of the main flow as well as the 
tangential ‘“‘two-dimensional”’ turning described above 

The three-dimensional ducts made of Lucite, shown 
in Fig. 8, were used to demonstrate some secondary 
flow patterns which may result from ‘‘three-dimen 
sional’ main-stream turnings. Smoke was _ injected 
rhe coer ren wer peace near the inlet of a duct [Fig. S(a)] which first acceler- 
ates and then turns the flow. The pressure surface was 
painted black (Fig. 9) to provide contrast for viewing 
the smoke. In the accelerating section, the main- 
stream direction near the end walls can be taken as 
roughly parallel to the curved end walls. The upper 
boundary-layer streamline along that region can be 
seen to overturn somewhat. A second boundary-layer 
340907 streamline, partially obscured in the picture, flowed 


_E1G. 6. Two-dimensional accelerating 90°-turning | duct. along in the corner formed by the end wall and the 
Top: (a) Main stream near suction surface at inlet. Bottom: . 


(b) Main-stream flow mid-span at inlet. pressure surface. Both streamlines fed into a sepa- 
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INVESTIGATION OF SI 


rated region along the pressure surface of the turning 
section. From this region came cross-channel flows 
on the end walls, resulting in the vortices seen near the 
suction suriace at the exit. A similar (qualitatively 
cross-channel flow and stagnant region can be seen 
more clearly in Fig. 10 for another duct. 

Fig 10 shows secondary-flow patterns in a three- 
dimensional duct [Fig. 8(b)] which first turns and then 
accelerates the flow. Here the suction surface has been 
blackened to provide background for viewing the smoke. 
In Fig. 10(a), the boundary layer along the suction 
surface overturns with respect to the main stream in 
of the duct. This kind of 


in a turbomachine duct 


the converging section 


end-wall can 


then be expected to give rise to radial overturning of 


convergence 


the boundary layer and radial flows on the pressure and 


suction surfaces. For Fig. 10(b), smoke was fed into a 
separated region on the pressure surface in the turning 
section. Cross-channel flow can be seen coming from 
this region, eventually to flow downstream along the 
suction surface of the duct. 

As pointed out in reference 10, there are two opposing 
actions occurring in the boundary layer as it moves 
from the inlet to the exit of a duct. The viscous action 
tends to thicken the boundary layer; the 


For the duct in Fig. 8(a) with a 


acceleration 
tends to thin it out. 
thin inlet boundary layer, the viscous effects would 
predominate (similar to no-spoiler case, reference 10) 
causing a net thickening of the boundary layer through 
the duct. This added to the increase due to separation 
on the rather small pressure surface at the turning sec- 
tion gives the total boundary-layer flow seen at the exit. 


In the duct of Fig. S(b), the early separation in the 
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Fic. 11 Coordinate axes for flow over surface 


turning section provides a large region of thick bound 
ary layer near the inlet. 
effects which tend to thin out the boundary layer pre- 
Which 
of two such ducts will produce less boundary layer at 


dominate throughout the accelerating section. 


the outlet depends, of course, on the relative magnitudes 
of the opposing effects and upon the relative sizes of the 
The duct of 
Fig. S(b) develops a large separated flow region, and 


separated regions in the turning sections. 


the boundary layer is then thinned down by acceler- 
ation. The duct of Fig. S(a) develops a boundary 
layer due to viscous effects and then a sizable separated 


flow over a small throat region near the discharge. 


ANALYSIS OF CROSS FLOWS 


As indicated in the Introduction, the assumption of 
laminar boundary-layer flow makes it possible to formu 
late solutions of the boundary-layer equations for cer- 
tain types of main-stream flow. In the following anal- 
ysis, solutions will be presented for boundary-layer 
cross flows similar to those shown in Figs. 2 and 3. 
Certain aspects of the solutions are checked experi- 
mentally by means of flow visualization. 

The equations for determining the steady incom- 


pressible laminar boundary-layer flow over a surface 


with coordinate axes oriented as in Fig. 11 are 
u(Ou/Ox) + v(Ou/OyY) + w(du/ds) — v(O°u/dy" 
X — (1/p) (Op/Ox) (la) 
u(Ow/Oxv) + v(Ow/Oy) + w(dw/Oz) — v(O°*w/Oy") 
Z — (1/p) (Op/Oz) (1b) 
(Ou/Ox) + (Ov/Oy) + (Ow/dz) 0 (Ie) 


The boundary conditions for the system of equations 
are 
uu v Ww 0 for y 0 
lim 1 U lim w WV 
The present investigation is concerned with two- 
dimensional main-stream flows defined by 


U Uy) (a constant) (2) 


For this configuration, the 
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Such flows are characterized by the fact that the: 


That 1S, th } 


entire streamline pattern can be obtained by translat 


streamlines form a system of ‘‘translates. 


ing any particular streamline parallel to the leadin 
edge of the surface (Fig. 12). The condition 
is considered a reasonable approximation for many | 
turbomachine configurations. 

For the type of main-stream flow described by E 

) 


2) and (3), the boundary-layer equations [Eqs 


Ib), and (le) | become 


u(Ou/ OX v(Ou,/ OV v(O-u/ OY () 
Ow Ow O-w te 
Uu + v — py > Uviain 
Ov Oy Ov 
(Ou/Ox) + (Ov/Oy 0 


Examination of Eqs. (4a) and (4c) and the boundan 
conditions reveals that these equations are identical | 
to the boundary-layer equations for flow over a flat 
surface in a uniform stream. The solution of thes 


equations is the well-known Blasius solution 


u UF’ (n 
r (1/2) V vUo/x (nF’ i | 
where n yV Uo/ wx 


and F(n) is the Blasius function (tabulated in Table | 
The boundary-layer velocity component w is now 


expressed as 

3. ax'P (n r 

where /; satisfies the ordinary differential equation 
a Lae, fo FP, + 1 Q \ 


with the boundary conditions 





PO 0) 
lim 2? (n | 

re-—ew 

pe lies aaial 
LEADING 
EDGE OF | 
SURFACE 

age 
Fic. 12. Streamline pattern as system of translates 


That is, th 


V translat 
he leading 
1) 


4c 
boundar 
identical 
a flat 
of thes 


ver 


Table | 


iS now 


ation 


WwW 
INES 


[INVESTIGATION OF 


rABLE | 
Values of Fand F 


It can then be verified by direct substitution that Eqs. 
5), (6), and (7) constitute a solution of Eqs. (4a), (4b), 
ind (4c Similar results were obtained in references 
OG; i. 


Numerical solutions of Eq. (8) were obtained for val 


7 and § for the special cases 7 
ues of 7 ranging from 1 to 10 on an IBM card pro- 

are 
(The 


function ?o(7) is not given as it is identical to the func- 


gramed electronic calculator. These solutions 


presented in Table 2 and plotted in Fig. 15(a). 
tion F’(m) listed in Table 1.) A description of the pro- 
cedure employed to obtain solutions to equations sim1- 
lar in form to Eq. (8) is given in reference 12. Part of 
the procedure for the present problem consists of replac- 
ing the boundary condition 

lim P,(yn) = 1 


Seal 


by a boundary condition at a finite value of 7. In con- 
formity with the findings of references 7 and 9 the 


boundary condition 
P(8 | 


was chosen. This is equivalent to assuming that the 
outer edge of the boundary layer is defined for practi- 
cal purposes by n = 8. 

Results similar to those given can be deduced from 
the analysis presented in reference 9. However, this 
particular analysis is of the small perturbation type and 
issumes small total turning of the main-stream flow 
ie, W< U, 


isSunnecessary. 


It can be seen now that this restriction 


Main-Flow Streamlines 


The streamlines of the main flow can be obtained 
rom the relation between streamline slope and velocity 


SEC 


ON 


DARY FLOWS I! BUECTsS 923 
components 
az We -_ 
> a l q 
dX l 
From Eq. (9), 
g constant + do*v + (a;" zZ 
( yi () 
where a,” a,/l 


As the coetticients in Eq 10) are arbitrary, a variety 


of streamline shapes can be constructed by the proper 


choice of a A particular shape (e.g., the contour 


of a channel wall) may be approximated by construct 


ing an approximating polynomial according to usual 
numerical procedures. In this respect, the class of 
flows becomes a useful tool in investigating the in 


fluence of various streamline configurations on bound 
ary-layer behavior. The tabulated values presented 
herein permit approximations by 11th-degree poly 


nomials. 


Boundary-Layer Streamlines 


The projection of a boundary-layer streamline in a 


plane normal to the leading edge of the surface (i.e 


parallel to the x,y-plane) may be obtained from the 
relation along the streamline 
dy/dx = v/u (1/2) V (vUo/x)/ UF’) (nF’ h 

1] 


Along a streamline, y is a function of x; hence, along 
a streamline, 7 is a function of x only. By differenti 
yv l to 


dx, there results a simple total-differen 


ating n vx and combining with Eq. (11 
eliminate dy 


tial equation 


dx/x —2(F’/F)dn (12) 


Thus, along a given streamline in the boundary layer, 


(13) 


[F(n 


By specifying the initial values x) and m, / (mo) can 


be obtained from Table 1, and the value of the constant 
stream function A can be calculated. For this value of 
K, F(n), and hence n and P;(n), can be obtained from 
the tables for all values of x along the streamline. 

With 
n for all points along the streamline, y can be com- 
puted from y nV vx/l 
the streamline in the x,y-plane, normal to the surface, 


a correspondence established between x and 
Thus, the projection of 


is easily obtained. 
The projection of a boundary-layer streamline in a 
plane parallel to the surface can be obtained from the 


relation 


12 WW 2 werd ee 
> : = > atx! Ti (9 14 
dx u Uol i 
P 
where J; PF’ 
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Values of P; and P;’ 
T T 1 
i ee T 4 PE | P; | Ps | Ps 
It ; | 1 7 
36/2.1971141]2.8579975] 3.4501246| 3.99534467] 4. 50568766] 4. 98663616)5.4492713) 5.89 +16.3 
2.09 2.708 13.25 |3.746 | 4.206 4.639 | 442 ale 
998 2.559 | 5-051 \3 497 | 3-908 4.292 4 4 4.99 324 
8968 }2.411 2.855 3.252 3.615 3. 95¢ 4.264 ‘ 4.838 
80 2.264 | 2.61) | 3.0 3.326 3.6 3.8683 4.13 4.3 } 
703 2.12 | 2.47€ lz 2s |} 3.04 3.29 3.513 3.716 9 
608 979 2.284 |2.545 2.772 2.974 3.156 $3 4 
514 841 2.103 }2.322 | 2.508 é 813 939 
1.421 1.706 1.927 |2.106 | 2.255 38 4a¢ ? F 
331 575 1.757 } 1.900 3 1 3 
243 448 1.594 }1 761 4 }1.684 914 333 
: 1.325 | 1.437 } 1.512 562 593 609 614 é 
3 206 1.286 } 1.355 36 }1.35 334 309 
392 093 1.143 }1.16 143 }1.114 | 7 33 
914 984 1.007 97€ 94] 894 64 
8 88 878 852 | .808 S4 69% | .624 53 
as 78) 757 71 6 58z 507 428 347 
‘ £86 643 58( SOS 424 339 251 63 
628 597 536 458 3 28 18 93 
é 513 436 346 48 4 
502 434 343 243 | .136 34 ? 
444 359 257 148 38 71 178 ~ 36 
388 290 178 062 53 64 27: 375 473 
336 25 105 016 33 4€ 354 4s ‘ 
286 65 038 86 : 318 ~.424 52 
9 23 149 268 3a¢ 484 8 
57 -.078 -.208 323 -.433 33 | - .62 1 
O 127 254 3 ‘ 574 |-.663 44 
§ 33 171 }-.297 42 4 607 691 é 
7 21 334 -.44° 544 | 633 -.7lé2 
46 8 366 “ts 568 }-.651 | é ’ 
ois -.393 - .49€ ses -.664 | 24 
13 |-.168 }-.415 |-.513 598 | .67% | 9 
39 | -.193 |--433 - .S2¢ 606 |-.674 | 84 
63 216 |-.447 1.58 61 673 | 3 
| + ee - 
8s 235 458 -.541 61 668 | 17 3 
104 252 | --465 }-.543 6 66 j-. 708 4 
2 ~ .266 | - -.469 [- 601 649 690 3 
37 .278 -.671 - .536 |-.593 63 673 ° 
2 47 |-.533 |- 582 é -.654 68 
163 29s 393 46 ‘ 60€ -.634 6: 
174 3 393 463 : 556 : 613 633 
183 308 392 456 SOS S42 ‘ é 
19 30 389 449 493 526 be 8S 
36 +30 | 385 44 | 46 $3 48 [: 
3 | -.379 43 4 493 12 S2¢€ 3 
e | -.306 | -.373 }-.419 Pee -.47 492 504 ‘ 
8 | -.303 365 | --4068 43 458 472 483 49 
9 |-.299 35 -.39 ry |-.453 4€ é 
| | cede Be , | al 
TABLE 2 (Concluded ) 
Values of P; and P;’ 
‘ — ——— 
n Py P, Ps P, Po Pe Po Ps Py Pio P} Ps Px P Pe Pe Ps Ps Pg Ps o 
2.45) 1.158) 1.279) 1.331] 1.356) 1.369) 1.376) 1.380) 1.381) 1.381) 1.381]-0.065/ -0.210 |-0.295 |-0.348/ -0.383/-0.407/ -0.422 /-0.433 |-0.440| 445 
2.50) 1.154) 1.266) 1.316) 1.359) 1.350) 1.356) 1.359] 1.360) 1.360) 1.359] -.070] -.210] -.290] -.339} -.371] -.391] -.405] -.414] -.420] -.424 
2.55} 1.151) 1.258) 1.302) 1.322) 1.332) 1.337) 1.339) 1.340] 1.339} 1.339] -.075] -.209] -.284 --329| -.358| -.376) -.388) -.395/ -.400/ -.403 
2.60) 1.147) 1.248) 1.288) 1.306) 1.315) 1.319] 1.320] 1.320] 1.320) 1.319] -.079] -.207] -.278 ~ +519) -.345|] -.361| -.371] -.377] -.381) -.383 
2.65] 1.143) 1.237) 1.274) 1.290] 1.298) 1.301) 1.302) 1.302} 1.301} 1.300] -.082] -.205] -.271] -.309] -.331] -.345] -.35 -.359/] -.362| -.363 
2.70) 1.139) 1.227) 1.261) 1.275) 1.282) 1.284) 1.265] 1.284) 1.284) 1.283] -.085|} -.202] -.263/] -.298) -.318] -.330] -.338 342 ~.344| ~.345 
2.75) 1.134) 1.217) 1.248) 1.261) 1.266) 1.268) 1.268) 1.268) 1.267} 1.266} -.087| -.199]| -.256] -.287} -.305| -.316] -.322 -325 |] -.326] -.327 
2.80] 1.1350] 1.207) 1.235) 1.246) 1.251] 1.253] 1.253] 1.252) 1.251) 1.250] -.089] -.195]| -.248] -.276| -.292] -.301| -.306 - 309 -310; -.310 
2.85} 1.126) 1.196) 1.223) 1.233) 1.237) 1.238) 1.238) 1.237) 1.236) 1.235] -.090] -.191] -.240] -.266] -.280] -.287| -.291 293 | -.294; -.293 
2.90) 1.121) 1.188) 1.211) 1.220) 1.223) 1.224) 1.223) 1.223) 1.222) 1.221] -.090} -.187]| -.232] -.255] -.267| -.274| -.277 27 - 278 | -.276 
2.95) 1.117] 1.179) 1.200} 1.207) 1.210) 1.210) 1.210) 1.209] 1.208] 1.207} -.091] -.182} -.223] -.244] -.255| -.260|] -.263 264 -.264| - .263 
3.00) 1.112) 1.170) 1.189) 1.192) 1.198) 1.198) 1.197)1.196/ 1.195) 1.194] -.091] -.177] -.215] -.234] -.243] -.247] -.249 250 | - 249 | ~.249 
3.10} 1.103) 1.153) 1.168) 1.1735) 1.174] 1.174] 1.174)1.173) 1.172) 1.171] -.090] -.166] -.196] -.213] -.220] -.223] -.224 +224] -.223) -.222 
3.20) 1.094) 1.137) 1.149) 1.153] 1.154) 1.153) 1.152)1.151} 1.151) 1.150} -.087] -.155] -.182] -.193] -.198] -.200/ -.201 -200 | -.199] -.198 
$3.30) 1.085] 1.122) 1.152) 1.134) 1.1355) 1.154) 1.133)1.133) 1.152) 1.131] -.084] -.143] -.166] -.175| -.178] -.179] -.179 178 |.- wie -.176 
3.40} 1.077) 1.108) 1.116) 1.118] 1.118) 1.117] 1.116/1.116) 1.115) 1.114] -.080] -.132] -.150] -.157] -.160/ -.160] -.159 159} -.158/ -.157 
$.50} 1.069] 1.095; 1.102/ 1.103) 1.103) 1.102] 1.101/1.101/ 1.100} 1.100} -.076/ -.121] -.136]| -.141/ - 142 | -.142| -.142 141] -.140] - 9 
3.60} 1.062] 1.084) 1.089} 1.090} 1.089) 1.089) 1.088/1.087) 1.087) 1.086] -.071] -.110] -.122] -.126; -.127| -.126] -.126 125] -.124/] -.123 
3.70} 1.055] 1.073) 1.077/ 1.078] 1.077] 1.077] 1.076/1.076} 1.075} 1.075] -.066| -.099] -.109] -.112} -.112] -.112] -.111 -110 | -.109] -.109 
$.80}] 1.049) 1.064) 1.067/ 1.067] 1.067] 1.066] 1.066]1.065/ 1.065} 1.065] -.061] -.089] -.097]| -.098/ - 099 | -.098/ -.098 -097 | - 096 | - .O96 
3.90] 1.043) 1.056) 1.058) 1.058/ 1.058} 1.057) 1.057/1.056/ 1.056) 1.056] -.056} -.080]| -.086/ -.087} - 087 | -.087| -.086 -085 - .085 | - .084 
4.00] 1.038] 1.048] 1.050] 1.050} 1.049] 1.049] 1.049/1.048] 1.048) 1.048} -.051] -.071] -.076] -.077!] - 076 | -.076| -.075 -O7S | -.074| -.074 
4.10) 1.033) 1.041) 1.043) 1.043} 1.042) 1.042) 1.042/1.041/ 1.041) 1.041] -.046] -.063 | -.067]| -.067|} -.067| -:066] -.066 -O65 | -.065| -.064 
4.20] 1.028] 1.035) 1.036/ 1.056] 1.036} 1.056) 1.035/1.035/ 1.035} 1.035] -.041] -.055 | -.058] -.059/ - 056 | -.058/ -.057 -0S7 | -.056| -.056 
4.30] 1.025) 1.030) 1.031/ 1.031] 1.031) 1.030] 1.030 /1.030] 1.030} 1.030] -.037] -.049 | -.051 | -.051] -.050/ -.050] -.050 -049 | -.049| -.049 
| | 
! 
4.40} 1.021/ 1.026) 1.026/ 1.026] 1.026} 1.026] 1.025/1.025/ 1.025} 1.025] -.033] -.042 | -.044] -.044]| - 046 | -.043|] -.043 -043 | -.042| -.042 
4.50} 1.018) 1.022) 1.022)1.022/ 1.022] 1.022)1.022/1.021]}1.021/1.021] -.029] -.037 | -.038 | -.038] - 038 | -.037| -.037 037 | -.037| -.036 
4.60} 1.015) 1.018) 1.019) 1.018) 1.018) 1.018/1.018/1.018]1.018/1.018] -.025] -.032 | -.033 | -.033| -.032/ -.032] -.032 032 ~-031 | -.031 
4.70) 1.013) 1.015; 1.016/1.015}1.015/ 1.015}1.015/1.015}1.015/1.015] -.022] -.027 | -.028 | -.028] -.026| -.027] -.027 -027 | -.027]| -.027 
4.60} 1.011] 1.013) 1.013) 1.013] 1.013} 1.013/1.013/1.013)1.013}1.013] -.019] -.023 | -.024] -.024/ - 028 | -.023| -.023 -023 | -.023] -.023 
+ 4.90) 1.009) 1.011) 1.011/1.011)1.011/1.011]1.010/1.010}1.010)1.010] -.016] -.020 | -.020 | -.020] -.020] -.020/ -.020 020 — -.019 
§.00/ 1.008} 1.009/ 1.009/ 1.009} 1.009/ 1.009}1.009 |1.009]1.009/1.009} -.014] -.017 | -.017 | -.017/ -.017] -.017] -.017 017 | -.016]} -.016 
5.10} 1.006} 1.007} 1.007/1.007}1.007/1.007}1.007 |/1.007]1.007/1.007] -.012/ -.014 | -.014] -.014] -.014] -.014] -.014 014 | -.014/ -.014 
§.20/1.005/ 1 .006/ 1.006/1.006/1.006/ 1.006] 1.006 |1.006}1.006}1.006 | -.010/ -.012 | -.012 | -.012] -.012] -.012] -.012 012 | -.012| -.012 
5.30] 1.004/ 1.005/ 1.005/1.005/1.005/1.005/1.005 |1 .005}1.005/}1.005 | -.009] -.010 | -.010 |] -.010] -.010] -.010] -.010 010 | -.010] -.010 
5.40} 1.004/ 1.004) 1.004/1.004/1.004]1.004/1.004 |1 .004/}1.004/1.004 | -.007/] -.008 | -.008 | -.008| -.008/ -.008| -.008 -008 | -.006 | -.008 
5.50] 1.003} 1.003] 1.003} 1.003/1.003}1.003/}1.003 |1.003}1.003)1.003 | -.006| -.007 | -.007 | -.007] -.007 | -.007] -.007 -007 | -.007 | -.007 
5.60] 1.002} 1.003} 1.003 /1.003/1.003/1.003/1.003 }1.003}1.003/1.003 | -.005] -.006 | -.006 | -.006/ -.006 | -.006/] -.006 -006 | -.006 | -.006 
5.70] 1.002} 1.002/ 1.002/1.002/1.002]1.002/1.002 |1.002/1.002/1.002 | -.004] -.005 | -.005 | -.005]| -.005 | -.005] -.005 -005 | -.005 } -.005 
5.80] 1.002] 1.002] 1.002/1.002/1.002}/1.002/1.002 |1.002/1.002/1.002 | -.003]| -.004 | -.004 | -.004] -.004 | -.004] -.004 -004 | -.004 | -.004 
§.90/1.001/ 1.001} 1.001/1.001/1.001/1.001/1.001 j1.001/1.001}1.001 | -.003] -.003 | -.003 | -.003] -.003 | -.003] -.003 -003 | -.003 | -.005 
6.00/1.001) 1.001} 1.001/1.001/1.001/1.001}1.001 /1.001}1.001/1.001 | -.002] -.002 | -.002 | -.002] -.002/ -.002] -.002 -002 | -.002 | -.002 
6.10/1.001/ 1.001} 1.001/1.001}1.001/1.001 /1.001 /1.001/1.001/1.001 | -.002/ -.002 | -.002 | -.002/ -.002 | -.002/ -.002 -002 | -.002 | ~-008 
6.20/1.001/ 1.001) 1.001/1.001/1.001/1.001/1.001 |1.001/1.001/1.001 | -.001 | -.002 | -.002 | -.002); -.002 | -.002] -.002 -002 | -.002 | -.002 
6.30} 1.000} 1.000} 1.000/1.000/1.000/ 1.000 /1.000 |1 .000/1.000/1.000 | -.001 | -.001 | -.001 | -.001] -.001 | -.001] -.001 -001 | -.001 | -.001 
| | | | | | | 0) ) ) ) ) ) ) ) Oo} o | 
6.40 | l | | | l | , 4 
6.50 | | 
6.60 | | | } | 
6.70 | | 
a See Senne es Pawnee: amide = 
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Numerical integration of Eq. (15) will then yield the 
shape of the streamline projection in the plane parallel 
to the surface. It is also possible to express the integral 


With a correspondence established between x and 7 
along a streamline by Eq. (13), the right side of Eq. 








e-pze FF 14) is expressible as a function of x alone, and it follows 
‘ ese , that in Eq. (15) as a function of 7 by substitution of Eq. (13 
5.800 t - . into the integrand. The equation for z then becomes 
051 Z > a;*x' I (n)dx + 2 (15 : ; 
és = Jn 
a1 72/ Eqs. (13) and (15a) can then be considered as para 
34) metric equations for the boundary-layer streamlines 
°01 - with 7 as the parameter. 
27 The actual three-dimensional boundary-layer stream 
47 line paths are then determined from a combination of 
5.6) the two projected views. The boundary-layer velocity 
22 components u, v, and w are computed directly from 
- } 4.8f Eqs. (5), (6), and (7). 
92 In experimental studies of three-dimensional bound 
ey ary flows over flat surfaces, it is evident that the 
boundary-layer streamlines near the surface approach 
Theoretical determination of this 


a limiting shape.” 
‘limiting streamline’ is possible for the present problem 


by integrating Eq. (15) at y 0. This vields 
m a.*yi4 
z, } J (0) + 2 16 
i 1+ | 
The value /,(0) P,(0)/F’(0 0/0 can be determined 
by use of L’ Hospital's rule 
lim [Pi(9)/F'( P'(0)/F"(0 


n—) 








The P,;’(O) values are obtained (see Table 2) during 


Eq. (10 The quantity 





the numerical solutions of 
F"(0) is 0.33206. 

By inspection of Eq. (16), the shape of the limiting 
deflection boundary-layer streamline depends only on 
the geometry of the main-stream flow and is inde 
The other boundary 


pendent of the inlet velocity | 
(13) and 


layer streamlines are functions of » [Eqs 
(15a) | and thus are not independent of / 


Normal and Tangential Velocity Components 


ae 


In order to obtain a better physical picture of bound- 
ary-layer flow in relation to main-stream flow, it is 
often desirable to resolve boundary-layer velocities into 
components tangent and normal to the main-stream 
Denoting the tangential component by 


flow direction. 
¢ and the normal component by 2, it can be shown 


(see Appendix) that 


(IS 


‘004 | ; 
: ‘ - where 7 Vl + |] 











“O02 0) 4 #6 |2 16 20 2.4 a ee 
J] Hi-10 
(B) H; (7) and * denotes division by The function // is 


Fic. 13. Function P;(y) and H,(n plotted in Fig. 13\b 








226 JOURNAL OF THE AERONAIT 


APPLICATION OF ANALYSIS 


Calculation of Functions 


Solutions of Eq. (S) for various values of 7 have been 
presented in several publications. References 7 and S 
| with a tabular listing of the 


Refer- 


present solutions for 2 
functional values given in the latter reference. 
ence 9 contains solutions for 7 0, 1, 2, 3, 4, and 
—1/2 in graphical form. In the present report, nega 
tive values of the exponent 7 have not been investi- 
gated because they require infinite velocities at the 
Where 7 0, Pi aide 


leading edge of the surface. 


Flow with Circular-Arc Streamlines 


As an application of the theory, the boundary-layer 
development for a flow with circular-are streamlines 
was investigated. The specific flow region analyzed is 
restricted to that shown in Fig. 14. As seen in the 
Figure, the flow is normal to the leading edge of the 
surface and crosses the downstream boundary of the 
region at 45 

Polynomial Approximation The equation of the 


family of circular ares depicted in Fig. 14 is 

(19) 

where c is an arbitrary constant and0 S XK S VY 2/2. 
The right side of Eq. (19) can be approximated by the 


polynomial 
(20) 


This polynomial then serves as an approximation to 
the streamlines of the main flow and is used to apply 
the theory to the determination of the boundary-layer 
characteristics. 
Boundary-Layer 
boundary layer are obtained from simultaneous solu- 
In the 
present example, three representative streamlines were 
0.05 and at 
positions corresponding, respectively, to 1/4, 1/2, and 
3/4 of the The initial 
value of x was chosen at 0.05 instead of at the surface 


Streamlines-—Streamlines in the 


tion of Eqs. (13) and (15) previously described. 
computed beginning at the position x 


boundary-layer thickness. 
leading edge because of the singularity in the flow equa- 


tions at x = 0. From Table 2 and Fig. 13 it is seen 
that a value of 7 which reasonably defines the ‘‘outer 
edge’ of the boundary layer is 6.4. Accordingly, the 
values n = 1.6, 3.2, and 4.8, respectively, correspond 
to positions at 1/4, 1/2, and 3/4 of the boundary-layer 
thickness. 

A plot of these streamlines including the limiting 
boundary-layer streamline determined from Eq. (16) 
is given in Fig. 15. These streamline paths are similar 
in form to the boundary-layer overturning on the end 
walls of the circular-arec sheet-metal cascades (as 
visualized and presented in reference 1) before passage- 
vortex roll-up occurs. 

Velocity Profiles 
ponents n* and ¢* at x = 0.2, 0.4, and 0.6 are presented 


It is interesting to note that the * profiles 


The profiles of the velocity com- 


in Fig. 16. 


ICAL 


SCIENCES MARCH, 


have maximum values for 7 < 1.5 that is, in the lowe- 
fifth of the boundary layer. An examination of ¢] 
[Fig. 13(b 


in general, if the coefficients in the polynomial appr 


curves for // and Eq. (17) indicates that 


mation of the flow streamlines a;* are positive, th 
n* will always have a maximum value for 7 1.3 

reference 6 This can be directly associated with t} 
discussion on vorticities of the cross-channel boundar 
layer flow in blade rows.' In this reference, the cross 


channel ‘‘double’’ boundary-layer flow is consider 


as divided into two regions. In the region near th 
wall, 2* rises from zero at the wall to some maximy 
value in the sublayer of the boundary layer. In th 
upper boundary-layer region, * declines from th; 
maximum value to zero at the main stream. Thes 
two distributions of velocity may be associated wit 
vorticities of opposite signs in the two parts of th 
boundary layer. It is further hypothesized that sue! 
vorticities generate the actual flow vortices observ 
in blade passages.‘ However, visualization studies ; 
various blade passages generally indicate only simp] 
passage-flow vortices. Usually there is either no ey; 
dence at all of the formation of a flow vortex fro 
boundary-layer cross flow in the sublayer adjacent t 
the wall or, at best, a formation that is small and poorh 
defined. It is therefore conjectured that the narrowness 
of the boundary-layer sublayer and the large viscous 
forces in that region tend to prevent the sizable roll-y 
characteristic of vortex formation in the upper bound 
ary-layer region. The position of the maximum valu 
of n* in Fig. 16 and the shape of the //; curves [Fig 
13(b) |] both substantiate this conjecture. 


Overturning in Thick and Thin Boundary Layers for 60 

Bends 

In a theoretical and experimental investigation, 
comparison was made between the overturning of 
relatively thick and a thin boundary layer for 60) 
main-stream turning. 

Theoretical Comparisons of Overturning 
cal comparison of thick and thin boundary-layer over 
turning was made by determining the overturning 0! 


A theoreti 


the boundary layer near the surface for two main 
stream flows [consisting of streamline translates de 
scribed by Eq. (10) ]. 

The flow field whose main-flow streamlines are repre 


sented by 


g (k*/6)xi + 2 O<-x 


lA 


in the thick boundary-layer case is shown in Fig. 17 
The value &* is taken equal to 7/3 for 60° main-stream 


turning. The boundary-layer development in_ th 
long nearly straight inlet section provides a_ thick 


boundary layer at x = 0.5, the region where appreciable 
turning begins. 

In the thin boundary-layer case, with streamlines 
represented by 
s = (k*/6) (x + 0.5)§ + 20 0<x<05 (2 


in Fig. 18, the appreciable turning begins right at the 
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15. Boundary-layer streamlines for circular-are flow 
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16. Velocity profiles of »* and ¢* for circular-are flow 
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17. Main-flow streamlines for thick boundary-layer flow 
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18. Main-flow streamlines for thin boundary-layer flow. 
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19. Boundary-layer streamlines for thick boundary-layer 
flow. 
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streamlines for thin boundary-lavet 


flow 


Boundarvy-laver 


leading edge with no chance for prior boundary-layer 


build-up. 
identical. The boundary-layer streamlines are pre 
sented in Figs. 19 and 20. 

The thicker boundary layer (Figs. 17 and 19 
turns more rapidly than the thinner boundary layer 
(Figs. 1S and 20). 


that thicker boundary layers with lower viscous shear 


over 
Physically, it would be expected 


forces would have sharper cross-channel overturning 
than thinner boundary layers with correspondingly 
The results 
Comparable results were ob 


higher viscous shear forces. theoretical 
agree with expectations. 
tained in experiments* in turbine nozzle blade con 
reference 2 demonstrate 


figurations. The results of 
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Smoke-flow visualization of main streamline in 60 
turning rectangular bend 


Fic. 21 


ICAL SCI 





The turning sections for the two cases are 
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the larger response of thicker turbine blade boundan 
layers to radial pressure gradients than thin blag 
boundary layers. The fundamental question remajgs 
as to whether or not the theory presented herein has 
any immediate demonstrable application to actu, 
flow behavior other than in a purely general sense 
This problem was investigated experimentally by meap; 
of flow visualization. 

Experimental Comparison of Overturning—The maip 
flow, which is a system of streamline translates, can by 
considered as the flow through a_ two-dimensional 
cascade of a large number of very thin blades having 
the main streamline shape with vanishingly small blag 
spacing. As in axisymmetric flow theories, the blades 
supply the required force field [Eqs. (1)] for the pre 
scribed flow turning. For such channels with ven 
small spacing, the turning is very close to the stream 
line-translates turning. Thus, an ideal one-dimen 
sional analysis indicates that, on the basis of continuity 
and turning considerations, the flow velocities should ap 
proximate the theoretical translate flow velocities quite 
closely. 

An experimental investigation was made to deter 
mine if the theory provides a reasonable approximation 
to the boundary-layer overturning in a channel wit 
spacing more nearly comparable to that found in tur 
bomachine configurations. The experimental deter- 
mination and comparison of thick and thin boundary 
layer overturning was made by means of smoke-flow 
visualization (apparatus and procedures described in 
reference 1) in two-dimensional rectangular bends 
with channel shapes the same as the previously de 
scribed main streamlines. 

With a maximum Reynolds Number of about 6 X 
104, the flow for these tests was well within the laminar 
shows the main-stream-flow paths 


As expected, boundary- 


range. Fig. 21 
traced by means of smoke. 
layer accumulations near the suction surface cause the 
from the desired 


to deviate somewhat 


Nevertheless, by keeping the solid 


main stream 
streamline shape. 
ity of the turning section high (approximately 2.0), 
sufficient main-stream guidance is provided to maintain 
the mean streamline close to the desired shape. 

The overturning of the end-wall boundary layers as 
visualized by means of smoke for both the thick and 
thin boundary layers is shown in Figs. 22(a) and 
2?(b), respectively. The smoke was introduced through 
static taps at the inlet to the turning sections of the 
bends. Because of the obvious difficulties in accu 
rately locating streamlines at various heights within 
the boundary layer and for the more important reason 
that only the limiting line is independent of inlet ve 
locity, attention was devoted mainly to obtaining the 
of the flow along the end wall 


deflection 


limiting line (y 0 
The experimental 
(Fig. 22) corresponds closely to the theoretical pre- 
boundary-layer 


cross-channel limiting 


dictions, except near the region of 
accumulation at the suction surface. 
Because the smoke traces always appear faint m 


overhead photographic shots (see reference 1), the pic 























boundan 


hin blac 


n remains 


1¢ rein has 


to actual 


ral Sense 


by means 


The maip 
eS, Can be 
mensional 
es having 
nall blad 
he blades 

the pre 
vith very 
e stream 
1e-dimen 
ontinuity 
hould ap 
ties quite 


to deter 
ximation 
nel with 
d in tur- 
il deter- 
yundary- 
oke-flow 
ribed in 
r_ bends 
usly de- 


ut 6 X 
laminar 
v paths 
undary- 
iuse the 
desired 
le solid- 
ly 2.0), 


laintain 


LVers as 
ick and 
a) and 
through 
; of the 
1 accu- 
within 
reason 
ilet ve- 
ing the 
d wall 
flection 
al pre- 
y-layer 


unt in 


he pic- 








rCONDARY F4L 






PRESSURE SURF ACE 








PRESSURE SURFACE 
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SUCTION SURFACE 


Fic. 22. Smoke-flow visualization of limiting (y = 0) cross 
channel deflection of boundary layer. Theoretical predictions 
1own as dotted lines Top: (a) Thick boundary-layer limit 


ing deflection Bottom: b) Thin boundarv-laver limiting 


deflection 


tures of the thin boundary layer presented in Fig. 25 


were taken looking upstream and can be compared 
The photo 


graphs indicate that, in the region near the wall, the 


with a simultaneous overhead photograph. 


theory predicts overturning very satisfactorily. 








Boundary-Layer Velocity Profiles. The velocity pro 


and ¢* for the thick and the thin boundary 


layers are presented in Figs. 24 and 25, respectively. 


file S of n* 


rhe values are plotted at equally spaced intervals of 
Within the turning sections of the bends and at corre- 
0.5, 0:6, 0.7, 0.3, 0.9, 
1.0 for the thick boundary-layer case and at x 
0.1, 0.2, 0.3, 0.4, and 0.5 for the thin boundary-layer 
The n* and the ¢* profiles peak at lower 7 values 


sponding locations (i.e., at x 
and 

















OWS IN DUCTS 220 









laver For 





for the thick than for the thin boundary 


the thick boundary layer (Fig. 24), »* becomes larger 
than | at stations far from the leading edge (. 0.9 
and 1.0), indicating that the component of the bound 


stream-flow direction 
For both 


the main 
the inlet v¢ 


arv velocity normal t 
has become larger than locity 
thick and thin boundary lavers (for a sizable portion of 
the tar il component ¢* at 


laver), genti 


the boundary 
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Smoke-flow visualization of 


thin 


Fic. 23 
cross-channel 


VWiddle: b 


limiting 
boundary-laver flow Tot i 


Downstream view f boundary 


view 
streamline near pressure surface at inlet Bottom: (c) Down 
stream view of boundary-layer streamline near mid-channel at 
inlet 
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n* t* 1.6 











Velocity profiles of 2* and ¢* for thick boundary-layer 
flow 


Fic. 24. 
sufficiently large x. actually becomes greater than the 
local main-stream velocity. This is illustrated by the 
t* profiles which overshoot and then approach the 
main-stream 7* values from above. 

This theoretical result, likewise observed in reference 
7, may be explained as follows: For flows such as are 
considered herein, where 


W = Day 


it follows [Eqs. (1) | that 


1 Op 
p O2 


Z— Uy > ta x' 
a0 
that is, there is a gradient of the external force field 
and the static pressure in the z-direction. With Uy 
constant everywhere along the leading edge, there then 
exists a total-pressure and energy variation along the 
leading edge. By convection, the more than main 
stream turning (overturning) of the boundary layer 
transports boundary-layer material from a region of 
high total energy into regions of lower main-stream 
total energy. large 
enough energy gradients at the inlet, it is possible for 


For sufficient overturning and 


regions to exist in which the boundary layer has a higher 
total energy than the local main stream does. Because 
the boundary layer is assumed to have the same static 
pressure the main stream, the boundary-layer 
velocities must exceed the local main-stream velocities, 
For these flows where veloci- 


as 


as illustrated in Fig. 24. 
ties in the boundary layer exceed the local main-stream 
velocity at some axial distance x from the leading 
edge, there will be boundary-layer velocities exceeding 
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the local main-stream velocities everywhere downstreap 


of that axial ocation. 


CONCLUSION 


Experimental verification of boundary-layer over 
turning for flows in a bend was obtained in this invest; 
gation for the laminar boundary-layer behavior as pr 
dicted by a theory based on similarity consideration; 
and main streamline translates. A main streamlin 
through the bend conforms quite closely to the desire 
flow path. Actually, with the great range of possibl 
flow paths describable by an 11th degree polynomial 
the requirement of main streamline translates (corre 
sponding to the assumption of axial 


turbomachine design) is not very restrictive. The 


Ssyimmetry in] 


For 


close agreement obtained between theoretical and ex.-| 


perimental boundary-layer overturning indicates that| 


the theoretical assumptions can be considered to be| 


reasonable for this laminar cross-channel flow theory 
In addition, experimental investigations similarly con 


ducted on other configurations very much like these} 


have provided qualitatively correct information con 
cerning the turbulent boundary-layer (overturning 
behavior for flows at operational speeds in turbomachin¢ 
configurations. These results indicate that the laminar 
boundary-layer theory may provide qualitative infor 
mation concerning turbulent boundary-layer behavior 

It is interesting to note that, while the theoretical 
development is based on boundary-layer velocity 
profile similarity in the x- and z-directions, the velocity 
profiles in the z-direction are similar only for a flow 


described by 
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Velocity profiles of n* and ¢* for thin boundary-laye! 
flow. 
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US INVesti-} avery For such flows the solutions have been ob yr wen tg 

10r aS pre-f yained here. The boundary-layer components in the 

iderations) direction are given by 

streamline) 

he desire w= dw > WP Ye ax'P W STREAMLINE 

1 possibk i zef(x 

lynomial For the flows considered in this report, the profiles 

= aes the boundary-layer velocity components in the : 

ey direction are all similar and are related by the Blasius Fic. 26. Resolutions of velocities into components normal and 
ve. The none : tangential to main flow 

l and ex ve , 

ates that | u UF '(n 5 ™ 

Ee uf i , ry " > * ‘ 
oe " nd, thus, are identical to the profiles obtained for ' T (7 t W* DO WAP l \Y 


Vv theory, | é : 
larly con. straight laminar flow over a flat plate. 
a Substitution of /7; + F’ for P;in Eq. (AQ) yields 


ike these | 
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on “—- ar (7 ] we Owen,)/T* (A10) 
> laminar Let 2 f(x) be the equation of a streamline in the ) 
ve infor w-plane. Then, from Fig. 26, ae — : 
P " ” Finally, substituting Eq. (AS) into Eq. (A10) gives 


ehavior ; 
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An Analysis of the Effect of Thermal Contact 
Resistance in a Sheet-Stringer Structure? " 
Antony D. Barber,* Jerome H. Weiner,** and Bruno 
Institute of Flight Structures, Columbia University, New Y 
September 21, 195¢ 
SUMMAR\ ( 
Ihe effect of thermal contact resistance in a typical sheet 1 
ture under transient heating conditions is analyzed, and numerica 
are presented. It is found that, for pulses of short duration, even a « 
resistance of small magnitude will produce a considerable effect 
perature and stress distributions 
DISCUSSION 
Qeverat recent experimental investigations have determi 
the magnitude of the thermal contact resistance in vari 
tvpes of aircraft joints (additional references are listed i 
reference 4 It is the purpose of this note to investigat« 
Contents cally the effect of such resistance on the temperature and stress IG. - 
Fi distributions in an aircraft structure under transient heating 
: } ; t This work is part of an investigation sponsored by the Wright Air D 
An Analysis of the Effect of Thermal Contact Resistance in velopment Center, U.S. Air Force, under Contract AF 33(616)-207 it 
a Sheet-Stringer Structure . * Research Assistant, Department of Civil Engineering and Engineer gatiol 
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cos v 
/ r } 7 
sinh 7 \ coth 4 + ¢ th ( \ ) 
h HW } 
LL 
rhe temperatures |; and J». are then determined by use of 
the inversion integral. For the special case of a step heat input 
f(r) = 0, 7 O; f(r) = 1,7 > 0] and b/a 1, the result is 
Contact h 
cturej Vi(g, tr ( ) r 
/ ? fl { 
py ‘ <> ™* COS ay l 5 
“= 1 an? COS a Y H] | l 
VE, 7) = 7 hy /he)\ t, 7 
67 T T T T 7 where H, = (La/RK 1/hy) + (1/h 5a 
neer 1 1 c " 
V ind a, are the roots of the equation 
na ? r - 
; id 4+ t + ay, tana, = H{l1 4 Ny / Ne H 5b 
T2.25 It is shown in reference 4 that the above equations simplify 
in the limit to the forms required for the special cases of perfect 
2: , + + t insulation, H, = O, and of perfect contact, H Results of 
oe computations for H; = 10 are presented for illustrative pur 
faster E poses in Fig. 2 
listed oe ° — . ~ 
sais 0! 1 | 1 J Ihe infinite series in Eq. (5) is slowly convergent for small 
re analy . values of +. An approximate solution for short times may be 
ind stres 1G.2. Temperature distributions for a step heat input, 7, = 10, ; , tee ; ; 
| Sail obtained by simplifying the Laplace transforms of the temper 
eee lane : iture in a manner valid for large value of p as follows (for b/a = 
l 
ht Air D - : ; 
71 itions. This note contains the principal results of this investi (¢, p) = [He Py pir 1. 4) Pe 
ngineer gation; further details may be found in reference 4 P(t. bp) = (1/p - (hi/ke) V ») | (6 
. . . . . sf f « ‘ . ’ ed 
rhe idealizations made in the formulation of the problem are 
similar to those introduced in the Appendix of reference 5 and rhe inverse transforms of the above expressions are then* 
n reference 6 and lead from the typical structure of Fig. l(a) to 7 
‘ — 7 , A s 
the one-dimensional model of Fig. 1(b). The temperatures 7 Vit. r) = erfc a 
ind 7, in the web and skin are then the solution of the following hy + he J 2% 
equations: t \ | 
ute , exp Mt + Hy?¢) erfe ( + Hi Ve | f(r a \4 
x (077 /Ox?) = O7)/Ot; box <0 la QV s 
0?7>/0x?) + [g(t)/hepc] = OT>/Odt; 0<z< 4 (ib Vt, r) = f(g )d¢ (hy /he) 1 I 
1 . ° = . « 
subject to the following boundary and initial conditions: 
— . : ie ; or the se heat i it given b 
Kh, |O7,(0, t)/Ox] = (L/R) [T2(0, t TO, t Qa For the pulse heat input given 
Kh, (07\(0, t)/Ox] = Khe[o7T.(0, t)/ox 2b ee Ree , 8 
OT;(—b, t)/Ox = OT>2(a, t)/Ox = 0 2c the result at & = Ois 
f, (x, 0) = 0 bx <0; T(x, 0 v7S25 -d ) j T Hy 703% Hy2r 
V4(0, r) = T ? ( ) : 
vhere A is the thermal conductivity, p the density, c the specific h, + ho | T H,?r l 
| leat, and x = K/(pc) is the thermal diffusivity. The contact a ‘ Hiro3!2 ‘ 
esistance is denoted by R and has typical units (ft.* hour °F of ) H i H x 
B.t.u., L is the effective length of contact, and g is the heat in j - ; 
ss I].2_,2 
put _ ‘ (9 
afi efi? erfe Hi Vr : 
itis convenient to introduce the following dimensionless quan T H,?r l 
tities H,?197( HM? 70? + 2 Hr Vr | 
i | T\Kh2/a?Qo; V2 = T2Khe/a?Q r = xt/a? om H,\?ro + 1)? V/ x (Hi2r9 + 1 { 
| & = x/a; H = La/Rkh,; f(r) = g/Qo : 
O,r) = mo }€ l r/tT li — (hy/h V,(0. + 
ere the constant Qo is an arbitrarily chosen heat input 
he solution of this boundary-value problem was obtained by where the function 
the method of Laplace transformation; the Laplace transform 
fa function with respect to time is denoted by a bar 
ir co! . tei : jm ; ri , ' * Transform pair #812, Campbell, G. A., and Foster, R. M., Fourter 
rhe transformed temperatures, V7, and Vs, are readily found ‘imal tp Meads Seatirations: DB. Ven Méectranl Commenr: tes.. Hn 
» be York, 1948 
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Fic. 3 (left 
contact resistance and rapidity of heating on maximum temperature difference across the joint, h;=h 


perature difference which would occur if 
of heating on the maximum stress, h; =h- 


mmay be expressed in terms of the error function for complex 
arguments. * 

Some numerical results based on Eqs. (9) are presented in Figs 
3-5, which may be used for an estimate of the importance of the 
that, 


resistance of small 


contact resistance. For example, it appears from Fig. 4 


for pulses of short duration, even a contact 


magnitude—i.e., large value of /1,;-—will produce a considerable 


effect on the temperature distribution 
l(a 


The axial thermal stress, ¢,, in the structure of Fig 


cor 
responding to the temperatures calculated above is, 


syvininetry, 


where 


(is | Tidy 4 ie . dx) (hy b+ hea 
= lef g(t) at | pc(bh, + ahz) (10a 


Inspection of Fig. 2 indicates that the maximum tensile stress 


occurs in the web (x = —) and the maximum compressive stress 


at the point on the skin midway between webs (x = a). For 


the special case of a = b and h, = he, these stresses will be 


numerically equal. This case is considered in Fig. 5, in which 
is plotted the variation with 7») and #7, of the ratio 
Ce o; = [Tar T3(a) T2 ins/2| maz (11 
mar marin = mar 


where the subscript ims refers to the case of perfect insulation at 
the that, for comparatively 
slowly a high thermal conductance at the 
On the other 


contact point 0 Fig. 5 shows 


varying heat inputs, 
joint is effective in reducing the maximum stress 


hand, for rapid pulses even a perfect thermal contact causes 
ilmost no stress alleviation 
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* Tables of the function have been calculated by Miller and Gordon, J 


Phys. Chem., Vol. 35, p. 2785, 1931 


because of 


! 


| 


Time variation of the temperature on each side of the joint for a pulse heat input, 4;=h 


the joint were perfectly insulated 
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Deflection of a Circular Section Cantilever 
Beam With All Fibers of Root Section at 
Elastic Limit Stress 


B. Saelman 


Design Weight Engineer, Lockheed Aircraft Corporat 
Burbank, Calif 
October 1, 195¢ 
SYMBOLS 
elastic limit stress, Ibs./sq.in 
deflection corresponding to a purely elastic stress distribut 


a vertical end load, Ibs 
section moment of inertia, in 
R outside radius of tube, in 
s of tube, in 
K section factor 


inside radiu 


length of beam measured from P, for which the stres tribut 


is purely elastic, in 
I Young's modulus, Ibs./sq.in 
distance from neutral axis te a point on section. in 


distance from neutral axis to fiber defining limit of elastic stre 


distribution, in 


p radius from center of tube to element of area on section, in 
1A element of area on section, sq.in 

VW bending moment at section, in.-lbs 

/ depth of rectangular section, in 


displacement of beam extremity as a result of deflection of tange 
at point A, in 
vertical displacement of point A, in 


D 
INTRODUCTION 


i REFERENCE | a comparison is made between the maximun 
deflections of cantilever beams of constant rectangular se 
tion and with concentrated end load when the root section has : 
purely elastic stress distribution and when all fibers of the root 
section are at the elastic limit stress; it is shown that the cor 
responding deflections are 180,/?/27Eh and 40¢,/?/27Eh, respec 
tively. These values indicate that the deflection corresponding 
to the plastic distribution is about 2.2 times as great as that for 
the elastic stress distribution 

Five methods for determining deflections of beams stressed 1 
the plastic range are discussed in reference 2. These are numeri 
cal integration of the moment-curvature diagram, mathematica 


integration of the idealized moment-curvature diagram, curva 





ture-area method, simple plastic theory, and plastic h 


method. This paper derives deflection formulas for flexure 1 


the plastic range for I-sections. The method used in the present 


paper is similar to that described on page 252s of this reference 
However, because of the formidable mathematics involved in this 


particular case, the strain-hardening effects as indicated wer 
neglected. 

In reference area-moment method fot 
range, the 
Deflections at 
loads are given with the test results obtained for rectangular 


rhe deflection at ultimat 


3 Meyer presents an 
the 
being made by graphical integration 


computing deflections in plastic calculations 


series ol 


square, and two tubular sections 


load is not given 
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TUBE SECTION 
(f) 


Region of beam from A to D has a purely elastic stress distribution. Region from A to B has a parti 
ind is fully plastic at B. »y, = f(x) defines the limit of the purely elastic stress distribution 


(@) be | 


il plastic distribution 


Seeley and Smith, in reference 4, derive a general formula (Eq $4) The stress-strain diagram can be idealized as in Fig. I(a 


645, page 547) for the aeflection in the plastic range in terms of 
the shape factor of the section. However, the equation is based strain 


and the elastic limit strain is small compared to the ultimate 


1a strain distribution (Eq. 640, page 546) which is valid only 
for rectangular sections. Hence, the general formula gives DISCUSSION 


exact” results, under the assumptions of the theory, for only Referring to Fig. 1(d), we have 
rectangular sections and approximate results for other sections ~e © R 
Reference 5 discusses methods for calculating the stresses be U= P(c+xu) =2 | o- (y?/y, dA 2 | o.ydA 
ond the elastic limit in rectangular and I-sections. However, ait by 
leflections are not discussed Ka.]/R 2 / i oC, a vd A I 
In view of the fact that references 1 and 2 do not discuss round see 
tubes in detail, that reference 3 does not present general results for Now, for thin tubes, 
round tubes nor the deflection at collapse moment for the two 7) "RO 
ises of round tubes discussed, and that reference 4 does not f, yd A I = p sin @ p dpda = R 
liseuss the accuracy of the authors’ formula for nonrectangular 
sections, a general formula for the deflection of thin round , cos 7 
tubes based on the simple plastic theory should be of interest rs "Ro — } 
rhe method used in the development of the formula for the a yA = | f, Yo Soo . , 
lastic bending deflection of a circular constant section canti ‘ = on 
ever beam with a concentrated end loaa is the same as that Ose ees 
presented by Van Den Broek in reference 1 and applied to a Substituting Eqs. (2) and (3) into Eq. (1), solving for ind 
rectangular section. The following assumptions are made: 1oting that c = o,//PR, 
1) The longitudinal strain varies linearly across the depth , , 
tthe beam - Ge to( Kk i 
2) The stress-strain relationship in bending is identical to . ' PR 3P 
that in simple tension tests 
3) The stress-strain relationship is the same in tension as in [ 40. R Zé oe R' ] cue a K i ; \ 
mpressior L 3P 2RP J 2RP in 9 
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then 
CO, ea dx CG, "0 o(R ee 
6. lim : = lim : 
gave C Ja: a—>ok = 2R2P 
to,( R8 r3) a(R r4 (= 0 6 cos ’) 2 
+ ire dé (0 
3PR 3R?P sin? 9 
2 o,°1 3 CO, I 
r= = Cae r r= (R4 ra 
3 EP?R? 8 EP?R3 
Now 
} J zig oep sin a@ dil 
0 
enr/2 LR . 
4 a J o. p sill a pdpda = (4/3)e,( R3 r = M (6 
where MR/I = Ka, (7 
Combining Eqs. (6) and (7), 
K = [4(R3 r3)R/3I] = (16/387) [(R r3)R/(R4 r4)| = 
(4/37) [(3n 6n? + 4n)/(n3 3n?2 + 4n 2 (8) 
I = (w/4) (R4 r4 (Ya 


Using the values of A and / from Eqs. (8) and (9a) in Eq. (5), 


5, = (3/32) (0,3/P2R3E) (R* — r4)?k 
(3 47/32) (o.3 P?R3E) (R4 yij2 nw 
(O.97K 0.92) [o,3(R4 r4)2/P2R3F (9 


From reference 1 and Eq. (6), 


; oe [°*'x dx a o£? BF BG 
6; = lim — = lim : : 
xa E PP a—oekE , a R R sin? 4 
BG 9 cos 9 F? FG cos 9 FG 9 cos? 9 FG 06 
cos 6 4 
R sin’ 6 R sin? 9 R sin® 6 R sin 6 
G? 9 G? & cos 6 i 
ey: : dd = (0.71 1L.19K O.61K?) > 
R sin’ 9 R sint ? 
o.*( R4 r4)2 
2 (10 
P2ER3 
where 
o, 1 to,( Rk? r oe/ 
B = (K 1) = 
RP oP RP 
m( R4 r4jo, 
tRP 
s to,( R® r3) o(R* — r4) o-l ; 2 
F = = K 
or 2RP PR T 
T 4 oO, 
kK (R4 r+) 
} T PR 
o-( ie r4 2 «J 
G = = 
2RP x RP 
We have eo, = PcR/I 11 
6. = Pc?/3EI] 12 


Combining Eqs. (9), (11), and (12), 


be = (7/48) (R4 r*)? (e,*/ER*P* 13 
Summing Eqs. (9), (10), and (13), 
6 =6. +6 + 6; = (0.61K? 0.22K 0.01) X 
[o,*(R4 r4)?/P2ER | (14 


For the root section, 


M = Pl = Ko,I/R = (#/4) (Ko./R) (R* r4 (15 
Solving for P in Eq. (15) and substituting into Eq. (14), 
6 = [(K? 0.36K 0.02)/K?] (¢,l?/ER) (16) 


For a thin tube, K = 1.27 and 


! 


i‘CAL SLIEANACES 


6 = V.70( o,/* ER 


Phe maximum deflection of the beam with a purel\ 


distribution is 


a1?/3ER 
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An Efficient Tripping Device* 


Francis R. Hama 

Assistant Research Professor, Institute for Fluid Dynamics and 
Applied Mathematics, University of Maryland, College Park, Md, 

October 15, 1956 


T’ HAS BEEN REVEALED in our recent investigations on transi 
tion from laminar to turbulent flow that the formation an 
development of a three-dimensional vortex loop are the esser 
tial features preceding the creation of a turbulent spot Asa 
application of this confirmation, an efficient turbulence-stimula 
tion or “‘tripping’’ device, which is commonly required for mode 
tests in a wind tunnel as well as in a towing tank, is proposed ir 
this note. ‘‘Efficient’’ here means that the device should provid 
a consistent and fixed stimulation with the least parasitic drag 
caused by the device 

Two-dimensional thin wire, a row of pins or spheres, or sand 
paper is conventionally used as an internal stimulation device 
whereas an external stimulation such as injection of air or water 
through a row of small holes is also employed in some applica 
tions. Any stimulation device, internal or external, is considered 
to produce vortices which eventually break down into randon 
motion. In order to shed vortices, a two-dimensional trip is far 
better than a three-dimensional one because the peak velocit) 
around a two-dimensional body is larger than that around one 
of comparable three dimensionality. Behind a two-dimensional 
trip, however, two-dimensional vortices are formed which might 
develop into three-dimensional configuration and eventual 
become turbulent, the course of the development depending 
upon Reynolds Number. This is the reason why the transitior 
approaches a trip wire rather slowly as the velocity increases.’ 
A three-dimensional tripping device, on the other hand, is mor 
direct in producing three-dimensional vortex loop, which is con 
firmed to be directly responsible for transition, as soon as the 
vortex shedding occurs, skipping over the two- to three-dimer 
sion process. A row of pins or spheres must be therefore superior 
in order to fix the transition point—i.e., a consistent stimulatiot 
is demonstrated in the comparison experiment performed | 
Klebanoff, Schubauer, and Tidstrom The superiority of three 
dimensional stimulation is also widely recognized among thos 
of the towing-tank circle. * 

Proposed here is a row of thin triangular patches pasted on 
flat plate as shown in Fig. 1. Streamlines near the triangle wil 
be as sketched in the Figure and will promote a spiral motion 10 
herent to a three-dimensional vortex loop. Furthermore, duet 
converging streamlines near the bottom corner of two neigh 
boring triangles, the fluid velocity adjacent to the surface will be 


increased, resulting in an efficient vortex shedding even bette! 


Air Force through the Of 
18(600)-1014 


* This research was supported by the U.S 


of Scientific Research, ARDC, under Contract AF 
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Comparative flow pattern behind triangular-patch and 
1.0 mm., x. = 30 cm., U; = 18.6 
420, k/56* = 0.46, Urk/v = 50) 


Fic. 2 
two-dimensional trips (k = 
em./sec., U;6*/py = 
produces 


two-dimensional trip. This arrangement 


three-dimensional vortex loops directly with more efficiency in 


than a 


shedding and therefore appears to provide a simple yet better 
way of tripping laminar boundary layers than any other known 
stimulation device. 

A comparison experiment is performed on a flat plate submerged 
in a towing tank, placing a two-dimensional trip and a row of 
triangular patches side by side at 30 cm. from the leading edge 
Each triangle is a regular triangle and is made of five layers of 
black insulating tape so as to form the trip thickness of approxi 
mately 1.0 mm., whereas the two-dimensional side is the same 


4 in width. Flow 


material and thickness of approximately 2 mm 
pattern is visualized by a dye issued two-dimensionally through 
i narrow slit cut across the plate at 20 cm. from the leading edge 
It is distinctly shown in Fig. 2 that the two-dimensional trip 
barely produces a weak two-dimensional disturbance which dis- 
ippears almost immediately, while the new device is already 
stimulating the flow almost right behind the trip, indicating 
the marked effectiveness of the device 

The effectiveness will be better illustrated by the Reynolds 


Number U;.k/» when the transition reaches a trip. U; is the ve 


FORUM 937 


locity at the trip thickness & and at the location of the trip calcu 
lated under the assumption that the trip does not disturb the 
kinematic vis 


smooth flat-plate bourdary laver, and »y 1s the 


It is reported that this Reynolds Number amounts t 


cosity 
approximately 570 for a row of spheres and at least 200 for a trip 
For the new device, on the other hand, the corresponding 
Needless to mention, this 


wire 
Reynolds Number is estimated at 45 
pronounced difference in figures implies that a lesser roughness 
thickness is sufficient to trip the flow with this device than with 
the other devices under equivalent flow conditions so that a less 


parasitic drag due to a trip can be expected. Moreover, a con 
sistent transition is assured because of the three-dimensional 


nature of disturbance 
Further systematic investigations searching for the best com 
now in progress 


bination of the dimensions of the triangles are 


ind will be rey wrted later 


REFERENCES 








Hama, F. R Long, J. D., and Hegarty, J. 4 ae ’ 
Laminar to Turbulent Flow, J. Appl. Phy in press 
rani, Itiro, and Hama, | R Some Experimen f Effe 
Single Roughness Element Boundary-Layer Transition, Readers’ Forum 
Journal of the Aeronautical Sciences, V« 0, N i, pp. 289, 290, Apr 
1953 
Klebanoff, P S., Schubauer, G. B., and Tidstrom, k. D., Me ren 
Effect of Two-Dimensional and Three-Dime ul R } I ? 
Boundary-Layer Transit Readers’ Forum, Journal of the Aeronautic 
Sciences, Vol. 22, No. 11, pp. 803, 804, November, 1955 
‘ Hughes, G., and Allan, J. I Tur Si S Ff 
frans. Soc. Nav. Arch. Mar. Eng., Vol. 59, pp. 281-314 
+ 


Limitations of Circulation Lift 


Heinrich B. Helmbold 

Research Scientist, Fairchild Aircraft Division, Fairc 
Airplane Corporation, Hagerstown, Md 

November 2, 1956 


w™ A SPANWISE elliptic distribution of circulation, IP 
MV 1 (2y/b)*, the trailing-vortex system far down 
stream of the wing will move like a rigid flat strip with the proper 
velocity wo = Io/b normally to itself, provided that its cros 
section is hypothetically endowed with the ability to resist thi 
deforming stress of the pressure field. The transversal momen 
tum carried with the vortex sheet increases per unit length b 
dJ»/dxo = p(xb?/4)wo, and the kinetic energy of the absolut: 


on.2 » 
(2297/2 


motion increases per unit length by dEo/dxo = p(2b?/4 
If now the cross-sectional rigidity is instantaneously destroyed 
ind the sheet is left free to yield to the deforming stresses of the 
pressure field, the sheet will start rolling itself in from its side 
edges and finally attain a stable configuration in the form of two 
symmetrical vortex cores of opposite and equal circulation Ty 
The distance between their centers of gravity is / th/4, and 
their proper velocity is 

w’ = T)/2al = 


it least approximately Hence wo = (77/2)u Since no outer 
forces are acting on the vortex system during the deformation, 
the transversal momentum and the total energy per unit length 
remain unchanged by the deformation; d/J;/dx = dJ »/dx»9 and 
dE/dx = dEo/dxp. (With the deformed vortex system the total 
energy consists of kinetic energy and work done by pressure 
forces: with the undeformed vortex sheet the kinetic energy 1s the 
total energy because then there are no perturbation velocities 
along the vortex sheet. ) 

The completely rolled-in vortex system is inclined by the angle 
5 toward the undisturbed relative velocity V. This angle ts 
The length of the final vortex 
The com 


determined by sin 6 = w’/V 
system increases per unit time by dx/dt = V cos 6 
ponents of the aerodynamic force on the wing are the lift 








238 JOURNAL OF T TRONAT 














Fic. 1. 


L = V+(dJ,/dx) — Dj; tgé 2 


cf. reference 1, Eq. (68)| and the induced drag 


D; = (dE/dx) cos 6 ° 
since D;V = dE/dt = (dE/dx) V 


According to preceding statements, 


L= 


cos 6 


) 


p (rb?/4) Vu p( 7b?/4) (wo?/2) sin 6 \ 
= p (7°/8) b* V? sin 6[1 (w?/4) sin? 6 
D = p wh? /4) (wo?/2) cos 6 . 
) 
p (2°/32)b? V2 sin? 6 cos 6 
w, after division by pV2b?/2, 
C,/A = $) {1 r?/4) sin? 6 6 
Cp;/A = (2*/16) sin? 6 cos 6 7 
is illustrated by Fig. 1 Here A = b?/S denotes the aspect ratio 


Phe reduced lift coefficient C,/A attains its maximum value, 


Cc; A r/(3V— 3 1.90, for sin 6 = 2/(47YV3 0.368 
with a corresponding value of the reduced drag coefficient 
Cp;/A = 2/12) VY x? (4/3) = 2.40 
REFERENCES 
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Flow Disturbances Induced Near a Slightly 
Wavy Contact Surface, or Flame Front, 
Traversed by a Shock Wave* 


G. H. Markstein 

Principal Physicist, Aerodynamic Research Department 
Cornell Aeronautical Laboratory, Inc., Buffalo 

October 18, 1956 


R“ NT OBSERVATIONS of the effects of shock waves on flame 
structure’ ? were originally interpreted in terms of Taylor 
instability However, it was later realized? that the concept of 
Taylor instability*® had to be modified in order to account for the 
effect of 


with passage of a shock wave 


sudden acceleration of very short duration associated 


In reference 2, this modification 
was only qualitatively discussed. Although the following anal 
vsis is based on fairly drastic simplifications, it should yield at 
least order-of-magnitude estimates of the effects to be expected 
Consider a plane contact surface separating gases of density 
and initially at rest 


p.\ for —2 <x < Oand p for0 <x < + 


* This research was conducted under the auspices of Project SQUID 
jointly sponsored by the Office of Naval Research 
Office of Scientific 


Department of the Navy 


Research, Department of the Air Force, and Office of 


Ordnance Research, Department of the Army 
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When this surface is traversed by a shock wave propagating 
the positive x direction, the same constant velocity iS induce 
} 


by the transmitted wave for x > 0 and by the combined act 


of the incident and reflected waves for x < 0. If, 


contact surface initially deviates somewhat from strictly pla 


shape, the reflected and transmitted waves will generat 


tional flow disturbances in its vicinity Instead of treating t 


ictual problem of reflection and transmission of a shock way; 


i slightly distorted contact surface, the following analysis ; 


places it by the much simpler one of computing the effect 


constant acceleration a of a very short duration 7 and of y 


magnitude that the relation 


ni ] 


is satisfied rhis simplification is equivalent to neglecting , 
pressibility effects or reg irding the speed of sound as infinite 


that the flow disturbances, instead of being generated by th 


flected and transmitted waves, appear simultaneously throu 


out the flow field 


s 


Taylor’s analysis* can be immediately applied to the 


fied problem. For a contact surface initially (¢ < 0 


of sinusoidal shape 


where & is the wave number and A the amplitude, one obt 


for the time interval 0 < ¢ r the expressions 


Ey, f fo(¥) cosh at 
for the shape, and 
u(x, y,t) = F w(x, y, t) = atoll ye ‘sinh af 
for the flow disturbances in the x and y directions, where t 


upper signs apply for x < 0 and the lower ones for x > 0 and wher 


_ on fone 
a= Ff p2 p (p tT p Raj 


is satisfied, the shape of the contact surface will remain practical 


unchanged during the acceleration, while a flow field 


u(x, vy, rT) = Fiv(x, y, 7) = a®rio(y)e*** 7 
will be established Substituting Eqs. (1), (2), and (5) into Eq 
7), one obtains finally 
u(x, y) = F iv(x, y) = UkRA|(p p po + p,)le*$F*T" (8 


for the disturbance fiow field present immediately after thea 
celeration has taken place. The result can be applied to t 
passage of a shock wave through a contact surface since it depen 
the 
ind unknown acceleration a and time interval r 


only on known velocity change Ul’, not separately on t 
nctitious 
can also be shown that it applies equally well to a flame front 
since, for large acceleration of short duration, the contributions 
from other terms arising in the analysis of flame-front stabilit 
become in comparison negligibly small 

the flow fiel 


This flow wil 


As no further acceleration takes place for ¢ 2 rf, 
u, v will persist for some time with little change 
cause distortions of the contact surface that grow approximatel 


linearly with Since the sign of « depends on that ‘ 


J 


time 
(p2 p,), the distortion will be in the same direction as th 
initial one if the shock wave traverses the interface from the 
thinner towards the denser medium, and vice versa The result 
is thus in agreement with the qualitative argument and the e& 


? suff 


perimental observations presented in reference 2. For 
ciently large initial distortion measured by 


number kA, which also represents the slope (d&/dy , and for 





1, the disturbance velocities may be 
Measure 


large density difference ps 
of the same order of magnitude as the velocity Ul’ 
ments of the velocity of the 
a flame front with a shock wave? showed indeed that it exceedet 
the flow velocity computed on the basis of one-dimensional theor 


by a factor in the range of 1.5 to 1.8 


the dimensionless 


“spike’’ formed during interaction 0 


| 
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Miniaturized Heat Meter for Steady-State 
Aerodynamic Heat-Transfer Measurements* 


Hartwig, C. A. Bartsch, and H. McDonald 
Aeronautical Laboratory, California Institute of 
xy, Pasadena, Calif 

195¢ 





TEADY-STATE aerodynamic heat-transfer measurements have 
S ilways presented severe difficulties. When pairs of thermo 
uples are used to measure the temperature gradient in a solid 
the differential millivolt output is usually quite low, unless 
But 


he heat flux is generally not normal to the surface because of 


ody, 


the two thermocouples are spaced reasonably far apart. 


ixial 
required to obtain the normal heat flux are often as large as 50 


temperature gradients in the material, and the corrections 


to 100 per cent of the measured values 


In an attempt to avoid these difficulties, a device has been 


mstructed which allows the distance between the thermocouple 


uir to be greatly reduced while the required millivolt output for 


given heat flow is maintained. The heat transducer or heat 
neter is really a thermopile on a miniature scale (Figs. 1 and 2 
[he meters were wound on a glass core 0.007 in. by 1/16 in. by 1/8 


Fifty turns of O.001-in. constantan wire spaced 0.001 in 


ipart were wound on the core. Then the long edge of the meter 
vas immersed in a silver-plating solution so that one half of each 
Thus an effective silver-constantan junction 


0D 


p was plated 


vas formed on each loop at the center of the two flat sides of the 
To prevent shifting or damage to the wires, a thin coat 
Adweld for 
) or Sauerisen for high tempera 


meter 
f high-temperature cement was applied, such as 
medium temperatures (500°F 
tures (1,000°F The 
).010 in. thick. For areas of high heat 
stagnation point, only 25 turns were used so that the meter meas 
by 1/16 in. by 0.010 in 


For our particular installations, a 


about 
the 


finished meter then measured 


transfer such as 


ired only 1/16 in 


hollow metal model was 


vated with 0.020 in. of porcelain into which a shallow indenta 


tion was ground. The meter was cemented into this indentation, 


the surface was smoothed over, and external connections were 


nade by means of silver paint to wires at the rear of the model 


* The work discussed in this paper was carried out under the sponsor 


f the Office 
Army 


ship and with the financial support Chief of Ordnance, and 





the Office of Ordnance Research, U.S 





1 (left) 


Fifty-turn heat meter 
Heat meter compared with ordinary pin head 


Fic. 
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The sensitivity of each meter is affected slightly by manufac- 


turing differences and even more by the method of installation 


and the thermal conductivity of the material in which it is in 


stalled. Therefore, it is necessary to calibrate the meters after 
installation In our case, the calibration was carried out in an 
oven capable of radiating 10,000 watts/ft.? continuously. The 


heater element wires were spaced sufficiently close together so 


that the model surface ‘‘saw’’ a nearly uniform source of radiant 


area where 
of a 


rest of the oven 


heat In order to prevent heat losses through the 


combination null 
The 


i stainless-steel reflector sheet surrounded by thre¢ 


the model support entered the oven, a 


heat meter and a bucking heater was used 


was covered by 


layers of asbestos sheet, outside of which were 3 in. of rock wool 
Commercial heat meters were used to measure the actual losses 


through tLe outer surface of the oven These losses never ex 


ceeded 5 per cent of the total input. Since the model was water 
cooled, its surface temperature was below the boiling point of 
water. Therefore, to prevent condensation, the atmosphere in 


the oven was replaced by dry nitrogen \ typical calibration 
a 50-turn meter is shown in Fig. 3 These meters have 
been used for heat-transfer the GALCIT 


7.8, and the accuracy and repeatability so far 


curve for 


studies in hyperson 


tunnel at JJ = 
obtained are very encouraging 


A Technique for Computing Natural 
Frequencies of Trapezoidal Vibrators 
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SYMBOLS 
Cartesian rectangular coordinate 
time 
Young's modulus 
b/ 12 moment of inertia of a cross section witl 
the neutral axis 
density 
1 Abell \ cros 
X 
thicknes 
1 X 
b breadths at 
} h 


sectional area 


breadth 
built-in and free ends 
length 
} zcos ft 
mode 
angular velocity 


l2pf I 
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v/vo Versus ¢ 


vie 
certain independent variable 
certain functions of £ 


certain constant 


DISCUSSION 
r | NHE DIFFERENTIAL EQUATION which governs the flexural 
vibration of a bar with variable cross section is 
(0?/O0X?) [EI(0? ¥/0X?)] + pA(o?Y/o?? 

which reduces to 
l (c/1)X](042/0X*) — 2(c/Ll) (0%2/0X*) 

(v4//4) [1 (c/)Xis = 0 (2 
in the case of the trapezoidal vibrator here considered. The 
boundary conditions at the free and built-in ends are 


(z)x=0 = (dz/dX)x =0 = (d*2/dX?)x =1 = (d*e/dX*)x =i; = 0 


With the introduction of a new independent variable — and the 


expansion of all variables in power series of c—i.e., » = w(1 4 
ch), y=(1 tere = £ +0(€), 2 = u(é) + cv(€) to the first 
order of c—Eq. (2) assumes the following form: 

(1 — c&) [(u'""’ + cv’) (1 + ey’) 


,_ erry 
} 


6cu'’’n’’ — 4Acu’’n 2c(u’’’ + cv’"’) 

wi(1 — cE) (ua + cv) (1 + 5en’) = 0 (83) 
where the primes indicate differentiation with respect to — On 
the other hand, the boundary conditions become 


(s\g=0 = (ds dé )¢ 0 = (d2z d&*); 1 


Cu’(0) + Du’’(0 


Bu'(0O) + Cu’’"(0) 
Bu(l) + Cu’(1) 


Cu(1) + Du'(1) 


Au(l) + Bu'(1 


We are here concerned only with the determination of (1), 
which, when solved, is found to be 
m1) = (1/04) [u’(1)0’’’"(O)/u(1)u’"(0)] = [cot (v/2)/vr0]? 
Thus the final result is, to the first approximation, 
9 


. 2) 
wi 1 + cleot (v/2)/v]?} 


numerical values of which are shown in Table 1 and Fig. 1. For 
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TABLE | 


v/vo Versus ¢ 


Ono's 
Fundamental - 

The present 
First overtone 
Second overtone 
Third overtone 
Fourth overtone 


From Eq. (3) we obtain, as the 


“ ‘4 = iF (Wu Jr=O 


which represent the flexural vibration of a uniform bar. The 
solution of the vibration is very well known but will be repro- 
duced below for reference: 
u = (const.) [COS mé cosh wyé COL wMy( SIN YoF sinh wok ) 
cosh vw = —1/cos », sinh vg = tan pv 
= 1.8751, 4.6941, 7.8548, 10.9955, 14.1372 


The first-order equations are 


votv “un 


a 


$ 


n( 0) = 


When solving these equations we choose 7 to satisfy 


i.e., in the same form as u Then 7 must satisfy 


d*/d&*) (u'n) votu'n 2u!"’ = 0 
HO) = "(1 = 7"(1 
If we put u'n H 
we obtain a differential equation and boundary conditions for H 


as follows: 


la wy tH 


(1/2)tu + . ! ua’? + Du'”’ 


where A, B, C, D are arbitrary constants. These constants are 


determined by 


u!’""(1)n(1)/» 
u'(1)/m! 

u’(1)n(1) u(1) 
the sake of comparison, Ono’s* results have been reproduced in the 
same table and figure. They show good agreement concerning 
their behavior near the origin, and we can conclude that the 
first approximation is valid up until c = 0.6 if the error of 5 pet 

cent is allowed. 

* Ono, Akimasa, Lateral Vibrations of Tapered Bars, Journal of the 
Society of Mechanical Engineers, Japan, Vol. 28, No. 99, p. 429, July, 1925 
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